
///' 


£> 


UAH Propulsion Research Center W 


Determination Of Uncertainties For The New SSME Model 


Final Technical Report 


NASA/Marshall Space Flight Center 
NASA Contract Number NAS8-38609 D.0. 140 


Prepared by 



Kendall K. Brown 


Submitted to: 

John P. Butas, EP-14 
Marshall Space Flight Center, AL 35812 


Dr. Hugh W. Coleman 
Principal Investigator 


Dr. Clark W. Hawk 
Center Director 


Propulsion Research Center 
The University of Alabama in Huntsville 


May 16, 1996 




Final Technical Report 
1 May 1995 to 16 May 1996 

Determination of Uncertainties for the New SSME Model 


National Aeronautics and Space Administration 
Marshall Space Flight Center 
Contract Number NAS8-38609 
Delivery Order #140 

by 

Kendall K. Brown and Hugh W. Coleman 
Propulsion Research Center 
University of Alabama in Huntsville 




Table of Contents 


Section Tonic Page 

1.0 Introduction 3 

1.1 Statement of Work 4 

1.2 Relationship of This Contract to Overall Effort 6 

1.3 Organization of This Report 6 

2.0 Test Data Uncertainty Considerations 8 

2.1 Development of the Baseline Data Set (BDS) 8 

2.2 Data Reduction Data Set 9 

2.3 Precision Uncertainty Estimates 9 

2.4 General Discussion of Conceptual Bias Uncertainties 10 

2.5 Discussion of Temperature Uncertainties 10 

2.6 Discussion of Pressure Uncertainties 11 

2.7 Discussion of Venturi Systematic Uncertainties 12 

2.8 Discussion of Valve Position Uncertainties 12 

2.9 Discussion of Speed Measurement Uncertainties 13 

3.0 Modeling Uncertainty Considerations 14 

3.01 Modeling Assumptions and Approximations 14 

3.02 Uncertainties from Previous Information 15 

3. 1 Discussion Of Physical Property Data Routine 16 

Uncertainties 

3.2 Discussion Of Combustion Routines 17 

3.3 Discussion Of Hardware Characteristic Uncertainties 18 

3.3.1 Discussion Of Turbine Map Uncertainties 19 

3.3. 1.1 HPFT Map Uncertainty 19 

3.3. 1.2 Other Turbine Performance Map Uncertainties 21 

3.3.2 Discussion Of Turbopump Pump Map Uncertainties 2 1 

3.3.3 Discussion Of Duct And Line Resistances 22 

3.3.4 Discussion of Heat Transfer Calculations 23 

3.4 Propagation Into Balance Relation Uncertainty 24 


1 




3.4.1 

Duct and Valve Flow Modules 

24 

3.4.2 

Turbine Modules 

25 

3.4.3 

Pump Modules 

26 

3.5 

Discussion Of Balance Relation Uncertainties 

26 

4.0 

Summary 

28 

4.1 

Uncertainty Estimate Improvements 

28 

A-l 

Appendix 1 - Uncertainty Analysis Methodology 

Al.l 

A-2 

Appendix 2 - Baseline Data Set And Data 
Reduction Data Set And Uncertainty Estimates 

A2.1 

A-3 

Appendix 3 - Uncertainty In Linear Regressions 

A3.1 

A-4 

Appendix 4 - Monte Carlo Simulation of 
HPFT Efficiency Map 

A4.1 

A-5 

Appendix 5 - Balance Relation Uncertainties 

A5.1 

A-6 

Appendix 6 - AIAA Technical Paper 95-3073 “Enhancing 
Rocket Engine Test Analysis And Performance Models 
With The Incorporation Of Uncertainties.” 

A6.1 

A-l 

Misc. Background Information 

A7.1 


2 




1.0 Introduction 

This report will discuss the methodology used and the development of 
uncertainties for use with the New SSME Model. The New SSME Model was 
developed by Dr. L. Michael Santi of Christian Brothers University through a 
separate contract effort. The new model requires information about the 
experimental uncertainties from the Technology Test Bed (TTB) test 
measurements and about the uncertainties within the component models of 
the new model. 

This effort was a follow-on and continuation of an engineering and 
research effort begun under delivery order 106, which ter min ated 30 April 
1995 1 . The primary products of that effort were an assessment of the 
experimental uncertainty in the determination of the venturi flowrates and 
the initial development of a methodology to assess the uncertainty in linear 
regressions. 

The primary products of this effort include information about the 
experimental uncertainty in additional TTB measurements (i.e. 
temperatures, pressures, turbopump rotational speeds, and valve positions), 
additional assessment of conceptual bias uncertainties, development of 
methodology to assess the uncertainty in linear regressions, and assessment 
of the uncertainty in the model hardware characteristics. 

A three month period at the beginning of this effort, June, July, and 
August was excluded from the contract performance period. This was done to 
allow Mr. Kendall Brown, Graduate Research Assistant working on this 
contract, to participate in the Air Force Office of Scientific Research 
Graduate Summer Research Program at the U. S. Air Force’s Arnold 
Engineering Development Center (AEDC). The purpose of his participation 
was to work with engineers at AEDC who were also interested in developing 
methods to incorporate uncertainties in the modeling and testing programs of 
advanced aerospace systems. While at AEDC, Mr. Brown made substantial 
progress in developing a new methodology to assess the uncertainty in linear 
regressions (curvefits) 2 . The products of the AFOSR Summer Research 
Program were immediately applicable to this effort and was a great 
leveraging of resources. 


1 Coleman, Hugh W., and Brown, Kendall K, "Impact of Uncertainty on Modeling and Testing,” 
Propulsion Research Center Report # 95-001, Final Report on NASA-Marshall Space Flight Center 
Contract NAS8-38609 D.O. #106, 10 January 1994 to 30 April 1995. 

2 Brown, Kendall K., “A Methodology for Assessing Experimental Uncertainties in Curvefits with 
Application to Compressor Map Characteristics,” Air Force Office of Scientific Research graduate 
s umm er Research Program Final Report, Arnold Engineering Development Center, Arnold AFB, TN, 
Aug 1995. 
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1.1 Statement of Work 

The Statement of Work for this contractual effort is provided. The 31 
March 1996 delivery order termination date was extended to 16 May 1996 
under a 45 day, no-cost extension 

The work identified under this statement of work is a continuation of work 
performed under UAH Contract NAS8-38609, Delivery Order No. 106. The 
primary focus of the work performed under Delivery Order No. 106 was to 
establish methods for generating uncertainty estimates for selected 
Technology Test Bed (TTB) experimental flowrate measurements and to 
assess the impact of these uncertainties on developing a new test integration 
strategy. The goal of this new effort is to generate or finalize uncertainty 
estimates for TTB pressure and temperature measurements and to evaluate 
these estimates within the new Space Shuttle Main Engine (SSME) steady- 
state performance model. The incorporation of these uncertainty estimates 
will enhance the capability to better support current performance analysis 
requirements such as, assessing vehicle/engine feed system interface flow 
characteristics, engine hardware design changes, evaluating engine 
hardware performance, predicting engine hardware operation, and 
supporting failure investigations. 

Rocket engine performance models consist of mathematical relationships 
which model physical processes within the engine. Experimental data is use 
to anchor these physical relationships. Performance models are used to 
derive engine component hardware characteristics from experimental data. 

Both experimental data and the physical relationships within a model 
contain various sources of errors. These sources of errors include calibration 
errors, signal processing and localized effects, uncertainties in both physical 
approximations and fluid property data. These sources of errors are 
neglected within the current SSME steady-state performance model. 
Experimental data is treated as absolute and computational predictions are 
forced to agree with the data at instrumented locations, often at the expense 
of physical consistency. This situation impacts the integration of 
experimental data within the model, thus reducing the accuracy of the 
computed hardware characteristics. 

A new SSME steady-state performance prediction program is currently 
being developed to evaluate and utilize experimental data derived from 
recent TTB tests. A modified test data integration scheme has been 
developed which reconciles uncertain experimental data with uncertain 
physical relationships. This strategy systematically transforms uncertain 
experimental data into a physically self consistent set of data. This is 
accomplished by forcing the minimum adjustment required in engine fluid 
measurements to satisfy prescribed uncertainty constraints. This 
reconciliation scheme was reformulated during this last year based on the 
method developed in estimating the experimental uncertainties of the few 
TTB measurements. The overall success of the reconciliation strategy is a 
function of determining a reliable method to both generate and incorporate 
these required uncertainty estimates. Thus, the incorporation of these 
uncertainty estimates will enhance the use of experimental data to anchor 
the physical relationships which in turn, will improve the overall accuracy of 
the performance predictions. 
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The research required to implement these uncertainty analysis concepts will 
be conducted within the SSME engine 3001 test program which is currently 
being conducted on the TTB test facility. Engine 3001 provides a 
significantly larger number of experimental measurements as compared to 
standard SSME flight engines. This test program provides a unique 
opportunity to assess how accurate we can predict SSME hardware 
performance based on selected instrumentation. 

A two phase research effort will be established to support the modified test 
data reconciliation strategy, and ultimately will improve the use of 
experimental data in generating performance predictions. Phase I involves 
generating or finalizing uncertainty estimates for TTB pressure and 
temperature measurements that were not generated under UAH Contract 
NAS8-38609, Delivery Order No. 106. Phase II involves estimating the 
uncertainties associated with incorporating previous experimental data and 
fluid property data within the model. A detail description of the specific 
tasks to support these phases are described below: 

Phase I - Estimate uncertainties for TTB pressure and temperature 
measurement. 

Task 1) Perform statistical analysis on existing pressure and 

temperature TTB test data. 

Task 2) Identify all significant sources of errors within 

instrumentation system for these measurements 

Task 3) Estimate both precision and bias uncertainties for the TTB 

pressure and temperature measurements 

Task 4) Estimate uncertainties associated with localized effects due 

to measurement locations. 

Phase II - Estimate uncertainties for incorporating previous 
experimental data and fluid property data within the new SSME model. 

Task 1) Identify the sub-models within the new SSME model which 

incorporate previous experimental data and fluid property data. 

Task 2) Estimate uncertainties associated with the use of previous 

experimental data (duct and valve characteristics, for example). 

Task 3) Estimate uncertainties associated with fluid property data. 

The following represent the deliverable products expected for each of the two 
phases: 

Phase I - Estimate uncertainties for TTB pressure and temperature 
measurement. 

1) Interim report documenting uncertainty estimates as completed for each 
of the TTB measurements selected. 

Phase II - Estimate uncertainties for incorporating previous experimental 
data and fluid property data within the new SSME model. 

1) Final report documenting research conducted in developing the 
methods for estimating uncertainties in test measurements and in the 
model. 
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1.2 


Relationship of this contract to overall effort 


A technical paper for the 32nd AIAA/SAE/ASME/ASEE Joint 
Propulsion Conference was jointly authored by the researchers on this 
contract, Dr. L. Michael Santi, and the COTR, Mr. John P. Butas 3 . This 
technical paper presented the overall effort to the propulsion community and 
was well received. A copy of this technical paper is provided as Appendix 4. 

As described in the technical paper, the overall goal of this effort is to 
develop a new model of the Space Shuttle Main Engine which incorporates 
the uncertainties in the experimental test program and some of the 
uncertainties in the modeling process. A key feature of the new model is that 
it is a physical model, it satisfies the conservation of mass and energy, 
whereas the existing model does not. As the developer of the new model, Dr. 
Santi, of Christian Brothers University, Memphis, TN, was responsible for 
developing the numerical solution strategy and how to incorporate the 
experimental and modeling uncertainties. As these parallel research efforts 
proceeded, this effort of assessing of the various uncertainties evolved into 
supporting Dr. Santi’s modeling effort. Extensive communication with the 
COTR and Dr. Santi ensured proper information. 


1.3 Organization of this report 

This report is broken into two main sections, first a discussion of 
uncertainties in the experimental data, and secondly a discussion of the 
modeling uncertainties. 

In order to support the new model, the assessment of experimental 
uncertainties in almost all types of TTB SSME measurements was required, 
the notable exceptions were strain gage and accelerometer measurements. 
An initial assessment of temperature, pressure, and flowrate uncertainties 
was conducted under the previous contract. Further investigation into the 
experimental uncertainty in these measurements, as well as investigation of 
the turbopump speeds and the valve positions sensors, was conducted under 
this effort. Conceptual bias uncertainties are often sig nifi cant components of 
the experimental uncertainty, and discussion of efforts to assess these 
uncertainties is included. 

The primary model uncertainty sources being included in this effort 
are the uncertainty in the physical property data, the uncertainty in the 
hardware characteristics, and the propagation of these uncertainties into the 


3 Brown, Kendall K., Coleman, Hugh W., Santi, L. Michael, and Butas, John P., “Enhancing Rocket 
Engine Test Analysis and Performance with the Incorporation of Uncertainties,” 31st 
AIAA/ASME/SAE/ASEE Joint Propulsion Conference, AIAA paper 95-3073, San Diego, CA, July 1995. 
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uncertainty associated with the balance relations within the model. The 
assessment of the uncertainty in the hardware characteristics and their 
propagation is a new topic and methods to assess these parameters had to be 
developed. Many of the SSME components, particularly the turbopumps 
utilize maps, or curvefits to represent their performance. These maps are 
generated using other information, and usually by conducting a sub- 
component test program and scaling the results to the engine operating 
conditions. The available documentation of how each specific hardware 
characteristic was generated is very poor, and thus the information presented 
in this report represents the best information that could be obtained. 

In some cases the information needed to make a more accurate 
uncertainty estimate was not available, so this report will present the 
methodology used to assess the uncertainties, so that when the appropriate 
information becomes available the uncertainties can then be estimated and 
included in the model. 
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2.0 Test Data Uncertainty Considerations 

The assessment of the experimental uncertainties of the data obtained 
from the Technology Test Bed (TTB) SSME Engine 3001 test program was 
conducted using the methodology described in Appendix 1. As described in 
Appendix 1, the determination of uncertainty is comprised of two main 
components, bias (systematic) uncertainties and precision (random) 
uncertainties. The manner in which the uncertainty value is used dictates 
how the uncertainty is estimated, and conversely, the way in which the 
uncertainty is estimated dictates how it may be used. Thus, a proper 
understanding of how the uncertainties would be used in the new model was 
necessary in order to make appropriate uncertainty estimates. This section 
describes how the uncertainty estimates were determined for use with the 
new model. 


2.1 Development of the Baseline Data Set (BDS) 

A baseline data set (BDS) was developed in order to specifically 
identify the measurements to be used to validate the new model and solution 
methodology. This data needed to represent the best information available 
about the operation of the engine, and the associated uncertainty. The new 
model contained hardware characteristics representing a Phase II engine, 
thus the data used to build the BDS must be taken from TTB tests with 
Engine 3001 configurations similar to a Phase II engine. This limited the 
available tests to TTB021-TTB038 and TTB052-TTB065. This group of tests 
was further limited because of a series of tests with the Pratt and Whitney 
ATD HPFTP. Since the performance of the Pratt and Whitney turbopump is 
different than the Rocketdyne turbopump, and since the new model is using 
HPFTP maps representing the performance of the Rocketdyne turbopump, 
this data had to be excluded from the BDS. The test profiles from the 
remaining tests were examined to find tests which were conducted at the 
same engine operating conditions. The primary engine control conditions 
used in the test program are the power level, the fuel/oxidizer mixture ratio, 
the propellant inlet pressures, and the propellant repressurization flowrates. 
It was desired to get sets from four different power level settings, 100% rpl, 
104% rpl, and 109% rpl, and a low power level (90% rpl or lower). There was 
insufficient data available at a low power level to create a data set for use in 
the model. 

After the four tests were identified, specific 5-second time slices within 
each power level of each test were selected. These 5-second time slices were 
chosen in regions of the test profile where the engine was operating in 
essentially a steady-state, regions sufficiently far from the engine start 
transient or power level changes where those effects would not be included in 
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the data. The data from within these 5-second time slices were averaged to 
obtain a single test point representing the operation of the engine during 
that test and at that power level. The data point from each of the four tests 
were then averaged to obtain the data point for the baseline data set. The 
data in the BDS is the best representation of the operation of the SSME at 
the given operating conditions. 

The TTB measurements included in the BDS were chosen based upon 
a sensitivity study performed by Dr. Santi and the personal experience of the 
COTR in performing SSME analyses. The precision uncertainty estimates to 
associate with the values in the BDS must reflect the test-to-test variation. 


2.2 Data Reduction Data Set 

One of the main purposes of the New SSME Model is test data 
reduction, incorporating the test data into the model to find anomalous 
readings. As such, a second set of data needed to be generated to use in the 
validation of the new model. The same base set of measurements were used 
and additional measurements were included, these additional measurements 
could be used to increase the number of data points used in data reduction or 
for comparison with output from the new model. The precision estimates to 
associate with the values in this data set must include the test-to-test 
variation and the variation of the data within tests. 


2.3 Precision Uncertainty Estimates 

The appropriate precision limit to use with a given set of data must 
reflect how that data was generated, the usage of the data, and the usage of 
the uncertainty estimates. Since the baseline data set is a single point which 
represents the average of data points from four separate tests, the 
appropriate precision limit is based upon the standard deviation of the four 
points. This precision limit represents the test-to-test variation of the engine 
operation. As discussed in Appendix 1, when using less than approximately 
ten points in a standard deviation calculation the large-sample 
approximation should not be used. The sample standard deviation calculated 
with the BDS is shown in Table 1 in Appendix 2. 

The precision uncertainty to use when using the new model for test 
data reduction is slightly different. While it still must include the test-to-test 
precision variation it must also include the within test variation for the 1- 
second averaged data. Hence the sample standard deviation is determined 
based upon five data points in the four tests, for a total of twenty data points. 
Since twenty data points are being used in the sample standard deviation the 
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large-sample approximation can be used. Table 2, in Appendix 2, shows the 
measurements included in the data reduction data set, the average value, 
and the uncertainties at 104% rated power level (rpl). It shows that the 
precision uncertainty is approximately equal to that determined for the BDS, 
indicating that the test-to-test variation dominates the within -test variation. 

A pooled sample standard deviation for the data reduction data set 
was also calculated and verified that the test-to-test variation is much 
greater than the within-test variation. Using the pooled standard deviation 
is not appropriate since it does not include the test-to-test variation. 


2.4 General Discussion of Conceptual Bias Uncertainties 

Potentially sig nifi cant bias uncertainty sources to consider in the TTB 
SSME measurements are conceptual bias uncertainties. The conceptual bias 
uncertainties in the temperature and pressure measurements are 
particularly important in this effort because of the interest in comparing the 
experimental results with the analytical predictions. In many of the SSME 
measurements, the fLowfield is highly complex due to the sharp turns and 
bends, valves, pump and turbine inlets and discharges, and other 
complicating factors. These factors accentuate the difference between the 
physical quantity at the sensor and the quantity for which the measurement 
is desired, typically an average value at a cross-section. These assessments 
require extensive review of the measurement, the sensor and its installation, 
the thermodynamic and fluid dynamic flowfield, and their interaction. 

The conceptual bias uncertainties can be estimated by reviewing the 
results of other analyses and reviewing the test data or computational model 
results. For example, data from the cold-flow testing of the HPFTP turbine 
shows a temperature profile at the turbine exit and this information can be 
used to estimate the temperature profile which might exist downstream at 
the sensor position. 

For the conceptual bias uncertainty associated with RTD temperature 
measurements a simple 1-D heat transfer analysis can be done to get a rough 
estimate of the temperature profile caused by the heat transfer through the 
duct walls. This analysis provides an order-of-magnitude estimate and helps 
determine if the conceptual bias uncertainty is significant or not. 


2.5 Discussion of Temperature Uncertainties 

The fluid temperature measurements used in this effort were obtained 
using RTD temperature probes. RTD probes typically provide very accurate 
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measurements and the primary systematic uncertainty sources are 
calibration, data acquisition, and conceptual biases. The precision 
uncertainty estimates were obtained as previously discussed and represent 
the within-test and test-to-test precision variation. 

The systematic uncertainty estimates for the calibration and data 
acquisition system were based upon the information in the TTB pre-test Data 
Review Documents (DRD). The temperature calibration and data acquisition 
system systematic uncertainties in the DRD were based upon a Sverdrup 
Technology study by Mr. James Fish 4 . The conceptual bias uncertainty 
added to the DRD listed uncertainty was based upon engineering judgment, 
since little information was available to make a more precise estimate. The 
RTDs typically used in the TTB instrumentation extend into the flowfield, up 
to approximately 2 inches, and are designed so that some of the spatial 
variation of the flow will essentially average out. Thus, the conceptual bias 
represents an estimate of the difference between the measured temperature 
and the one-dimensional average temperature at that location. Thus, a total 
systematic uncertainty of 2% of the measured value is used for the majority 
of the temperature measurements. This value is probably conservative and 
when experimental or computational results are obtained which provide 
additional information about the temperature profiles at the instrumented 
locations, the conceptual bias estimate can be updated. 


2.6 Discussion of Pressure Uncertainties 

The pressure measurements used in this effort were obtained 
using various pressure transducers, however usually strain-gage type 
transducers were used. The primary systematic uncertainty sources Eire 
calibration, data acquisition, and conceptual biases. The precision 
uncertainty estimates were obtained as previously discussed and represent 
the within-test and test-to-test precision variation. 

The systematic uncertainty estimates for the calibration and data 
acquisition system were based upon the information in the TTB pre-test Data 
Review Documents (DRD). The pressure calibration and data acquisition 
system systematic uncertainties in the DRD were based upon the Fish 
study 4 . The conceptual bias uncertainty added to the DRD listed uncertainty 
was based upon engineering judgment, since little information was available 
to make a more precise estimate. The pressure transducers are typically 
mounted along a duct wall. The location of the pressure measurements and 
the geometry of the engine components and ducts provides the potential for 


4 Fish, James, E., “NASA/MSFC Test Area Measurement System Uncertainty Study,” Sverdrup 
Technology, Inc., MSFC Group, Report no. 335-002-92, October 1992. 
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substantial pressure variations in the flowfield. Connected directly to the 
discharge of the High Pressure Oxidizer Pump (HPOP) is a venturi to 
determine the flowrate. The pressure is measured at the entrance of the 
venturi and another pressure transducer is located on the side of the duct in 
the middle of the duct turn. A pressure profile will exist at both of the 
pressure measurement locations. Thus, the conceptual bias represents an 
estimate of the difference between the measured pressure and the one- 
dimensional average pressure at that location. Thus, total systematic 
uncertainties between 1.5% and 2% of the measured value are used for the 
majority of the pressure measurements. These values are probably 
appropriate, however when experimental or computational results are 
obtained which provide additional information about the pressure profiles at 
the instrumented locations is obtained, the conceptual bias estimate can be 
updated. 


2.7 Discussion of Venturi Systematic Uncertainties 

The previous contract report 1 discussed in detail the estimation of the 
venturi flowrate uncertainties. The systematic uncertainties associated with 
the venturi flowrate determinations are between 2% and 3%. The primary 
uncertainty source was found to be the uncertainty associated with the 
venturi discharge coefficients. The venturis were calibrated with water at 
room temperature at maximum Reynolds numbers less than 10% of the 
Reynolds number during engine operation 5 . The difference between using 
water to simulate cryogenic propellants and the Reynolds number 
extrapolation for discharge coefficient are sources of systematic uncertainty. 

It was also observed during the previous effort that the precision 
uncertainty of the flowrate measurements is a function of how the engine 
balances during the specific test. Meaning the precision uncertainty of a 
flowrate measurement represents the variation in the engine operation and 
not a measurement related variation. 


2.8 Discussion of Valve Position Uncertainties 

The systematic uncertainty estimate for the valve position 
measurements was used directly from the pre-test data review documents, 
0.5%. The RVOTs were manufactured by Moxon, Inc., however attempts to 
contact Moxon for uncertainty information were unsuccessful. The precision 
uncertainty estimates based upon the baseline data set and the data 


5 Lepore, Frank A., Rocketdyne Division, Space Shuttle Main Engine No. 3001, Technology Test Bed, 
Differential Flowmeters Calibration Final Report, Contract No. NAS8-27980, March 1980. 
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reduction data sets indicates that the valve position test-to-test variation is 
the dominant uncertainty characteristic. 


2.9 Discussion of Speed Measurement Uncertainties 

Little information was obtained upon which to base turbopump speed 
measurement uncertainty estimates. A nominal 1% systematic uncertainty 
was estimated based upon engineering judgment. This is probably a very 
conservative value and should be updated when more information is obtained 
to support a better estimate. 
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3.0 Modeling Uncertainty Considerations 

The numerical optimization and numerical solution strategy developed 
by Dr. Santi for the New SSME Model requires the uncertainty associated 
with each balance relation. The balance relations are used in the numerical 
solution algorithm to ensure that the conservation of mass and energy are 
within a certain preset tolerance at each node of the system. The numerical 
solution strategy in the new model uses uncertainty estimates to replace the 
numerical tolerances used in the basic ROCETS model. This section 
discusses how uncertainty estimates for the engine subcomponent models 
were obtained. A total of 92 balance relations are used in the current version 
of the model, the balance relations and the uncertainty estimates for each 
balance relation are shown in Appendix 5. 

No accepted methodology currently exists for the assessment of 
uncertainties from analytical models. Within this research effort 
considerable progress has been made towards addressing this issue. 

When comparing output of a model with experimental data, the 
uncertainties that should be associated with the model predictions must be 
considered for proper conclusions to be drawn. In the past, most of the work 
reported in this area has simply considered the sensitivity of the model output to 
uncertainties in the input data. This obviously does not include any 
uncertainties in the model itself and thus is not a satisfactory approach. In this 
research effort, we have divided the sources that cause uncertainty in the model 
output into three categories: (1) uncertainties due to assumptions and 
approximations in the model, (2) uncertainties due to the incorporation of 
previous experimental data into the model, and (3) uncertainties due to the 
numerical solution algorithm. Uncertainties due to the numerical solution 
algorithm are not considered in this effort because the magnitude of those 
uncertainties are much, much less than the uncertainties in the test data and 
the uncertainties in the subcomponent models. 

3.0 1 Modeling Assumptions and Approximations 

When a model of a physical system is developed, assumptions and 
approximations about the system are made to simplify the system to one 
which mathematical expressions can describe. By making these 
simplifications an error is introduced and the model cannot exactly describe 
the physical system. Some of the primary assumptions and approximations 
made within the subcomponent flow models include: 1 -dimensional, fully 
developed, steady-state, adiabatic, ideal gas, invisdd, etc. 

If the uncertainty to associate with a particular assumption or 
approximation can be estimated, then sometimes the model should be 
improved to include this estimate instead of trying to estimate the 


14 




uncertainty. For example, Dr. Santi determined that a turbine exit 
temperature was being predicted using an ideal gas, constant specific beat 
approximation, which for the specific temperature range of interest was a 
poor approximation. Instead of trying to estimate an uncertainty to associate 
with that approximation, the model was altered to include a better 
thermodynamic description of the process. 

The first category, uncertainties due to assumptions and approximations 
in the model, does not include the installation and/or conceptual bias sources 
since those uncertainties are associated with the measured value. Consider the 
temperature at a particular position in the flow. The uncertainty associated 
with the measured value of the temperature includes the effect of making a 
point measurement but desiring a cross-sectional averaged value. The inability 
of the model to calculate a correct average temperature at a particular location 
because the one-dimensional flow approximation has been made results in an 
uncertainty in the predicted temperature. (Stated another way, if the model 
predicts the correct average temperature at a particular location, then the one- 
dimensional flow approximation has caused no uncertainty in the model output.) 

If the individual engine subcomponent models are developed based 
upon the soundest assumptions and approximations available, assessing the 
uncertainty due to these assumptions and approximations would require an 
effort beyond the scope of the current research program. 


3.02 Uncertainties from Previous Information 

Uncertainties due to the incorporation of previous experimental data in 
the model arise when physical property data is used, when valve resistance 
characteristics are used, when turbopump performance maps are used, etc. 
These are all instances in which previous experimental data has been used by 
replacing the data with curvefits. The original data contained uncertainties, but 
the curvefit equations used in the predictive models have been treated as the 
"truth" in most previous considerations of uncertainty in model outputs. Adding 
further complication, there is no accepted way of estimating the influence of 
systematic uncertainties on the uncertainty associated with a regression. The 
methodology developed as a part of this program to assess regression 
uncertainties is discussed later. 

In all of the sub-component modifies, information from previous 
testing is used. For example, the model of the liquid oxygen flow through a 
duct or through a valve is based upon its component testing, which provides 
an equation for the resistance through the duct as function of the flowrate. 
This test information is often reduced to the form of a line or curve. The 
values of the polynomial constants in the thermodynamic property routines 
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are also examples of using previously obtained uncertain test information in 
a model. 

The methodology to assess the uncertainty in the coefficients of a 
linear regression were developed as part of this effort. The uncertainty 
analysis methodology presented in Appendix 1 was applied to the expressions 
for the regression coefficients to develop the technique. The details of this 
methodology were presented at the 1995 AIAA Aerospace Sciences 
conference 6 . The work presented in that paper demonstrated this technique 
provides the uncertainty in linear regression coefficients, and properly 
incorporates the effects of correlated systematic uncertainties. The research 
conducted during the AFOSR research program extended the methodology to 
provide the appropriate uncertainty interval for the predicted value using the 
regression model equation. The application of the regression uncertainty 
methodology to the types of models and previous experimental information 
found in the SSME program is shown in Appendices 3 and 4. 

This chapter will discuss specific aspects of the assessment of the 
balance relation uncertainties. The only uncertainty sources being 
considered in the balance relation uncertainties are those due to the 
uncertain physical property data and the uncertainty due to the hardware 
characteristics, no uncertainties due to modeling assumptions and 
approximations were considered. The only exception to that statement is 
that an overall uncertainty for the heat flux correlation model is being 
estimated instead of propagating uncertainties for the individual constants. 


3.1 Discussion of Physical Property Data Routine 

Uncertainties 

When the model calls the physical property data subroutine to obtain 
the value of a state variable, density, enthalpy, entropy, etc., an uncertainty 
is introduced because the property data model is not the truth. The property 
data subroutines use data generated using curvefit equations to model 
experimental data. The thermophysical property package used in ROCETS 
is a property routine developed by Pratt & Whitney called PROP05. 7 
PROP05 contains thermophysical property subroutines for Para- 
Hydrogen(H2), Oxygen (02), steam (ST), Nitrogen (N2), and Methane (ME). 
The P&W letter introducing the availability of PROP05 states that the 02 


6 Brown, Kendall K., Coleman, Hugh W., and Steele, W. Glenn, “Es timating Uncertainty Intervals for 
Linear Regression,” AIAA Paper 95-0796, 33rd Aerospace Sciences Meeting and Exhibit, Reno, NV, 
January 9-12, 1995. 

7 Long, R.A., and Perry, Michael, J., Internal Correspondence: Introduction of Property Package 

PROP05, August 16, 1988. 
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and ST tables are based on NASA Lewis database, which is provided in the 
programs GASP 8 and WASP 9 , respectively. It also states that the “H2 tables 
are based on NBS database up to about 600 degree R and GASP database 
from 600 degree R to 5400 degree R.” 

The GASP documentation shows the uncertainty in density for Oxygen 
to be within ±0.2% outside the critical region and within approximately 
±2.0% near the critical region. The WASP documentation shows the 
uncertainty in density for Steam to be within ±0.25% for the conditions of 
interest. The GASP documentation shows the uncertainty for Hydrogen 
outside of the critical region to be within ±0.2%. The NBS data 10 for low 
temperature hydrogen shows the uncertainty for density to be within ±0.1% 
outside the critical region and within approximately ±6.0% near the critical 
region. However, it is interesting to note that the GASP documentation 
shows that the uncertainty for the H2 property data near the critical region 
to be within ±2.5%. McCarty indicates the GASP routine uses a modified 
Benedict Webb Rubin (MBWR) equation of state and “the main disadvantage 
of the MBWR is that it is functionally incorrect in the critical region, i.e. 

p = p c ± 20% p c . 


3.2 Discussion of Combustion Routines 

Related to the property routine uncertainty is the uncertainty in the 
combustion routines. Combustion occurs in the fuel and oxidizer prebumers 
and the main combustion chamber. The combustion occurs at very high 
pressures, over 5200 psia in the prebumers and 3126 psia in the main 
combustion chamber. Much of the combustion process is fuel-rich. The 
documentation available on the combustion routines in ROCETS is limited, 
the following co mm ents are in the computer code: 

1. THESE COMBUSTION CURVES WERE GENERATED FROM BRINKLY PERFECT 
GAS COMBUSTION DECK. 

2. AN UNCOMBUSTED H2/02 MIXED TEMPERATURE IS CALCULATED USING 
CONSTANT VALUES OF CP’S. THESE CP'S ARE CONSISTENT WITH THE VALUES OF 
CP’S USED TO GENERATE THE COMBUSTION PROPERTIES. 


8 Hendricks, Robert C., Baron, Anne K., and Peller, Ddiko C., GASP - A Computer Code for Calculating 
the Thermodynamic and Transport Properties for Ten Fluids: Parahydrogen Helium, Neon, Methane, 
Nitrogen, Carbon Monoxide, Oxygen, Fluorine, Argon, and Carbon Dioxide, NASA TN D-7808, 
February 1975. 

9 Hendricks, Robert C., Peller, Ildiko C, and Baron, Anne K., .WASP -A Flexible Fortran IV Computer 
Code for Calculating Water and Steam Properties, NASA TND-7391, November 1973. 

10 McCarty, Robert D., Hydrogen Technological Survey - Thermophysical Properties, NASA SP-3089, 
Washington, DC, 1975. 
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3. NEXT THE H2/02 COMBUSTION PROPERTIES ARE LOOKED UP AS A FUNCTION 
OF OFR AND MIXED TEMPERATURE. THESE COMBUSTION PROPERTIES WERE 
GENERATED AT 2500 PSI, AND WERE GENERATED FOR 5 MIXED H2/02 
TEMPERATURES OVER AN O/F RANGE OF .25 TO 100. A PRESSURE CORRECTION 
IS THEN APPLIED TO THE COMBUSTION TEMP. THE MAP INPUTS AND RESPECTIVE 
RANGES ARE: 

P - 100. TO 5000. PSI 

TMIX - 200 TO 1000. DEG R (FOR BURN OPTION) 

TC - 200 TO 6998. DEG R (FOR UNBURN OPTION) 

OFR - 0. TO 1.0 

4. THE COMBUSTION PROPERTIES ARE THEN DILUTED FOR HELIUM. 

An uncertainty for the combustion properties will exist based upon the 
assumptions and approximations and the property data used to build the 
combustion routine. One of the primary assumptions upon which the 
combustion routines are based is Dalton’s partial pressure model, at high 
pressures and temperatures the Dalton model is not very accurate because of 
real gas effects. The constant specific heat assumption is also questionable 
at high temperatures and pressures. The injector characteristics, such as 
atomization, and heat transfer effects also introduce diffi culties in the 
combustion modeling. From the information available, it could not be 
determined how well the other aspects of the combustion routine account for 
the real-gas and other real effects in high pressure, high temperature, fuel 
rich combustion processes. Thus, it is difficult to provide uncertainty 
estimates for the combustion properties, with the likely uncertainty in the 
10% to 20% range. 


3.3 Discussion of Hardware Characteristics Uncertainties 

In general, the hardware characteristics used in the New SSME Model 
were taken from the Rocketdyne SSME Power Balance Model, version PBM 
91a. The hardware characteristics are specific to each component and in 
most cases are based upon sub-component level testing. For example, the 
turbines were cold-flow tested to determine their non-dimensional 
performance and this information, after adjusting to engine operating 
conditions, is used in the model. Since this information was experimentally 
obtained it contains experimental uncertainties and obtaining or 
independently assessing the uncertainties was necessary for this effort. The 
following sections discuss the information obtained and the assessment of 
hardware characteristic uncertainties. 
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3.3.1 Discussion of Turbine Map Uncertainties 

The performance of the turbines are represented in the models by 
using performance maps. These maps are based upon the available 
subcomponent test data or other analyses, and as such they are not perfect 
representations of the turbine performance. The uncertainty associated with 
each map must be estimated based upon how that particular map was 
generated. 

3.3.1. 1 HPFT Map Uncertainty 

The HP FTP turbine efficiency map used in the ROCETS SSME model 
was obtained from the Digital Transient Model (DTM) and the SSME Power 
Balance Model, PBM 91a. Where the portion of the turbine efficiency 
parameter versus speed parameter map for the nominal operating points of 
the SSME (i.e. 100%, 104%, 109%) was taken from the Power Balance Model, 
version 91a The remaining portion of the map was taken from the DTM and 
is only used during transient simulations, engine starts and shutdowns. As 
discussed more thoroughly in Appendix 4, the high pressure fuel turbine was 
tested in the Turbine Test Equipment facility at MSFC in 1991. 11 The 
results of the cold-flow testing were analyzed by Rocketdyne engineers and 
the cold-flow turbine efficiency map was decreased by 3.2% to adjust to 
engine operating conditions 12 . The conversion of the air test results for 
SSME engine operation were made to account for thermal effects, Reynolds 
number effects, disk cooling, and platform seal leakage effects. Since each of 
these corrections is based upon a number of assumptions and approximation 
they introduce additional uncertainties, so an uncertainty source to account 
for the uncertainty in the correction terms must be included. 

When the results from this testing, after applying the corrections, are 
plotted in the same form as the HPFT map in ROCETS, turbine efficiency 
divided by speed parameter versus speed parameter, Figure 1 is obtained. 
Since the map obtained from the cold-flow testing of the HPFT model 
provides essentially the same results at the mainstay operating conditions, it 
is proposed that the uncertainty associated with this map be used as the 
uncertainty in the HPFT balance relation. Details of how the uncertainty 
associated with a value from the air-test HPFT map are given in Appendix 3. 


11 Hudson, Susan T., Gaddis, Stephen W., Johnson, P., Dean, and Boynton, Janies, L., “Cold Flow 
Testing of the Space Shuttle Main Engine High Pressure Fuel Turbine Model,” AIAA Technical Paper 
91-2503, AJAA/SAE/ASME/ASEE 27th Joint Propulsion Conference, June 24-26, 1991, Sacramento, 
CA. 

12 Boynton, J., and Daumann, A., “Revision of SSME HPFTP Turbine Performance Based on 
NASA7MSFC Baseline Air Test Results,” July 23, 1991. 
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The total uncertainty associated with the HPFT map is the 
combination of the uncertainty from the cold-flow testing and the uncertainty 
in the corrections to engine conditions. This data reduction equation is 

ETA map = ETAjjgt +X cor x ETA test 

= ETA test x(l+X cor ) () 


The uncertainty expression for this equation is 
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seta^ 


\2 


seta^; 
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( 2 ) 


and after determining the partial derivatives and algebraically reducing the 
above equation the uncertainty expression can be written as 

U£TAmp = J Uet ^ + - U ^ / (3) 

ETAnup [ETAj^ (l+X cor ) J 

Using the uncertainty in the cold-flow efficiency map determined in 
Appendix 3 of 0.75% and making a conservative estimation that the 
corrections applied to adjust to engine conditions is within 25%, the 
uncertainty in the map efficiency is 1.1%. 



Figure 1 Comparison of HPFT Maps, map existing in ROCETS and map based upon cold-flow 
testing of the HPFT. 


This uncertainty estimate represents the uncertainty in the 
performance on a high pressure turbine with the same characteristics as the 
model used in the cold-flow tests. This uncertainty does not include any 
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uncertainty due to manufacturing variations or any uncertainty due to wear- 
related performance degradations. By comparing this uncertainty estimate 
with the uncertainty estimates from the turbopump pump maps, which 
include these additional uncertainty sources, the uncertainty for the HPFT 
performance maps is increased to 4.0%. 

3.3.1.2 Other Turbine Performance Map Uncertainties 

Similar information was obtained for the high pressure oxidizer 
turbine; however, the HPOT test program was conducted in 1974 and while 
the cold-flow tests were conducted similarly to the HPFT tests, the 
experimental uncertainty would be expected to be greater. 13 Uncertainty 
estimates for the HPOT performance maps using this information and 
including additional uncertainty to account for manufacturing variability 
provides performance map uncertainties of 5%. 

No information was obtained upon which to base uncertainty 
estimates for the low pressure fuel and oxidizer turbines., LPOT, and LPFT. 
Thus, uncertainty estimates of 5% were used for the turbine performance 
maps. 


3.3.2 Discussion of Turbopump Pump Map Uncertainties 

Information concerning the SSME turbopump pump maps was 
requested from the Turbomachinery Branch, Rocketdyne Division. Pump 
performance maps with ±2o bands representing the performance variation 
caused by manufacturing variability were provided, however the 
nomenclature and scale of the maps were not the same as the performance 
maps in the ROCETS model. The maps for the HPFP and the LPFP provided 
by Rocketdyne show “Normalized Head” versus “Normalized Flow,” and the 
maps used in ROCETS show “Head Coefficient” versus “Flow Coefficient.” 
The Rocketdyne provided maps for the LPOP and the HPOP use the same 
nomenclature of “Head Coefficient” versus “Flow Coefficient,” but the scales 
and units are different. The relationship between the Rocketdyne provided 
performance maps and the ROCETS maps could not be determined prior to 
the end of this contract. 

In order to make uncertainty estimates for the pump hardware 
characteristics, it is being assumed that the relationships between the two 
sets of maps are linear transformations that account for unit conversions and 
design point normalizations. Making this assumption allows the percent 
uncertainty from the Rocketdyne maps to be directly transferred to the 
ROCETS maps, and that uncertainty estimate will be used to determine the 


13 Boynton, J., Rocketdyne Division, Personal Communication, 27 March 1996. 
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balance relation uncertainty. It is felt that this is a good assumption, with 
the information available If future information shows the relationship is 
non-linear then the hardware characteristic uncertainty for the pump 
performance maps will have to be estimated using the uncertainty 
propagation methodology provided in Appendix 1. 


3.3.3 Discussion of Duct, and Valve Resistances 

The duct and valve resistance values used in the ROCETS model 
originated from a variety of sources, and for most cases the exact origin of 
the resistance or flow coefficient values are not known. The uncertainties to 
associate with the values of the flow resistance hardware characteristics in 
the model must be estimated. Very little information is currently available to 
base these estimates. The ideal way to assess these uncertainties is to 
independently test the particular duct section at engine operating conditions 
and determine the flow resistance. This however is impractical, so the 
resistance values and their uncertainties must be otherwise estimated. 

The resistance values for the ducts which contain the venturi 
flowmeters were taken from the venturi calibration report. 5 The venturi 
flowmeters were calibrated within the duct as installed in the TTB engine as 
an assembly. The resistance for the duct assemblies, with the venturi 
flowmeter, are provided in the venturi calibration report, however the 
uncertainties associated with the experimental determination of the 
resistance values are not provided. The origin of the resistance values of the 
other ducts and for the valves is unclear. The resistance values in the 
Rocketdyne Power Balance Model were obtained by analyzing engine hot-fire 
data and from separate analytical estimates While some of the ROCETS 
resistance values could be obtained from the PBM, the PBM resistances 
cannot be used directly in the ROCETS model because the model flow circuits 
are significantly different 

A review of how uncertainty estimates for the resistance would be 
estimated would be helpful. The flow resistance is defined by the expression 

R= APxRHO 

w 2 

where W is the total mass flowrate, AP is the pressure difference across the 
duct section, and RHO is the density of the fluid. The uncertainty expression 
is then 




and the expression can be written in terms of percentage uncertainty as 



Now the resistance uncertainty can be examined by considering the 
uncertainties in the possible test rigs or test data. As Eq. (6) shows, the mass 
flowrate uncertainty has twice as much influence as the pressure difference 
or the density. As with the venturi calibration testing, the tests would be 
conducted with water at ambient conditions and the data would be 
extrapolated to engine operating Reynolds numbers. As was extensively 
discussed in the previous contract report, additional uncertainties may be 
introduced by not performing the tests with the actual cryogenic fluid. 

3.3.4 Heat Transfer Calculations 


The subroutines QH2 and Q02 calculate the convective heat flux to, or 
from, the fluid in a duct. This flux is multiplied by an effective area to obtain 
the total energy added to or removed from that volume. For example, the 
nozzle coolant heat transfer is calculated by using QH2 to determine the heat 
flux between the hot wall of the nozzle and the liquid Hydrogen coolant. This 
value is then multiplied by an effective heat transfer area to obtain the total 
energy added to the nozzle coolant volume. The heat flux calculation in QH2 
is 


Q' = 
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where W is the weight flowrate, D is the hydraulic diameter, TF is the fluid 
temperature, TM is the metal temperature, and Cp , p, and k are fluid 
specific heat, viscosity, and thermal resistance, respectively. A different heat 
flux equation is in QH2 if coolant boiling occurs. The total heat transfer Q is 
then found with 


Q = Q'*A eff (8) 

The origin of the heat flux correlation is not known, there is no 
documentation in the QH2 subroutine. The uncertainty for the heat flux 
equation can be estimated in two ways. First, uncertainties could be 
estimated for each of the parameters and constants in Eq. (6) and then 
propagated to obtain an uncertainty in Q\ But, probably the most 
appropriate way is to estimate the uncertainty for the overall heat transfer 
correlation by examining how the heat flux equation was developed. Since 
no information is available for these heat flux correlations, uncertainty 
estimates must be estimated based upon engineering judgment and 
knowledge of the accuracy of other heat transfer correlations. It is obvious 
that the uncertainty will be rather large by noting that the uncertainty in C P , 
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|i, and k would be large for super-critical Hydrogen. Thus, heat flux 
uncertainties of 20% will be used. 

The value of the effective area is also subject to uncertainty. Again 
looking at the nozzle coolant heat transfer, since the nozzle is made up of 
over a thousand individual tubes an accurate estimation of the actual heat 
transfer area would be difficult. The uncertainty to associate with the 
effective area must be estimated based upon how detailed the area 
calculation was. 


3.4 Module Calculations and Uncertainty Expressions 

The uncertainties from the physical property data and the hardware 
characteristics are propagated to obtain the uncertainties for the result of 
each module using the methodology presented in Appendix 1. This section 
will briefly present the equations used in each of the major subcomponent 
modules and how the hardware characteristics are propagated to obtain the 
uncertainty for the module results. The uncertainty expressions and 
calculations for each subcomponent are shown in Appendix 4. 


3.4.1 Duct and Valve Flow Modules 

The equation used in the FLOWOO module for the flowrate, W, through 
a duct is 


w = 


I 


AP xRHO 
R 


(9) 


where AP is the pressure difference across that section of duct, RHO is the 
density of the fluid in the duct, and R is the hardware characteristic for the 
flow resistance of the duct. 

The expression for the uncertainty in the flowrate due to uncertainties 
in density from the property routine and uncertainty in the flow resistance 
value. 


u 


w 



( 10 ) 


and the expression can be written in terms of percentage uncertainty as 



The calculations and the uncertainty calculations for the fuel and 
oxidizer valves are performed in the same manner as for the ducts and lines. 
The primary difference is that the resistance for the valve is a function of the 
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valve position and the resistance of the valve is added to the resistance value 
for the duct in which the valve is located. 

Some of the ducts and valves use a slightly different model, they use a 
flow coefficient, CF, instead of a resistance. The equations for these ducts 
and valves have the form 

W = ^AP x RHO x CF 2 x FLOCON (12) 

where FLOCON is defined as the flow constant, and considering 
uncertainties in density and the flow coefficient the uncertainty propagation 
for this is 



Comparing Eq. (13) to Eq. (11) shows that the uncertainty in flow coefficient 
has a direct relationship on the flow uncertainty, while only half of the 
uncertainty in resistance is propagated into the flow uncertainty. 


3.4.2 Turbine Modules 


The mass flowrate from the turbine is obtained directly from a map 
and thus the uncertainty is that of the map. 

The energy from the exit of the turbine is found from the expression 

HTOUT = HTIN - DH (14) 


where DH is the predicted enthalpy drop across the turbine and obtained 
from the turbine maps for ETAmap and DHIDLmap and the expression 

DH = (ETA map )(DHIDL map ) (15) 


where 


DHIDL map = HTEN-HTOUTI 

= HTIN - f(SIN,PTOUT) 


(16) 


and HTOUTI is the ideal exit enthalpy and is determined as an isentropic 
turbine process as a function of the exit pressure, inlet entropy, and the 
thermophysical property routine. Thus the energy balance relation 
uncertainty is a function of the uncertainty in the efficiency map and the 
physical property routine uncertainties. The uncertainty expression in 
percentage form is 
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(17) 
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Using the efficiency map uncertainty of 1.1% and an uncertainty of 0.5% 
from the thermophysical property package, the total uncertainty in DH is 


1 . 2 %. 


3.4.3 Pump Modules 


The equation for the exit pressure from a pump is 

PTout=PTin+PSI x RHO x SNRAD 2 (18) 


where PSI is the value from the head coefficient map at the flow coefficient, 
PHI, determined from the expression 


PHI = 


Win 

RHOx SNRAD 


(19) 


and where Win, is the flowrate into the pump, and SNRAD is the pump speed 
in radians per second. The power required by the pump, QDOT, is calculated 
with 

QDOT = -TORQ x SNRA %j (20) 

where TORQ is the pump torque. The pump torque is determined in the 
model with the equation 

TORQ = -TAUxRHOxSNRAD 2 (2 1) 

where TAU is from the torque map parameter versus flow coefficient 
performance map. 


The flows through the pumps are iteration variables, so the pump 
discharge pressure and the power required by the pump are the results used 
in the balance relations. 


3.5 Balance Relation Uncertainties 

The ROCETS model uses conservation of mass and conservation of 
energy at each flow circuit node to balance the system in the numerical 
solution algorithm. The 92 balance relations are shown in Appendix 4, as 
provided by Dr. Santi. In the data reconciliation optimization strategy 
developed by Dr. Santi the balance relation numerical tolerances are 
replaced by balance relation uncertainties based upon the hardware 
characteristic uncertainties and the physical property data uncertainties. 
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The balance relation can be rewritten by subtracting the left-hand side 
from the right-hand side, for the continuity and energy balance relations can 
be rewritten as 


BAL# = 2]Wo l rr-EW IN 


( 22 ) 


and 


B AL# = 2 Energy our - £ Energy^ 


(23) 


The balance relation uncertainty is determined by using the balance 
relations in this form and applying the uncertainty propagation methodology 
as presented in Appendix 1. In the solution algorithm certain parameters, 
pump flowrates, pressures, and enthalpies are iteration variables. No 
uncertainties are being considered for iteration variables. 

Appendix 4 is the output from a MathCAD V.6.0+ file which contains 
the equations and calculations for the module parameters and uncertainties 
and the balance relation uncertainties. The balance relation equations and 
the balance relation uncertainties are summarized in Table 1 in Appendix 4. 
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4.0 


Summary 


This report has discussed the primary areas of uncertainty in the 
testing and modeling of a new Space Shuttle Main Engine model and 
performance analysis tool. The uncertainty estimates in this report will 
support the development of the new data integration and reconciliation 
model. The overall uncertainty estimates for the test data indicate 
systematic uncertainties are very important, and are often the dominate 
uncertainty characteristic. The development of the data reduction data sets 
indicated precision uncertainty estimates must be performed carefully to 
ensure proper comparisons. The work conducted to assess the uncertainties 
in the ROCETS subcomponent models was a unique extension of the 
uncertainty analysis methodology. It is now apparent that significant 
uncertainties exist in the hardware characteristics due to the uncertainty in 
the information used in the development of the characteristics. A 
methodology for assessing the uncertainty in linear regressions was 
developed and used to assess the uncertainty in the HPOT efficiency map. 


4.1 Uncertainty Estimate Improvements 

The uncertainty estimates provided in this report represent the best 
information available at the time of this report. They should not be 
considered final uncertainty estimates, they should be updated as more and 
better information becomes available. The uncertainty analysis 
methodologies presented in the appendices should be used to improve these 
estimates. 

The following items are areas that need further investigation, or are 
areas where information was obtained too late in the performance period to 
allow proper analysis and inclusion in this report. 

1. The balance relation uncertainty estimates for balance relations 50- 
69 and 79-93 need to be completed. These balance relations involve the 
combustion routines and involve the hot gas ducts and engine 
components. The complexity of many of these relations and time did not 
allow detailed investigation of these balance relations. 

2. The pump maps provided hy Rocketdyne do not correlate with the 
maps in the ROCETS model. The exact correlation of these maps and the 
map uncertainties must be obtained and the performance map 
uncertainty estimates updated, if necessary. The information provided by 
Rocketdyne is included in Appendix 7. 
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3. No uncertainty information was obtained for the LPOTP turbine 
and the LPFTP turbine. This information should be obtained and the 
balance relation uncertainty calculations updated. 

4. The only uncertainty information obtained for HPOTP and HPFTp 
turbine performance maps was from cold flow testing of a single turbine 
model. Information describing the turbine performance variation based 
upon manufacturing variability is needed to update the estimates used in 
the balance relation uncertainty calculations. 

5. The conceptual bias estimates for the static pressure measurements 
a the inlet and exits of the high pressure fuel and high pressure oxidizer 
pumps need to be re-examined. Information was provided by the MSFC 
Aerophysics Branch, Experimental Division, Structures and Dynamics 
Laboratory upon which better estimates of the conceptual bias 
uncertainties can be made. However, this information was obtained too 
late for inclusion in this report. 
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Appendix 1 
Uncertainty Analysis 


The use and application of uncertainty analysis in engineering has 
evolved considerably since Kline and McClintock's classic paper 1 in 1953. 
Developments in the field have been especially rapid and significant over the 
past decade, with the methods formulated by Abemethy and co-workers 2 that 
were incorporated into ANSI/ASME Standards in 1984 3 and 1986 4 being 
superseded by a more rigorous approach 5 . Publication in late 1993 by the 
International Organization for Standardization (ISO) of the Guide to the 
Expression of Uncertainty in Measurement S 6 in the name of ISO and six other 
international organizations has, in everything but name only, established a new 
international experimental uncertainty standard. 

The approach in the ISO Guide deals with "Type A” and "Type B" 
categories of uncertainties, not the more traditional engineering categories of 
bias and precision uncertainties, and is of sufficient complexity that its 
application in normal engineering practice is unlikely. This issue has been 
addressed by AGARD Working Group 15 on Quality Assessment for Wind 
Tunnel Testing and by the Standards Subcommittee of the AIAA Ground Test 
Technical Committee. The documents 6 ” 7 produced by these groups present and 
discuss the additional assumptions necessary to achieve a less complex "large 
sample" methodology that is consistent with the ISO Guide, that is applicable to 
the vast majority of engineering testing (including most single-sample tests), and 
that retains the use of the traditional engineering concepts of bias and precision 
uncertainties. (The chapters on uncertainty methodology in the AGARD 6 and 
AIAA 7 documents were authored by the Principal Investigator of this research 
program.) 


1 Kline, S. I, and McClintock, F. A., "Describing Uncertainties in Single-Sample Experiments," Mechanical 
Engineering , Vol. 75, 1953. 

2 Abemethy, R. B., Benedict, R. P., and Dowdell, R. B., "ASME Measurement Uncertainty," J. Fluids 
Engineering, VoL 107, 1985. 

3 American National Standards Institute/American Society of Mechanical Engineers, Measurement Uncertainty 
for Fluid Flow in Closed Conduits, MFC-2M-1983, ASME, 1984 

4 American National Standards Institute/American Society of Mechanical Engineers, Measurement Uncertainty, 
PTC 19.1-1985 Part 1, ASME, 1986. 

5 International Organization for Standardization, Guide to the Expression of Uncertainty in Measurement, ISO, 
ISBN 92-67-10188-9, 1993 

6 Quality Assessment for Wind Tunnel Testing, AGARD- AR-304, 1994. 

7 American Institute of Aeronautics and Astronautics, Assessment of Wind Tunnel Data Uncertainty, AIAA 
Standard S-071, 1995. 
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2.1 Overview 

The word accuracy is generally used to indicate the relative closeness of 
agreement between an experimentally-determined value of a quantity and its 
true value. Error (d) is the difference between the experimentally-determined 
value and the truth, thus as error decreases accuracy is said to increase. Only in 
rare instances is the true value of a quantity known. Thus, one is forced to 
estimate error, and that estimate is called an uncertainty, U. Uncertainty 
estimates are made at some confidence level ~ a 95% confidence estimate, for 
example, means that the true value of the quantity is expected to be within the 
±U interval about the experiment^y-determined value 95 times out of 100. 

As shown in Figure 1(a), total error d can be considered to be composed of 
two components: a precision (random) component e and a bias (systematic) 
component b. An error is classified as precision if it contributes to the scatter of 
the data; otherwise, it is a bias error. Tt. is assumed that corrections have been 
made for all .wsfo anatic errors whose values are known. The remaining bias 
errors are thus equally as lilcelv to be positive as negative. 





Magnitude of X 

(b) 

Figure 2.1 Errors in the Measurement of a Variable X:(a) two readings; (b) infinite 
number of readings. 


Suppose that we are making a number of measurements of the value of a 
variable X that is absolutely steady. The k and k+1 measurements are shown in 
Figure 1(a). Since the bias is a fixed error, it is the same for each measurement. 
However, the precision error will have a different value for each measurement. 
It then follows that the total error in each measurement will be different, since 
the total error is the sum of the bias error and precision error in a measurement. 

If we continued to take measurements as previously described until we 
had a sample of N readings, more than likely as N approached infinity the data 
would behave as shown in Figure 1(b). The bias error would be given by the 
difference between the mean (average) value m of the N readings and the true 


Al-3 




value of X, whereas the precision errors would cause the frequency of occurrence 
of the readings to be distributed about the mean value. 

As an estimator of (3, a bias limit B is defined 8 . A 95% confidence estimate 
is interpreted as the experimenter being 95% confident that the true value of the 
bias error, if known, would fall within ±B. A useful approach to estimating the 
magnitude of a bias error is to assume that the bias error for a given case is a 
single realization drawn from some statistical parent distribution of possible bias 
errors. For example, suppose a ther mis tor manufacturer specifies that 95% of 
samples of a given model are within ±1.0 C of a reference resistance-temperature 
(R-T) calibration curve supplied with the thermistors. One might assume that 
the bias errors (the differences between the actual, but unknown, R-T curves of 
the various thermistors and the reference curve) belong to a Gaussian parent 
distribution with a standard deviation b=0.5 C. Then the interval defined by ±B 
= ± 2b = ±1.0 C would include about 95% of the possible bias errors that could be 
realized from the parent distribution. (The bias limit is sometimes called the 
"systematic uncertainty".) 

As an estimator of the magnitude of the precision errors (the width of the 
distribution of readings in Figure 1(b)), a precision limit P is defined 8 . A 95% 
confidence estimate of P is interpreted to mean that the ±P interval about a 
single reading of X should cover p 95 times out of 100. (The precision limit is 
sometimes called the "precision uncertainty".) 

In nearly all experiments, the measured values of different variables are 
combined using a data reduction equation (DRE) to form some desired result. A 
good example is the experimental determination of mass flow rate using a 
venturi meter as discussed in the previous contract report 9 . Functionally, the 
mass flow rate is given as 

We = W e (P, T, AP,d, D, a, C D ) (1) 

One can envision that errors in the values of the variables on the right hand side 
of Eq. (1) will cause errors in the experimental result W e . 

A more general representation of a data reduction equation is 

r = r(Xi, X 2 , Xj) (2) 

where r is the experimental result determined from J measured variables X. 
Each of the measured variables cont ains bias errors and precision errors. These 
errors in the measured values then propagate through the data reduction 


8 Coleman, H. W., and Steele, W. G., Experimentation and Uncertainty Analysis for Engineers, Wiley, New 
York, 1989. 

9 Coleman, Hugh W., and Brown, Kendall K., "Impact of Uncertainty on Modeling and Testing,” 
Propulsion Research Center Report # 95-001, Final Report on NASA-Marshall Space Fli ght Center 
Contract NAS8-38609 D O. #106, 10 January 1994 to 30 April 1995. 
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equation, thereby generating the bias and precision errors in the experimental 
result, r. 


If the "large sample assumption" is made 67 then the 95% confidence 
expression for U r becomes 


Ur = 2>?5f + 2 2 Z 

+ ■? Z Z OiOkP* 

i=; i~l k=i+l 


(3) 


where 


Oi = 


dr 

dXi 


(4) 


and where the 95% confidence precision limit for a variable X is estimated as 


Pi = 2 St N>10 


(5) 


and the sample standard deviation is calculated using 


Si 


1 


N-l 


t t(X‘h - x,f 


( 6 ) 


where the mean value is defined as 

Z (Xi)k 

k=l 



(7) 


and Pik is the 95% confidence estimator of the covariance of the precision errors 
in X and Xk, and Bik is the 95% confidence estimator of the covariance of the bias 
errors in X and Xk. 


If we define the bias limit (systematic uncertainty) of the result as 

j j-i j 


Bi = Z +2 Z z 


j=i i=i *=*+; 

and the precision limit (precision uncertainty) of the result as 


j-i j 


P^r = Z + 2 Z Z 


i=l 


i—l k=i+l 


then Eq. (3) can be written as 


Ur = Bi+ Pi 


( 8 ) 

(9) 


(10) 
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and Eqs. (8) and (9) can be viewed as propagation equations for the bias limits 
and precision limits, respectively. 

2.2 Determining Precision limits 

Single Test. When the result is determined from a single test -- that is, 
at a given test condition the result is determined once using Eq. (2) 

r = r(X t , X 2 , Xj) (2) 

and when the Xi's are considered single measurements, then Eq. (9) is used to 
find the precision limit of the result. This situation is often encountered in large 
scale engineering tests in which measurements of the variables are made at a 
given set point over a period that is small compared to the periods of the factors 
causing variability in the experiment. A proper precision limit (one indicative of 
the dispersion of the variable over several cycles of the factors causing its 
variation) cannot be calculated from readings taken over such a small time 
interval. For such data, the measurements) of a variable Xi should be 
considered a single reading -- whether the value of X is the average of 10, 10 3 or 
10 6 readings taken during the short measurement time. In such a test, the value 
for the precision limit to be associated with a single reading would have to be 
based on previous information about that measurement obtained over the 
appropriate timp interval 10 . If previous readings of a variable over an 
appropriate interval are not available, then the experimenter must estimate a 
value for Pi using the best information available at that time 6 ' 7 . 

For single tests in which some of the variables (X 2 and X, for instance) 
can be determined as averages from multiple readings over an appropriate time 
period but the other variables cannot be, then 

r = r(X,, X 2 , X 3 Xj) (11) 

and Eq. (9) is used to find the precision limit of the result as follows. For the 
variables that are single readings, the Pi's are the precision limits determined 
from previous information or estimated from the best available information. For 
the averaged variables when N2 and N3 are equal to or greater than 10, P2 and 
P3 should be taken as precision limit s of means, (2S2)/(N2) V2 and (2Ss)/(N3) 1/2 , 
with the S's calculated using Eq. (6). When N2 and N3 are less than 10, it is the 
authors' recommendation that the precision limits used in Eq. (9) for the 
averaged variables be taken as (P^Na) 1 * and (Pa)/^) 172 , where P2 and P3 are 
determined from previous information, as is done for the single reading 
variables. 


10 Steele, W. G., Taylor, R.P., Burrell, R. E., and Coleman, H. W., “The Use of Data from Previous 
Experience to Estimate the Precision Uncertainty of Small Sample Experiments,” A1AA Journal , Vol. 31, 
No. 10, 1993. 
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For tests in which multiple readings of all of the variables can be obtained 
over an appropriate period, the following method is recommended. 


Multiple Tests. If a test is performed so that M multiple sets of 
measurements (Xi, X 2 , ... Xj)k at the same test condition are obtained, then M 

results can be determined using Eq. (2) and an average result r can be 
determined using 

1 M 

r = X7 Z (12) 

If the M sets of measurements were obtained over an appropriate time period, 
the precision limit that should be associated with a single result would be 

Pr = tSr (13) 

where t is determined with M-l degrees of freedom and is taken as 2 for M>10 
and Sr is the standard deviation of the sample of M results 


Sr = 


1 


M-l 


X (*"* 


k=l 


-, 1/2 


rf 


(14) 


The precision limit that should be associated with the average result is given by 


Pr 


Pr 

4m 


(15) 


with P r given by Eq. (13). Using the large sample assumption, the uncertainty 
that should be associated with a single result would be 


U] =B?+(2S r ) 2 

(16) 

and with an average result r 


U\ = B 2 r+(2Sr />[M) 2 

(17) 

with Br given by Eq. (8). 



Correlated Precision Uncertainties. The Pa terms in Eq. (3) take into 
account the possibility of precision errors in different variables being correlated. 
These terms have traditionally been neglected 13 - 4 ' 5 - 7 , although precision errors in 
different variables caused by the same uncontrolled factor(s) are certainly 
possible and can have a substantial impact on the value of the precision limit 11 . 
In such cases, one would need to acquire sufficient data to allow a valid 
statistical estimate of the precision covariance terms to be made if using Eq. (3). 


11 Hudson, S. T., Bordelon, W., and Coleman, H. W., “Effect of Correlated Precision Errors on the 
Uncertainty of a Subsonic Venturi Calibration,” AIAA-95-0797, 1995. 
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Note, however, that the multiple tests approach using Eq. (14) implicitly 
includes the correlated error effect - a definite advantage when multiple sets of 
measurements over an appropriate time period are available. 


2.3 Estimating Bias Limi ts 

Bias Limits of Individual Variables. When attempting to estimate 
the bias limits Bi of the individual variables in Eq. (8), one might, separate the 
bias errors which influence the measurement of a variable into different 
categories: calibration errors, data acquisition errors, data reduction errors, test 
technique errors, etc. Within each category, there may be several elemental 
sources of bias. For instance, if for the Jth variable, Xj, there are M elemental 
bias errors identified as significant and whose bias limit s are estimated as (Bj)i, 
(Bj) 2 , ..., (Bj)m, then the bias limit for the measurement of Xj is calculated as the 
root-sum-square (ESS) combination of the elemental limits 


Bj 


t(Bj k 


(18) 


The elemental bias limits, (Bi)k, must be estimated for each variable Xi 
using the best information one has available at the time. In the design phase of 
an experimental program, manufacturer's specifications, analytical estimates 
and previous experience will typically provide the basis for most of the estimates. 
As the experimental program progresses, equipment is assembled, and 
calibrations are conducted, these estimates can be updated using the additional 
information gained about the accuracy of the calibration standards, errors 
associated with the calibration process and curvefit procedures, and perhaps 
analytical estimates of installation errors. 


As Moffat 12 suggests, there can be additional conceptual bias errors 
resulting from not measuring the variable whose symbol appears in the data 
reduction equation. An example would be a point temperature measurement 
interpreted to be indicative of a cross-section averaged temperature, but there 
may be a cross-sectional variation of temperature, which may or may not have a 
predictable profile, causing the “average” value to be different than the point 
value. Hence, the inclusion of an elemental bias term for the conceptual error 
would be appropriate. 


Correlated Bias Limits. Correlated bias limits are those that are not 
independent of each other, typically a result of different measured variables 
sharing some identical elemental error sources. It is not unusual for the 
uncertainties in the results of experimental programs to be influenced by the 
effects of correlated bias errors in the measurements of several of the variables. 


12 Moffat, R. J., “Describing the Uncertainties in Experimental Results,” Experimental Thermal and Fluid 
Science, Vol. 1, 1988. 
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A typical example occurs when different variables are measured using the same 
transducer, such as multiple pressures sequentially ported to and measured 
with the same transducer or temperatures at different positions in a flow 
measured with a single probe that is traversed across the flow field. Obviously, 
the bias errors in the variables measured with the same transducer are not 
independent of one another. Another co mm on example occurs when different 
variables are measured using different transducers all of which have been 
calibrated against the same standard, a situation typical of the electronically 
scanned pressure (ESP) measurement systems in wide use in aerospace test 
facilities. In such a case, at least a part of the bias error arising from the 
calibration procedure will be the same for each transducer, and thus some of the 
elemental bias error contributions in the measurements of the variables will be 
correlated. 

The Bik terms in Eq. (8) must be approximated -- there is in general no 
way to obtain the data with which to make a statistical estimate of the 
covariance of the bias errors in X and the bias errors in X- The approximation 
of such terms was considered in detail in Ref. 12, where it was shown that the 
approach that consistently gives the most satisfactory approximation for the 
correlated bias limits was 


Bm - £ (BJJBO, (19) 

a~l 

where L is the number of elemental systematic error sources that are common 
for measurements of variables X and X. 

If, for example, 


r = r(Xi, X 2 ) (20) 

and it is possible for portions of the bias limits Bi and B 2 to arise from the same 
source(s), then Eq. (8) gives 

Bi = e 2 ,Bi+0 2 2B2+20 1 0 2 Ba (21) 

For a case in which the measurements of Xi and X are each influenced by 4 
elemental error sources and sources 2 and 3 are the same for both Xi and X, Eq. 
(18) gives 

Bi = (Bi? 1 + (Bif 2 +(Bif 3 + (Bi? 4 (22) 


and 


13 Brown, K. K., Coleman, H W., Steele, W. G., and Taylor, R. P., "Evaluation of Correlated Bias 
Approximations in Experimental Uncertainly Analysis," AIAA 94-0772, 1994. 
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( 23 ) 


while Eq. (19) gives 


Bi = (B 2 f 1 + (B 2 f 2 +(B 2 f 3 + (B 2 / 4 
B 12 = (B,) 2 (B 2 ) 2 +(Bi) 3 (B 2 ) 3 (24) 
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Appendix 2 

Baseline Data Set, Data Reduction Data, 
and Uncertainty Estimates 


This appendix contains tables for the baseline data set (BDS), the data 
sets to use for data reduction and verification of the new model, and the 
associated uncertainty estimates. These data sets were generated by 
reviewing the test profiles of all Engine 3001 tests performed on the 
Technology Test Bed facility. The tests and time slices chosen and for 
inclusion in these sets have the same mixture ratios, inlet conditions, 
repressurization flowrates, and similar hardware configurations. 
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| Engine Control Conditions: 104% RPL j 







PTMCHB=PSMCHB? 

63 

MCC Pc Avg 

3126.25 









mixture ratio 

6.01 









02 repress flowrate 

1.10 









H@ repress flowrate 

0.20 







TTHTNK 

1021 

LPFP inlet temp lev 6B 

37.02 







TTVL18 

1058 

LPOP inlet temp lev 6A 

166.27 







PTHTNK 

9910 

LPFP INLET PR 

38.43 







PSVL18? 

9911 

LPOP INLET PR 

97.56 







ROC ETS variable 

PID# 

NAME 

avg 

S 


B 

DEI 

U 

OH 

PTFPRB=PSFPRB? 

8500 

FPB Pc (thru liner) 

5278.81 

74.86 

1.4 

105.58 

2.0 

183.20 

3.5 

PTOPRB=PSOPRB? 

8458 

OPB Pc (thru liner) 

5244.05 

11.21 

0.2 

104.88 

2.0 

107.25 

20 

PSPBSO 

341 


7302.45 

33.95 

0.5 

146.05 

2.0 

161.06 

22 

PSVL02 

86 


237.72 

4.14 

1.7 

4.75 

2.0 

9.55 

4.0 

PSVL19? 

209/210 avg 

HPOP inlet PR avg 

337.20 

14.97 

4.4 

6.74 

2.0 

30.69 

9.1 

SNFL 

754 


15668.87 

28.52 

0.2 

313.38 

2.0 

318.53 

20 

SNFH 

734 


5225.15 

15.91 

0.3 

104.50 

20 

109.24 

21 

TTHTFD 

231-232 AVG 


1801.09 

26.15 

1.5 

36.02 

2.0 

63.51 

3.5 

jjj 

233-234 AVG 


1345.57 

28.47 

21 

26.91 

20 

62.97 

4.7 

UJ 

659 


93.35 

0.22 

0.2 

1.87 

2.0 

1.92 

21 

7TVL12 

18 


429.56 

3.56 

0.8 

8.59 

20 

11.15 

26 

WFOPB? 

8805 


35.14 

0.46 

1.3 

0.70 

2,0 

1.15 

3.3 

WLPFT 

8801 

LPFT inlet flow (calc) 

29.51 

0.34 

1.2 

0.59 

20 

0.90 

3.1 

XRFPOV 

42 

FPOV act posA(FPVI) 

82.94 

1.29 

1.6 

0.83 

1.0 

272 

3.3 

XROPOV 

40 

OPOV act pos A (OPV1 ) 

67.23 

1.23 

1.8 

0.67 

1.0 

! 256 

3.8 


Table 1 Baseline Data Set for 104% rated power level. 
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RocetsVar 

PM# MoasurwTMnt Nan* 


Test 21 


Test 31 


9903 SLICE START TIME 


67 

68 

89 

90 

91 


54 

55 

56 

58 

58 


9904 SLICE END TIME 


88 

89 

90 

91 

92 


55 

56 

57 

58 

59 


9906 MCC PC CONTROL 


3126.25 

3126.25 

3126.25 

3126.25 

3126.25 


3126.25 

3126.25 

3126.25 

3126.25 

3126.25 


Data Reduction Manurtnonti 













PTFPRB 



5292.16 

5286.79 

5286.98 

5773.48 

5315.98 


5245.12 

5238.32 

5244.72 

5243 

5240.67 

PTOPRB 



5227.96 

5224.66 

527216 

5303.1 

5268.42 


5239 

5224.06 

5227.35 

5224.67 

■g~HI-.ll 




7344.22 

7336.53 

7341.87 

7356.28 

7335.17 


7291.94 

7279.96 

7290.5 

7289.63 

7290.13 




231.89 

231.41 

231.37 

23Z46 



239.66 

239.66 

236.93 

236.94 

239.36 


20GU210 *vg HPOP inlet PR 


342.295 

342.7 

343.375 

343.43 



344.51 

343.965 

344.945 

345.175 

344.97 

SNFL 

754 LPFTP Speed A 


15706.7 

15687.64 

1569Z53 

15710.43 

15681.86 


15660.63 

15658.16 

15650.69 

15650.94 

15652.05 

SNOl 

734 LPOTP speed NFD 








5207.08 

5204.37 

5210.57 

5206.67 

5209.42 

TTHTFO 



1770.755 

1767.335 

1765.095 

1764.82 

1763.485 


1814.205 

1810.11 

1807.675 

1810.45 

1806.91 

TTHTOD 



1313.36 

1312.705 

1319.19 

1317.66 

1321.58 


1343.405 

1349.415 

1346.855 

1344.93 

1341.765 

T7VL03 



83.72 

93.67 

83.65 

93.74 

93.63 


93.24 

93.22 

93.23 

93Z3 

93.21 

TTVL12 



433.39 

432.72 

43Z99 

433.26 

433.78 


*25 

425 

425 

425 

425 

WFOPB 



35.93 

35.68 

35.56 

3554 

35.68 


36.01 

35.02 

35.04 

34.98 

34.93 

WLPFT 



30.07 

30.01 

30.01 

30.05 

29.96 


295 

29.27 

29.28 

29.3 

29.31 




65.75 

85.74 

6559 

65.44 

65.35 


67.16 

67.14 

67.25 

67Z8 

67.34 

tXRFPOV 



91.1 

81.05 

81.07 

80.96 

61.04 


83.38 

83.42 

83.43 

63.29 

83.26 

1 Control Point M*w*mant* 





PTMCHB 

63 MCC Pc Avg 


3129.9 

3126.9 

3127.43 

3130.65 

3125.09 


3127.84 

3123.97 

3127.23 

3125.86 

3125.24 

W01 

9906 FACILITY FUEL FLOW 








154.57 

154.6 

154.57 

154.8 

154.58 

WLPFP 

9909 FACtUTY OXYGEN FLOW 








925.99 

925.66 

926.12 

924.73 

926.22 















T7HTNK 



36.98 

36.97 

36.99 

36.98 

36.98 


37.04 

37.04 

37.04 

37.04 

37.04 

nvue 



167.02 

167.03 

167.02 

167.02 

167.02 


164.67 

164.72 

164.69 

164.7 

164.72 

PTHTNK 



35.64 

35.67 

35.74 

35.64 

35.68 


39.3 

39.24 

39.27 

39.28 

39.36 

PSVL187 



97.63 

97.87 

97.77 

97.68 

97.95 


95.89 

95.81 

95.77 

95.67 

95.99 

1 AddXJonal MMSuramants 





TTVL02 

15 HPFP inlet temp avg 


4Z64 

4Z63 

4Z64 

42.65 

4Z64 


4Z75 

42.75 

42.74 

■KESZl 



638 LPFP OIS TEMP 


47.77 

47.76 

47.76 

47.79 

47.78 


43.62 

43.62 

43.63 

43.62 

43.63 


458 HPFP DISCH PR NFD 


6192.66 

6165.23 

6161.99 

6194.29 

6176.48 


6184.88 

6178.84 

6181.06 

6180.82 

6178.73 

PSVL21? 

6751 HPOP CHS VOITURI IN PR 


4158.62 

4153.15 

4154.17 

4157.48 

4150.93 


4123.63 

4118.52 

4123.85 

412Z26 

4123.45 

PSVL21? 

334 HPOP (fecharge Pr NFD 


4100.28 

4096.31 

4098.29 

4100.07 

4096.4 


4095.83 

4091.09 

4096.82 

4095.43 

4096 36 

PSVL13? 



3519.41 

3618.44 

3516.66 

3520.99 

3514.13 


3512.16 

3509.51 

3511.79 

3511.2 

3510.76 

PSVL21? 



3633.41 

361Z77 

3599.74 

3602.39 

3583.06 


3709.49 

3692.87 

3689.46 

3686.24 

3685.91 




4502.18 

4496.26 

4497.19 

450Z7 

4494.02 


448Z86 

4478.68 

4481.04 

4480.04 

*479.76 




432.78 

433.25 

433.23 

433.25 

433.66 


42851 

428.25 

428.58 

428.18 

428.18 

TTPBSO 



208.4 

208.37 

208.4 

208.46 

206.39 


207.17 

207.11 

207.17 

207.16 

207.18 

PSVL04 



5810.34 

5803.32 

5798.79 

5815.16 

5794.43 


5827.87 

5822.04 

5623.7 

5823 

5821.86 

PSVL12 



4458 92 

4454.12 

4455.92 

4461.32 

4448.7 


447Z38 

4467.55 

4469.98 

4469.38 

4466.34 

WCCV 



73.95 

73.06 

74.17 

73.41 

7357 


7159 

71.95 

72.17 

7Z12 

72.06 




14.26 

14.25 

14.21 

14.21 

1455 


14.33 

14.24 

14.35 

14.26 

14.21 




11.63 

11.79 

11.99 

11.7 

115 


1Z58 

12.62 

12.59 

1Z6 

12.56 




14.38 

14.33 

14.79 

1459 

14.61 







WF2 



40.47 

40.37 

40.99 

40.3 

40.66 







WLPOT 



190.14 

191.71 

191.68 

190.73 

189.8 


19Z33 

192.15 

192.32 

192.31 

192.25 

WHPOP-... 



897.49 

896.23 

894.75 

697.59 

697.96 


895.92 

895.26 

895.72 

895.33 

894.77 

WLPOV 



26.72 

26.79 

26.65 

26.62 

26.88 


26.69 

26.65 

26.65 

26.66 

26.7 

WFPOV 



76.06 

75.66 

76.66 

76.17 

75.1 


74.76 

74.65 

74.63 

74.68 

74.62 


Table 2. Data Reduction Data Sets, Tests TTB021 and TTB031, 104% RPL 
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Table 3. Data Reduction Data Sets, Tests TTB034 and TTB035, 104% SPL 


A2-4 























Rocets Var 

PM* 

Measurement Name 


9903 

SLICE START TIME 


9904 

SLICE END TIME 


9906 

MCC PC CONTROL 


Data Reduction Measurements 


B (%) U 


PSVL03 

PSVL21? 

PSVL21? 

PSVL13? 

PSVL21? 


PSVL04 

PSVL12 


HPFP inlet temp avg 
LPFP CHS TEMP 
HPFP DISCH PR NFD 
HPOP CHS VENTURI IN PR 
HPOP dscharge Pr NFD 
LPFT discharge Pr 
MINJ LOXINJ PR 2 
LPFT VENTURI PR 
LPFT VENTURI TEMP 
PSP discharge temp ch A 
MFV discharge Pr 
MCC coolant discharge Pr A 
CCV inlet flow (calc) 

Nozzle clnt inlet flow-1 (calc) 
Nozzle clnt inlet flow-2 (calc) 
Nozzle clnt inlet flow-3 (calc) 
TOTAL NOZZLE FLOW 
LPOT inlet flow (calc) 

HPOP dscharge flow (calc) 
OPB LOX flow (calc) 
FPBLOX flow (calc) 


15668.87 

5225.15 

1801.09 

1345.56 



Control Point Measurements 

PTMCHB 

63 

MCC PC Avg 

W01 

9908 

FACILITY FUEL FLOW 

WIPFP 

9909 

FACILITY OXYGEN FLOW 
mixture ratio 

02 REPRESS FLOWRATE 
H2 REPRESS FLOWRATE 

TTHTNK 

1021 

LPFP inlet temp lev 6B 

TTVL18 

1058 

LPOP inlet temp lev 6A 

FTHTNK 

9910 

LPFP INLET PR 

PSVL18? 

9911 

LPOP INLET PR 


Additional Measurements 



6342 

20 

3.10 

20 

18.57 

20 

0.74 

20 

3.84 

23 

3.36 

8.7 

36.27 

41.5 


Table 4. Data Redaction Data Sets, Average Values and Uncertainties, 104% RPL. 
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ROCETS vAR 


PTFPRB 

PTOPRB 

PTPBSO 

PTVL02 

PTVL20 

SNFL 

SNOL 

TTHTFD 

TTHTOO 

TTVL03 

TTVL12 

WFOPB 



PSVL03 

PSVL21? 

PSVL21? 

PSVL13? 

PSVL21? 


TTPBSO 

PSVL04 

PSVL12 

WCCV 


PW # Measurement Name 

SLICE START TIME 
SUCE END TIME 
MCC PC CONTROL 

Data Reduction Measurements 


8500 FPB Pc (thru liner) 

8458 OPB Pc (thru liner) 

341 PBP discharge Pr NFD 
86 HPFP inlet Pravg(LFDP) 
209/210 avg HPOP Mat Pr avg 
754 LPFTP Speed A 
734 LPOTP speed NFD 
231/232 avg HPFT discharge temp avg 
233/234 avg HPOT discharge temp avg 
659 HPFP discharge temp 
18 MCC dnt discharge temp B 

8805 OPB fuel flow (calc) 

LPFT inlet flow (calc) 
OPOV act pos A (OPV1) 
FPOV act pos A (FPV1) 


Control Point Measurements 


63 | MCC Pc Avg 

FACILITY FUEL FLOW 

9909 FACILITY OXYGEN FLOW 
mixture ratio 

02 REPRESS FLOWRATE 
H2 REPRESS FLOWRATE 
1021 LPFP inlet temp lev 68 

1058 LPOP inlet temp lev 6A 

9910 LPFP INLET PR 

9911 LPOP INLET PR 




7712.7 

247.9 

335.9 
161582 

5381.6 

1851.6 
1420.2 

96.1 
424.0 

36.1 
312 

72.1 
86.3 


5.9 

03 

22.4 

1.6 


15 

HPFP inlet temp avg 

639 

LPFP DIS TEMP 

459 

HPFP DISCH PR NFD 

8751 

HPOP DIS VENTURI IN PR 

334 

HPOP discharge Pr NFD 

437 

LPFT discharge Pr 

8684 

MINJLOXINJPR2 

8014 

LPFT VENTURI PR 

8016 

LPFT VENTURI TEMP 

93 

PBP discharge temp ch A 

356 

MFV discharge Pr 

17 

MCC cooiarrt discharge Pr A 

8818 

CCV inlet flow (calc) 

8815 

Nozzle dnt inlet flow-l (calc) 

8816 

Nozzle dnt inlet flow-2 (calc) 

8817 

Nozzle dnt inlet flow-3 (calc) 
TOTAL NOZZLE FLOW 

8802 

LPOT inlet flow (calc) 

8819 

HPOP discharge flow (calc) 

8804 

OPB LOX flow (calc) 

8810 

FPB LOX flow (calc) 



6.7 

2.0 

161.6 

1.0 

53.8 

1.0 

37.0 

2.0 

28.4 

2.0 

1.9 

2.0 

8.5 

2.0 




0.0 

0.0 

0.1 

02 

9.1 

0.1 

6.6 

0.2 

4.2 

0.1 

33 

0.1 

27.5 

0.7 

4.7 

0.1 

1.9 

0.4 

1.6 

0.8 

6.5 

0.1 

8 2 

02 

0 2 

0.3 

0.3 

1.9 

0 2 

1.8 

0.1 

0.3 

0.1 

0.1 

03 

0.2 

13.7 

1.4 

0.5 

1.8 


u 

U(%) 

84.1 

1.5 

892 

1.6 

157.9 

2.0 

5.6 

2.3 

42.4 

12.6 

170.1 

1.1 

61.9 

12 

38.9 

2.1 

53.0 

3.7 

1 3 

2.0 

9.1 

2.1 

0.9 

2.5 

0.7 


1.6 


0.7 

0.8 



Table 5. Data Reduction Data Sets, Average Values and Uncertainties, 109% RPL. 
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Measurement Name 

9903 

SLICE START TIME 

9904 

SUCE END TIME 

9906 

MCC PC CONTROL 

1 Data Reduction Measurements 



Control Point Measurements 

PTMCHB 

63 

MCC Ft Avg 

W01 

9908 

FACILITY FUEL FLOW 

WLPFP 

9909 

FACILITY OXYGEN FLOW 
mixture ratio 

02 REPRESS FLOWRATE 
H2 REPRESS FLCWRATE 

TTHTNK 

1021 

LPFP inlet temp lev 68 

TTVL18 

1058 

LPOP inlet temp lev 6A 

PTHTNK 

9910 

LPFP INLET PR 

PSVL18? 

9911 

LFOPINUET PR 


Additional Measurements 


PSVLD3 

PSVL21? 

PSVL21? 

PSVL13? 

PSVL21? 


HPFP inlet temp avg 
LPFP DIS TEMP 
HPFPDiSCH PRNFD 
HPOP DIS VENTURI IN PR 
HPOP discharge Pr NFD 
LPFT dsehargePr 
MINJLOXINJPR2 
LPFT VENTURI PR 
LPFT VENTURI TEMP 
PSP discharge temp ch A 
MFV discharge Pr 
MCC coolant discharge Pr A 
CCV inlet ftow (calc) 

Nozzle dr* inlet flow-1 (calc) 
Nozzle cint inlet flow-2 (calc) 
Nozzle clnt inlet flcw-3 (calc) 
TOTAL NOZZLE FLOW 
LPOT inlet flow (calc) 

HPOP discharge ftow (calc) 
OF© LOX ftow (calc) 

FPB LOX flow (calc) 


4966.36 

4963.12 

6970.15 

230.76 

342.61 

15228.81 

5096.70 

1766.23 

1296.59 

90.99 

431.60 

3409 

28.11 

6454 

80.85 


74,79533 

74.44679 


0.590342 

0.322703 

0.404248 


85.55 

1.7 

78.04 

1.6 

162 45 

2.3 

8.27 

3.6 

7.66 

22 

167.69 

1.1 

53.29 

1.0 

61.65 

3.5 

35 48 

2.7 

1.86 

2.0 

11.37 

2.6 

0,96 

2.8 

0.83 

3.0 

1.07 

1.7 

2.04 

25 


60.86 

20 

3.00 

20 

18.00 

20 

0.74 

20 

3.85 

23 

3.14 

8.1 

275 

28 


4256 

0.05 

0.1 

200 

0.9 

0.86 

20 

44.48 

1.90 

4.3 

200 

0.9 

3.90 

8.8 

5873.75 

5.50 

0.1 

200 

117.5 

117.99 

20 

3946.85 

1106 

0.3 

200 

78.9 

81.99 

21 

3911.84 

5.09 

0.1 

200 

78.2 

78.90 

20 

3363.04 

280 

0.1 

200 

67.3 

67.49 

20 

3489.10 

80.16 

23 

200 

69.8 

174.85 

5.0 

4280.82 

7.32 

0.2 

200 

85.6 

86.86 

2.0 

43280 

240 

0.6 

200 

8.7 

9.90 

23 

206.42 

1.12 

0.5 

200 

4.1 

4.69 

2.3 

5530.18 

9.40 

0.2 

200 

110.6 

112.19 

20 

4256.53 

6.74 

0.2 

200 

85.1 

86.19 

20 

70.17 

0.60 

0.9 

3.00 

210501 

2.42 

3.5 

13.22 

0.41 

3.1 

210 

0.277526 

0.86 

6.5 

11.42 

0.27 

2.4 

210 

0.239862 

0.59 

5.2 

13.65 

0.14 

1.0 

210 

0.286706 

0.40 

29 

38.13 

0.54 

1.4 

210 

0.800632 

134 

3.5 

187.61 

0.84 

0.4 

2.10 

3,939884 

4.28 

2.3 

870.76 

1224 

1.4 

210 

18.28594 

30.56 

3.5 

25.14 

0.66 

2.7 

210 

0.5 

1.45 

5.8 

70.777 

0.384052 

0.542622 

- 

> 1.41554 

1610508 

2.275468 


Table 6. Data Reduction Data Sets, Average Values and Uncertainties, 100% RPL. 
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Appendix 3 

Uncertainty in Linear Regression 


Introduction 

When experimental information is used in a model in the form of a 
curvefit an uncertainty due to the original experimental program is 
introduced. Currently there is no accepted methodology to assess this 
uncertainty. The work reported here is a continuation of a previous effort 1 to 
develop a methodology to assess the uncertainty in linear regression 
analysis, or curvefits, when the experimental uncertainties contain precision, 
systematic, and correlated systematic uncertainties.. The work reported in 
reference [1] developed a methodology to assess the uncertainty in the 
regression coefficients, slope and y-intercept, for a 1 st order linear regression. 
The most common forms of regression models used are straight lines and 
polynomial curvefits where both the X and Y axes are functional 
relationships containing uncertain experimental information. While 
multivariate regressions are also important models, they were not studied in 
this effort. 

In using the basic regression models with typical engineering data 
many of the underlying assumptions used in the statistical development of 
the regression models are violated. The underlying assumptions assume that 
the errors “ (1) are unbiased; (2) have constant variance; (3) are uncorrelated, 
and (4) are normally distributed.” 2 In a typical experimental program the 
first three assumptions will usually be violated. Seber 2 and Montgomery 3 
discuss these assumptions and provide methods to deal with some of the 
violations, however their methods are mathematically rigorous and do not 
readily lend themselves to the propagation of experimental uncertainties. 
Violation of these assumptions does not preclude using the regression models 
to obtain a model for the data, but the model must be used with the 
realization that it only represents the “best fit” of a curve through the data. 
Any confidence interval then used with that curvefit must properly account 
for any violation of the underlying assumptions. 

The methodology presented in this report is an extention of accepted 
uncertainty propagation techniques to the linear regression equations. The 
effectiveness of this new methodology was evaluated using Monte Carlo-type 
simulations. This report will briefly discuss the uncertainty analysis 
techniques and then present its application to linear regression analysis, first 
for 1 st order (straight line) regressions and then for higher order (polynomial) 
regressions followed by a summary of the Monte Carlo simulation technique 
employed to evaluate the methodology. The report concludes with the 
application of the curvefit uncertainty methodology to compressor map 
characteristics. 

Regression Analysis Uncertainty Methodology for Slone and v-Intercept 
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The uncertainties in the slope, m, and the Y-intercept, c, are functions of 
the uncertainties in the determinations of the X and Y variables. The values 
of m and c are obtained by minimizing the sum of the squares of the 
deviations between the line and the data points, commonly known as the 
method of least squares. The development of the equations to calculate m 
and c can be found in statistics books 2 ’ 3 9 and only the equations will be 
presented here. For N (X, Yi) data pairs, the slope of the line, m, is 
determined from 


N 


m = 


i=l 


i— 1 i=l 


N f N N 

V 1=1 > 


i=l 


and the intercept, c, is determined from 

N N N N 


C = 


i* 1 


i=l 


r=l t=I 


N f N 


1=1 


i=l 


( 1 ) 


( 2 ) 


In this paper we present a methodology to determine the uncertainties in 
linear regression coefficients. The effectiveness of this methodology was 
analyzed using a Monte Carlo-type simulation assuming the true 
relationship between the X and Y variables is in fact linear. 

The approach presented in this paper is an application of the 
uncertainty analysis methodology presented above to the regression 
equations for slope and intercept 


Considering Eq.s (2) and (3) to be data reduction equations of the form 

m=m(X i ,Xi,...,X N ,Y i ,Y 1 ,...,Y N ) (3) 


and 


c=c{X x ,X 1 ,...,X n XY 2 ,...J n ) (4) 

and applying the uncertainty analysis equations, Eq.s (5)-(13), the most 
general form of the expression for the uncertainty in the slope of the line, m, 
is 
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(5) 




-Zl^l *WZ I 


dm 

dX 


dm 

dY, 

r 


P 2 
r r f 


dm 

VdX~ t 


+ 2ZZ 


7 ' f dm' 


i= 1 * = 1 


dm 


B 


x t r k 


\dxj\dY,, 

where B r< is the systematic uncertainty for the Y t variable, B x< is the 
systematic uncertainty for the X t variable, B rji is the covariance estimator 
for the correlated bias uncertainties in the Y x and Y k variables, B XfXt is the 
covariance estimator for correlated bias uncertainties in the X i and X k 
variables, P Tj is the random uncertainty for the Y i variable and P x is the 
random uncertainty for the X t variable, and B XJi is the covariance estimator 
for the correlated bias uncertainties between X i and Y t . 

A similar expression for the uncertainty in the intercept is 


J-l J 


f If dc 




J 7 dc' 


V dY t j 
f 


dc 


dXj 


P 2 

r x. 


(6) 


j j ( 


2 S z 


dc 


dc 


ft ft \dX t )\dY k 

The partial derivatives are 


B X t Y k 


dm 


NX.-^X, 


(7) 


r N { ” Y 

' "2KHZ* 

/=! V t=\ / 


dc 


— '=1 


tffl-XiZx, 

/si 


7 N { N Y 

1 Z* 

i=\ ^ /=! J 


( 8 ) 


NY - * 
dm i 

V 

-l 


( N N N ' 

> /-I i-1 t-l 


/-i y 

flX N 

' *E(*, S )-| 


2 

( N ( 

: 

^ i«l 'j 

M \ 2 

L*,) 

r«i / 

r 


( 9 ) 


A3 -3 



and 


dc 

ax, 


2 X.fl -i x.r.-r.f x. [f(x?)fr, -fx.fx.rXiNx, - 22 ;*,) 

i-i i-l 1-1 ^i-i 1-1 /-i 1-1 S\ i-i / 






( 10 ) 


Equations above show the most general form of the equations for the 
uncertainty in the slope and the y-intercept, allowing for correlation of bias 
errors among the different Xs, among the different Fs and also among the 
Xs and Fs.. If none of the systematic error sources are common between the 
X variables and the F variables, the last term of the equations, the X-Y 
covariance estimator, is zero and the equations reduce to 


Ul = 't( 3n ) B 2 + 22 +f t (* H ) 2 p T 3 

m tfUrJ r ‘ is & XarJ larj rjt iilorj r - 

rAdxJ x > ^ttAdx.Kdxj x ' x > j?Adx 


(id 


and 


Ul = I 


i=] 


dc 


dY, 


j - i j 


iJ 


K + 2 Z I 


i—i Jt=i+: 


dc 






j-\ 


B 


x t 


ill 

1=1 Jfe=i +1 





( 12 ) 


Regression Analysis Methodology for Predicted Value; Straight Line 

The expression for the uncertainty in the predicted value, Y, at a given 
X variable 


7(X) = mX+c 


(13) 


(where X is assumed to have no uncertainty) is determined in a similar 
manner by applying Eq.s (2) through (8) to equation (21) to obtain 


i=l \0*tJ i= 1 Jk=i+1\^, 


(dY^(dY\ B 

YJk Z-J 

i=l 


+ 


m 


/-i j 


Y| “ I B 2 x+ 2 yy(^ 


.ar* 




dY. 


B x,x t + 4 


i=i 


iJ 

dY 

dX, 


P 2 

r Y, 


+ 2 


dY 
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and 


(16) 


dY 3n v dc 

cY, m, 

and where the partial derivatives for m and c are as defined in Eq.s (15) 
through (18). The partial derivatives can also be determined numerically, as 
will be discussed later. 


Calibration Problem 

If the value of X being used in the regression model is measured with 
the same apparatus as the Xi data, and thus shares the same systematic 
error sources the uncertainty expression must be expanded to include these 
uncertainties and correlations. This problem is often encountered in 
calibration processes and the uncertainty expression is 



If no errors sources are common between the X and Yi variables, the X-Yi 
correlation terms are omitted, these terms comprise the last line of the above 
expression. 


Numerical Partial Derivative Determination 

The partial derivatives can be approximated directly using a finite 
difference numerical technique, often referred to as a “jitter-routine.” [4] A 
simple forward- difference partial derivative, which has accuracy of order h is 

/(A 1 +/»,Y 2 ,...,Jr B )-/(A;,A 2 ,...,Y > ,) 

dX x h y J 


It was determined during this effort that round-off errors could be 
si gnifi cant, so double precision arithmetic was used. Accuracy was also 
enhanced by using a second order accurate central-differencing scheme 9 , such 
as: 


/(X, + h,X at ...,X m )-f(X l -A,Y 2 ,..., X„) 
dX l 2 h 


( 19 ) 


Traditional Statistical Regression Uncertainties: Straight Tine 
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To demonstrate the value of the uncertainty propagation methodology 
presented above a comparison to the confidence intervals provided by 
traditional statistical techniques is needed. The form of the statistical 
intervals used in this study was taken from Natrella 10 , more discussion of 
these intervals and additional forms of the intervals can also be found in 
references [2] and [3]. The variance of the slope is given as 



( 20 ) 


with the 95% confidence interval represented by 

m±Xs m (21) 


where the 95% coverage factor, X, is determined using the two-tailed 
Student’s t-distribution, with n-2 degrees of freedom, as 

* = *.«(»- 2) (22) 


Similarly, the variance for the y-intercept is given by 


1 X 2 


and the 95% confidence interval is determined as 


(23) 


c±As c 


(24) 


A confidence interval about any point along the regression line, what has 
been referred to as the predicted value is obtained with the interval Wi as 


W i =Xs I \ 


i (*--*)’ 


— + 
n 


Yi 


(25) 


and the 95% confidence region is expressed as 

y t ±W t (26) 

Where the sum of squares statistics used in the determination of these 
intervals are defined as 


and 


. (z* 

I n 


$»= 'EX?-!!*—! 

m n 

i= 1 

(27) 

( » ^ 

2 


• M 

n 


« n 

i=l 

(28) 
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and 


(29) 


and 




T.x.U, 


1=1 1=1 

n 


=i(x,-m-y i 



(30) 


X and Y as Functional Relations 

In many, if not most instances, the test data will be expressed in some 
form of functional relationship. In these cases, the data will not be the 
variables used in the regression, instead the result from the functional 
relationships will be used. Examples of common functional relations are 
Reynolds number, flow coefficient, pressure coefficient, specific fuel 
consumption, etc. This can be expressed for a general example as 

X i =f(VARl l ,VAjR2 1 ) (31) 

and Y t = f(VARl : , VAR4 t , VAR5 i ) (32) 


so the expression for the uncertainty in the predicted value becomes 


n 5 


<?= zi 


dY 


ZftVdVARj^ 

n-1 n 5 f 

*ZZZ 


[ VAR h 


dY 


n 5 

'sfw ySVASj,. 


B 


VARU 




dY 


f 


dY 


\B V 


n n 4 5 

+2 ZSE S 


dY 


VdVARjJ VARj ‘ VARJ * 
dY 


B v 


i=jti V dVARj i ) v dVARl k J yARj '™ 1 ' 


(33) 


where the second term accounts for correlated systematic uncertainty sources 
within each variable and the third term accounts for systematic uncertainty 
sources common between variables. Similar expressions for the uncertainty 
in the slope and y-intercept are readily obtained. The complexity introduced 
by the transformation necessitates determining the partial derivatives 
numerically. 

Another important form of transformation is when one or both of the 
variables are represented in logarithmic space. A similar technique is used 
whereby the functional relationship as expressed in Eq.s (39) and (40) 
become the logarithmic function desired, such as 

X i =\n(yABl i ) and Y^XniVARl,) (34) 
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The partial derivatives then account for the transformation and the 
uncertainties in the measured variables are propagated in the units of the 
measurand without transformation. 

Methodology Development and Monte Carlo Simulations 

The methodology for assessing the uncertainty in linear regressions 
was developed using Monte Carlo-type simulations. Monte Carlo-type 
simulations are often used in uncertainty analysis to determine the 
effectiveness of a particular uncertainty model. For this work, what is 
referred to as a Monte Carlo-type simulation simply means generating 
numbers to represent experimental data with some amount of error randomly 
obtained from predefined error distribution populations. 

The Monte Carlo simulations were conducted in the following manner. 
“True” values for data from a relationship with specified coefficients were 
determined. The word true is emphasized to indicate that it represents the 
actual physical quantity of the parameter if it could be measured without any 
bias error or precision error, which is always an unobtainable value. The 
two-sigma (2 standard deviation or 95% confidence) bias limits and precision 
limits for each variable were then specified. The errors in each variable were 
assumed to come from these normally distributed error populations with the 
specified standard deviations. A random value for each bias error and 
precision error was found from a Gaussian random deviate generator 
subroutine using the specified standard deviations. The Gaussian deviates 
have a mean of zero and an equal probability of being positive or negative. 
The bias errors within each variable were assumed to be from the same 
source, and were assigned the same random deviate from the Gaussian 
distribution. In some cases these bias errors were a fixed amount and in 
some cases they were a percent of reading type error. Precision errors were 
obtained by sampling the precision error populations repeatedly to obtain 
independent random deviates for each variable. The individual error values, 
bias errors and precision errors were then summed and added to the true 
value to obtain a data point with errors from the specified error populations. 
These data points were then used in the linear least squares equations to 
obtain the value of the regression coefficients and the curvefit model. These 
coefficients and the curvefit model represent the “best fit” of the experimental 
data when the bias and precision errors are present. 

A 95% confidence uncertainty interval for the result was calculated 
from the uncertainty propagation analysis equations for m, c, and Y. A ±Um 
interval was placed around the slope coefficient value, m, and if the true 
value of the slope was found to be within the interval a counter was 
incremented. A similar procedure was used with the y-intercept and the 
predicted value from the curvefit expression. This procedure was repeated 
10,000 times and the percent coverage, or number of times the true result 
was within the estimated interval, was determined. Using this procedure, 
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the effectiveness of the uncertainty propagation equation could be 
investigated by checking whether or not the true value is within the 95% 
confidence uncertainty interval about the measured result 95% of the time, 
we can determine how well our uncertainty approximation equations work. 

In the previous work 1 the precision uncertainties used in the 
propagation equations were determined using the “large sample 
approximation” and calculating a standard at each setpoint. The effect of 
sample size, i.e. the number of setpoints and the number of data taken at 
each setpoint, was thus studied. The conclusion drawn from this analysis 
was that when the “correct” precision limit is used in the propagation 
equations the desired confidence interval is obtained. Thus, all the work in 
this report will assume that the correct precision uncertainties have been 
estimated based upon prior information and experience with the test rig so 
that the simulations will evaluate the methodology only. 

A useful statistic from the simulation is the uncertainty ratio, the ratio 
of the average uncertainty intervals for the regression coefficients from the 
10,000 iterations divided by the true 95% confidence intervals. The true 95% 
uncertainties are calculated as twice the sample standard deviations, Sm, Sc, 
and Sy, from the 10,000 samples of the regression coefficients {not to be 
confused with the standard deviations of the curvefit, s m , Sc, and W, as 
defined in Eq.s (28), (31), and (33)}. The sample standard deviations from 
the 10,000 sample population can be expected to be good representations of 
the actual standard deviations of the infinite population with the elemental 
uncertainty sources as defined. An uncertainty ratio near unity shows that 
the uncertainty methodology works for the particular case, with values 
greater than one meaning an over-prediction and values less than one 
meaning an under-preduction. 


1 st Order Regression Uncertainty Simulation and Results 


Case 1: Dominant Systematic Uncertainties 
Bxi= 2.5 Bx2— 2.5% Bxs~By3“ 2.5 

Byi= 2.5 Byj= 2.5 

Pr= 0.0 Py= 0.1 


1 

Uncertainty Propagation 

Classical Statistics j 

ra 

Xi 

Yt 

UY(«f) 

Cov (%) 

Ratio 


Cov (%) 

Ratio 


80.0 

90.0 

4.58 

95.4% 

0.99 


3.9 

0.02 


90.0 

100.0 

4.81 

95.3% 

0.99 


2.5 

0.02 


100.0 

110.0 

5.06 

95.4% 

0.99 


1.8 

0.01 


110.0 

120.0 

5.32 

95.3% 

0.99 

0.11 

2.3 

0.02 


120.0 

130.0 

5.59 

95.4% 

0.99 

0.15 

3.0 

0.02 

I Coefficient Uncertainties 

Um 

Cov (%) 

Ratio 

ts 

Cov (%) 

Ratio 


true slope, m=1.0 

0.036 

95.6 

1.01 

0.005 

21.1% 

■m 

I true y-intercept, c=10.0 

3.56 

95.3 

0.99 

0.47 

20.2% 

— 


Table 1. 1st Order Regression Simulation and Results; Dominant Systematic 
Uncertainties. 
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Case 2: Comparable Uncertainties 
Bn= 2.5 Bx2= 0.0% Bx3=By3= 0.0 

Byi= 2.5 Byi= 0.0% 

Px =2.5 Py =2.5 





Uncertainty Propai 

ration 

Classical Statistics I 


Xi 

Yi 


Cov(%) 

Ratio 


Cov (%) 

Ratio 


80.0 

90.0 

4.47 

95.3 

0.99 

3.42 

85.3 

0.89 


90.0 

100.0 


95.4 

0.99 

2.41 

77.6 

0.70 


100.0 

110.0 

3.87 

95.4 

0.99 

1.98 

71.4 

0.60 


110.0 

120.0 


95.4 

0.99 

2.42 

77.4 

0.70 

16 

120.0 

130.0 

4.47 

95.6 

0.99 

3.43 

84.9 

0.89 

| Coefficient Uncertainties 

Um(t*f) 

dm 

Ratio 


Cov (%) 

Ratio 


true slope, m=1.0 

0.112 

95.1 

1.00 

0.164 

94.4 

1.464 

| true y-interoept, c=10.0 

11.84 

95.0 

1.00 

16.5 

94.0 

1.394 


Table 2. 1st Order Regression Simulation and Results; Comparable Systematic 
and Precision Uncertainties. 


Case 3: Dominant Precision Uncertainties 
Bn= 0.0 Bi2= 0.0 Bx3=By8= 0.0 

Byi= 0.0 Bys= 0.0 

Px =2.5 Py =2.5 




Uncertainty Propai 

fl 

0 

I 

Classical Statistics 


Xi 

Yi 


Cov(%) 

Ratio 

IE9HH 


Ratio 


80.0 

90.0 

2.74 

95.6 

1.00 

3.15 

95.0 

1.46 


90.0 

100.0 

1.96 

95.3 

1.00 

2.25 

94.8 

1.46 


100.0 

110.0 

1.59 

95.3 

1.00 

1.84 

94.6 

1.46 


110.0 

120.0 

1.98 

95.4 

1.00 

2.28 

94.9 

1.46 

\LJ 

120.0 

130.0 

2.76 

95.4 

1.00 

3.17 

94.8 

1.46 

1 Coefficient Uncertainties 

Un(«rf) 

Cov (%) 

Ratio 

tS terri 

Cov (%) 

Ratio 



0.112 

95.4 

1.01 

0.163 

94.94 

1.48 


11.31 

95.5 

1.01 

16.5 

94.79 

1.47 


Table 3. 1st Order Regression Simulation and Results; Dominant Precision 
Uncertainties. 


Calibration Problem Results 

The calibration problem was simulated, wherein the systematic 
uncertainty associated with the X variable in the curvefit expression was 
correlated with the systematic uncertainty in the X data used in the curvefit. 
Equation (18) was used with the X-Y correlation terms omitted. The same 
true linear expression was used as in the previous simulations, however the 
curvefit was evaluated at calibration points between the curvefit data. The 
input data, uncertainty information, and simulation results are shown in 
Table 4. 
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Calibration Problem 

Case 1: Dominant Systematic Uncertainties 
Bx= 2.5 Px= 2.5 

By = 2.5 Py= 2.5 

l 

Uncertainty Propaj 

Ration 

wm 

Xi/Xi(cal) 

Y, 

Uy 

Cov (%) 

Ratio 


80.0/85.0 

90.0 

4.24 

■ 

1.00 


90.0/ 95.0 

100.0 

3.91 

■ 

0.99 

IS 

100.0/ 105.0 

110.0 

3.87 

95.0 

1.00 

IX 

110.0/115.0 

120.0 

4.12 

95.0 

0.99 

Efl 

120.0/ 125.0 

130.0 

4.68 

95.1 

0.99 


Table 4. Calibration Problem Simulation and Results; Comparable Systematic and Precision 
Uncertainties. 


Table 4 again demonstrates that the appropriate uncertainty interval is 
obtained with the proposed methodology. The regression coefficient results 
and the traditional statistics results are similar to Tables 1, 2, and 3. They 
have been omitted because they are not affected by the calibration problem. 

Polynomial Regr ession Uncertainty Simulation Results 

A 2 nd order polynomial curve was used to generate “true” (X,Y) data 
pairs and the uncertainty propagation methodology was evaluated using the 
Monte Carlo simulation technique. The methodology provided the 
appropriate coverage for the curves simulated. Several different curves have 
been modeled and simulated, and the methodology appears to work in most 
curvefit applications. A few areas need further work before the methodology 
can be fully recommended, these areas are discussed in the next section. 

In lieu of presenting the results of the polynomial simulations 
performed, a polynomial curvefit for the High Pressure Fuel Turbopump 
turbine efficiency map will be discussed in Appendix 4. 

Future Research 

A number of areas of continued research are identified as a result of 
this effort In all of the Monte Carlo simulations performed, the precision 
limits used in the propagation equations were the exact (2a) values. Work 
similar to that presented in reference [1] needs to be performed and extended 
to study the effect of sample size, the large sample approximation, data 
location, the use of prior information, and the use of other precision limi t 
estimates. The polynomial curvefit results indicated that the methodology 
provides the appropriate uncertainty intervals in a well defined and well 
behaved problem. However, if the data and the curvefit using that data do 
not represent the “true” phenomena being investigated incorrect uncertainty 
intervals will be obtained. This could happen due to over-fitting or under- 
fitting the data. In general, the limitations of using this method to assess the 
uncertainties in polynomial curvefits must be investigated, and any 
underlying assumptions identified. 
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Since the primary motivation of this work is the assessment of 
uncertainties in propulsion system models, and since multivariate regression 
models are often used in engine system models a means of assessing the 
associated uncertainties is needed. Thus extension of this methodology to 
multivariate regressions is necessary. 

Another important area for future research is the engineering 
application of this methodology. For example, how should it be integrated 
into data acquisition and test data reduction systems. How should the 
uncertainty information be documented so that others reading a report have 
the necessary information. And finally, an algorithm or flowchart describing 
how to implement this methodology would be very useful. 

Conclusion 

The work presented in this appendix provides a new methodology for 
assessing the uncertainty associated with curvefits of experimental 
information. This methodology is based upon the propagation of 
experimental uncertainties using accepted techniques. In the analysis it was 
shown that the methodology provided the appropriate uncertainty intervals 
when correlated systematic uncertainties are present, and that traditional 
statistical confidence intervals do not provide the appropriate intervals in 
these situations. The methodology was shown to extend to cases where the 
curvefit variables are functional relations of experimental data. 
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Appendix - 4 

Monte Carlo Simulation of HPFT Efficiency Map 

Appendix 3 introduced the proposed methodology to assess the 
uncertainty of n th order linear regressions when the regression parameters 
are functional relations. To verify the methodology was applicable to the 
type of component performance maps used in the SSME model, a Monte 
Carlo simulation of the turbine efficiency map for the High Pressure Fuel 
Turbopump Turbine (HPFT) was conducted. This simulation is also 
important as another demonstration of the regression uncertainty 
methodology. 

As discussed in Section 3.3.1, the uncertainty to associate with a value 
from the HPFT efficiency parameter versus speed parameter map is being 
estimated from the cold-flow testing of the HPFT in the MSFC air flow 
Turbine Test Equipment facility 1 and the adjustments to engine operating 
conditions. 

The information used to obtain the HPFT turbine efficiency map was 
obtained from a July 1991, Rocketdyne Internal Letter 2 . This memo only 
contains the reduced data, and not the measured inlet and exit total 
pressures, total temperatures, and mass flowrates. Using information from 
this memo, and discussions with MSFC personnel involved in the HPFTP 
model testing 3 , information to set up the Monte Carlo simulation was 
generated. 

For the purposes of the Monte Carlo simulation a fourth order linear 
regression was obtained from the data in the Rocketdyne memo. The fourth 
order curvefit used to model the efficiency map was 

eta„ -<un(%J ' + 6.609 (i) 

Information about the inlet test conditions and the test profile from Ms. 
Hudson was used to generate a set of data which yields values on this curve. 
This information is shown in Table A-4.1. She also stated that the overall 
uncertainty associated with the temperature measurements is 1.0°R, 
consisting mostly of bias uncertainty, the overall uncertainty in the pressure 


1 Hudson, Susan T., Gaddis, Stephen W., Johnson, P., Dean, and Boynton, James, L., “Cold Flow Testing 
of the Space Shuttle Main Engine High Pressure Fuel Turbine Model,” AIAA Technical Paper 91-2503, 
AIAA/SAE/ASME/ASEE 27th Joint Propulsion Conference, June 24-26, 1991, Sacramento, CA. 

2 Boynton, J.,and Daumann, A., “Revision of SSME HPFTP Turbine Performance Based on NASA/MSFC 
Baseline Air Test Results,” July 23, 1991. 

3 Personal Communication with Ms. Susan Hudson, Experimental Branch, Structures and Dynamics Lab, 
NASA- Marshall Space Flight Center. 
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measurements is 0.25 psia, consisting mostly of precision uncertainty, and 
the overall uncertainty in speed measurement is approximately 15 rpm, 
consisting mostly of precision uncertainty. Using this information the 
uncertainties to use in the Monte Carlo simulation were chosen. 


| T01 

T02 

P01 

P02 

N 

PR 

U 

u/co 

ETA 

ETA/U/C0 



100 

62.5 

2000 

1.6 


0.096 

0.440 

4.562 

549.69 

524.96 

100 

71.4 

2000 

1.4 

87.870 

0.113 

0.491 

4.346 

549.69 

533.25 

100 

83.3 

2000 

1.2 

87.870 

0.152 

0.589 

3.883 

549.69 

508.37 

100 

66.7 

4000 

1.5 

175.739 


0.687 

3.324 

549.69 

499.32 

100 

62.5 

5000 

1.6 

219.674 

0.241 

0.729 

3.025 

549.69 

504.84 

100 

66.7 

5000 

1.5 

219.674 

0.258 

0.746 

2.887 

549.69 

496.52 

100 

62.5 

6000 

1.6 

263.609 

0.289 

0.770 

2.661 

549.69 

502.67 

100 

66.7 

6000 

1.5 

263.609 

0.310 

0.782 

2.522 

549.69 

494.85 

100 

62.5 

6982 

1.6 

306.753 

0.337 

0.794 

2.359 

549.69 

501.48 

100 

66.7 

6982 

1.5 

306.753 

0.361 

0.802 

2.222 

549.69 

517.67 

100 

76.9 

6000 

1.3 

263.609 

0.382 

0.807 

2.113 

549.69 

508.95 

100 

71.4 

6982 

1.4 

306.753 

0.394 

0.809 

2.051 

549.69 

517.47 

100 

76.9 

6982 

1.3 

306.753 

0.444 

0.812 

1.828 

549.69 

527.06 

100 

83.3 

6000 

1.2 

263.609 

0.455 

0.811 

1.782 

549.69 

517.74 

100 

76.9 

8000 

1.3 

351.478 

0.509 

0.805 

1.582 

549.69 

527.36 

100 

83.3 

6982 

1.2 

306.753 

0.530 

0.800 

1.511 

549.69 

528.15 

100 

83.3 

8000 

1.2 

351.478 

0.607 

0.772 

1.272 

549.69 

531.65 

100 

83.3 

10000 

1.2 

439.348 

0.759 

0.647 

0.852 


Table A-4.1. Input “Data” used in HPFT Efficiency Map Monte Carlo Simulation 


The x-axis of the turbine efficiency map is the velocity ratio, defined as 
the disk tangential speed, U, divided by the isentropic spouting velocity, Co. 
Where the velocity ratio is determined from 


U = ^ pitch X ^ ft/ 

229.18 /sec 


( 2 ) 


and the isentropic spouting velocity is found from 

C 0 = ^2g c J AH T _ T(isen) y/ ec 
and the delta-enthalpy is calculated as 


AH 


T-T(isen) 



(3) 


(4) 


where gc is the gravitational constant, J is the Joules Constant, and the 
229.18 in the denominator of the U is a unit conversion constant. The data 
reduction equation for the turbine model efficiency is 
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T - T 
^01 l 02 


T 1 - f ^ 02 / 1 7 

1 oi 1 l /p J 
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However, the map efficiency, the y-axis, is plotted in terms of the efficiency 
divided by the velocity ratio, 


Using these data reduction equations and the generated test measurements 
the true values of map efficiency and velocity ratio were calculated, shown in 
Table A-4.1. 

The uncertainty to associate with the fourth-order regression is 
obtained by applying Eq. (33) in Appendix 3. The resulting expression is 




?-f 5i i 


SP 2 (i) 


TT Pp, (0 




i=ivaPj(o; £iisw(iv 

+2V V V d’Hmap V +2V V f Y I p 

it 1 5 T, (i) J l T, (k) J Tk Tlk hikXd T 2 (i) J U T 2 (k) J T2 ‘ T2 

,2Y vf ^ T)m>p Y d*W 1 B +2Y y 
hh Ut, © J l a t 2 (k)J Tii TJt h Jr + : Up, oJUp, ooJ b pii pit 


aT 2 (i)AaT 2 (k) 


+2Y Y f ^ r ' m * p I f 

&k&.Up 2 (i)JUp 2 (k),r~'~‘ -£!£'AdP l (i)AdP 2 (k)r ,, ‘“’ 1 ' 

|2 Y Y f £H=l¥. I fi +1 00116,341011 tennsforB TlPl» B T 2 P|, l 

i=i k=i+i < aN(i)J LaN(k)> ' 1 1 B TIN ,B T2Ni B P1N ,etc. are neglected herej (7) 
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u/co 




c 

RATIO 

0.096 

4.7119 

0.053 

95.54 

0.0266 

0.995 

0.113 

4.4896 

0.0435 

95.51 

0.0219 

0.994 

0.152 

4.0116 

0.0337 

95.28 

lilM 

0.993 

MEHM 

3.4338 

0.0282 

95.34 

0.0142 

■iLvilMI 

0.241 

3.1246 

0.0231 

95.21 

0.0116 

0.996 

0.258 

2.9824 

0.0202 

95.14 

0.0102 

0.996 

0.289 

2.7489 

0.0155 

95.23 

0.0078 

0.997 

0.310 

2.6055 

0.0134 

95.33 

0.0067 

0.999 

0.337 

2.4365 

0.0127 

95.51 

0.0064 

1.002 

0.361 

2.2957 

0.0137 

95.51 

0.0068 

1.002 

0.382 

2.1834 

0.0149 

95.37 

0.0074 

1.001 

0.394 

2.1191 

0.0156 

95.29 

0.0078 

1.001 

i 0.444 

1.8881 

0.0167 

95.28 

0.0084 

0.996 

0.455 

1.8411 

0.0167 

95.22 

0.0084 

0.994 

0.509 

1.6342 

0.0176 

95.22 

0.0089 

0.986 


WmE-Smi 

0.0194 

95.12 

0.0098 

r 0.984 

0.607 

13141 

0.0322 

95.21 

0.0163 

0.989 

0.759 

0.880 

0.036 

95.79 

0.018 

1.000 


Table A-4.2. Results of HPFT Efficiency Map Monte Carlo Simulation 

The average uncertainty ranges from 1.1% at U/C0=0.096 to 4.1% at 
U/C0=0.759. At the general mainstage operating regime of the SSME the 
uncertainty in the efficiency parameter is approximately 0.5% - 0.75%. 
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Appendix - 5 

Balance Relation Uncertainty Calculations 

Table 5.1 contains the balance relation equations and the balance 
relation uncertainty estimates. The remainder of this appendix is the output 
from a MathCAD version 6.0+ file which performs the calculations for the 
balance relation uncertainty estimates. This file is being provided to the 
COTR, Mr. John P. Butas, NASA/MSFC, Propulsion Laboratory, EP-14. As 
better estimates for the hardware characteristic are obtained this file can be 
quickly modified to obtain the updated balance relation uncertainty 
estimates. 

Balance relation uncertainties for balance relations 50-69 and 79-93 are 
essentially default values, these uncertainties are only slightly better than 
pure guesses. These value are intentionally over-estimated so that when the 
balance relation uncertainties are used in the solution methodology the 
solution will work to resolve the balance relations for which more 
information is known. When additional information is obtained to allow the 
estimation of uncertainties for these subcomponent models these balance 
relation uncertainties should be updated. 
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Table S.l, continued. Balance Relations and Balance Relation Uncertainty Estimates. 








































































































*GAMAFSF/CPFSF=WHG5*(GAMAFSF-1)»TTFSF 

RFSF*WHG5+HFROSF*WHG6+HFRVL16*WFSLV 
=HFRMFI t (WHG5+WHG6+WFSLV) 




















































FRFSF*WHG5+OFROSF*WHG6+OFRVL16*WFSLV 

=OFRMFI*(WHG5+WHG6+WFSLV) 



Table 5.1, continued. Balance Relations and Balance Relation Uncertainty Estimates. 
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These balance relation uncertainties are essentially default values, these uncertainties are only slightly better than pure guesses. 
These value are intentionally over-estimated so that when the balance relation uncertainties are used in the solution methodology 
the solution will work to resolve the balance relations for which more information is known. 
























































Balance Relation Equations, Calcuations, and Uncertainty Calculations 


Constants and information used throughout the equations. 
Joules Constant RJ =9337.92 


CLEN := 1.0 
CMASS := 1.0 
CFORCE =1.0 
CTEMP := 1.0 
GC =386.088 
FLOCON =0.53531 

The turbopump speeds are: 

SNRADFL = SNFL — 

30 

SNRADFH = SNFH — 

30 

SNRADOL = SNOL — 

30 

SNRADOH = SNOH-— 

30 


SNRADFL = 1.444* 10 3 
SNRADFH = 3.524* 10 3 
SNRADOL = 537.135 
SNRADOH = 2.885* 10 3 


SNFL = 
SNFH : 
SNOL 
SNOH : 


13785.63 
33655.01 
= 5129.258 
= 27553.1 
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Balance Relation 1 : WLPFP=WF1 

* WLPFP Is an iteration variable 

* WF1 is the flow in fuel line 1 , calculated in FLOWOO 


First, look at the calculation of WF1 and it’s uncertainty 


WFlss |aP RHQVL ° 1 
v RF1 

Consider uncertainties in the resistance and 
the density from the property routine. 


U WFl = J( 0 RFl U RFl) 2+ ( e RHOVL01 U RHOVLOl) 2 

Evaluating the partial derivatives and dividing by WF1 to obtain 
the uncertainty as a percentage. 


WLPFP := 149.0299 
WF1 := 149.0299 
RF1 =3.81 10" 6 
U rfi : = 0.10 RF1 

RHOVLOl :=2.56- 10" 3 
U RHOVLOl ' = 0. .001 RHOVLOl 


UH TO1 := 

ri <J\ 2 RF1 / \2 RHOVLOl / 

U% wfi = 5*% 

U WF1 -WFlU'/owp! 

U \ypj =7.452 


Now rewrite the balance relation as 
U BAL1 =U WFI 


BALI := WFI - WLPFP 
BALI =0 


U BALI = 7452 


Now write the uncertainty as a percentage of the right-hand side of the balance relation: 


U% 


u 


BALI 


BALI 


WFI 


U% BAL1 = 5 * % 
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Balance relation 2: WLPFP*HTHTNK+QDOTLPFP=WF1*HTVL01 

* WLPFP is an iteration variable 

* HTHTNK is the H2 tank enthalpy and is an iteration variable 

* QDOTLPFP is the power required by the pump, calc'd in SPUMP 

* WF1 is the flow in fuel line 1 

* HTVL01 is the volume 1 enthalpy and is an iteration variable 

Now rewrite the balance relation and determine the uncertainty 
BAL2=WF1HTVL01 - WLPFP- HTHTNK- QDOTLPFP 

The LPFP torque is calculated in SPUMP and uses a map for TAU vs PHI 

TORQLPFPs- TAULPFP- TORQMLPFP- RHO VLO 1 SNRADFL 2 

TORQLPFP =- 10521.79 


HTVL01 :=- 96.92439 
HTHTNK =-107.8393 
WLPFP := 149.0299 
TAULPFP := 1.971925 

SNFL = 1.379* 10 4 
SNRADFL = 1.444- 10 3 


TRQPFL := TORQLPFP 


Considering an uncertainty associated with the torque coefficient vs flow coefficient map, and an 
uncertainty in density due to the property routine. 


U% 


TORQLPFP - 


u TAULPFP V 2 / u RHO VLO 1 
TAULPFP / + \ RHO VLO 1 


2 


U tauLPFP :: 0.05- TAULPFP 


U% TORQLPFP =5001 * % 

U TORQLPFP | U% TORQLPFP T0R Q LPFP | 
U TORQLPFP = 526195 


QDOTLPFP :=- TORQLPFP- 


SNRADFL 

RJ 


QDOTLPFP = 1.627- 10 3 

a -SNRADFL 

0 TORQLPFP = — 


U 


QDOTLPFP ( 0 TORQLPFP' U TORQLPFp) 2 
U QDOTLPFP = 81 349 


.... _ U QDOTLPFP 

u/ « QDOTLPFP - QDOTLPFP 

U% QDOTLPFP = 5001 * 0/o 
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Returning to the balance relation uncertainty, 

and noting that the uncertainty sources are U WF1 AND U QDOTLPFP 


BAL2 := WF1-HTVL01 - WLPFP HTHTNK - QDOTLPFP 
6^, -HTVL01 BAL2 =-0.004 

0 QDOTLPFP : = - 1 


U BAL2 := J( e WFl' U WFl) 2+ ( 0 QDOTLPFP U QDOTLPFp) 2 
U BAL2= 726834 

Now writing the uncertainty as a percentage of the RHS of the balance relation: 
TT0/ U BAL2 

U70 DAT O • 

| WF1HTVL01 1 
U % BAL2 = 5032*% 
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Balance Relation 3: WF1=WHPFP 

* WHPFP is an iteration variable 

* WF1 is the flow in fuel line 1 , calculated in FLOWOO 


WHPFP := 149.0299 


Recalling the WF1 uncertainty 

U % WF1 =5 * % 

Uwn =7.452 


Now rewrite the balance relation as BAL3 := WHPFP - WF1 BAL3 = 0 


U BAL3 :=U WF1 
U BAL3= 7 - 452 


Now write the uncertainty as a percentage of the RHS of the balance relation" 


U% BAL3 := 
U% BAL3 


U BAL3 
WHPFP 

= 5 •% 
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Balance relation 4: WF1 *HTVL01 +QDOTVL02=WHPFP*HTVL02 

* WHPFP is an iteration variable 

* HTVL01 is the volume 2 enthalpy and is an iteration variable 

* QDOTHPFP is the power required by the pump, calc’d in SPUMP 

* WF1 is the flow in fuel line 1 

* HTVL01 is the volume 1 enthalpy and is an iteration variable 
Now rewrite the balance relation and determine the uncertainty 

BAL4 := WHPFP HTVL02 - WF1 HTVL01 - QDOTVL02 
BAL4 =0.001 
0 Wpi = - HTVL01 


HTVL02 :=- 96.92438 
QDOTVL02 =0 


U BAL4 : -J( e WFr U WFl ) 2 
U BAL4= 722268 


Now write the uncertainty as a percentage of the RHS of the balance relation" 


U% BAL4 := - 

U% BAL4 


U BAL4 

| WHPFPHTVL02 
= 5 •% 
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Balance Relation 5: WHPFP=WMFV+WFTC 

* WHPFP is an iteration variable 

* WMFV is the main fuel valve flow, calculated in FLOWOO 
Now rewrite the balance relation as: WMFV+WFTC-WHPFP 

Looking at the flow and uncertainty calculations for the MFV 


WMFV= 


AP- 




RHOVLQ3 

RMFV 


Consider uncertainties in the resistance and 
the density from the property routine. 


4 


U WF1 (®RMFV' U RMFV 


0 RHOVL03 U RHOVL03 


Evaluating the partial derivatives and dividing by WF1 to obtain 
the uncertainty as a percentage. 


1 U RHOVLQ3 \ 2 

2 RHOVL03 / 


u WMFV - 7 255 

Looking at the flow and uncertainty calculations for the FTC 

WFTC=V APRHO VL03CFFTC 2 FLOCON 

Consider uncertainties in the resistance and 
the density from the property routine. 

U WFT crjfi CFFT C U CFFTc) 2+ ( e RHOVL03 U RHOVL03) 2 

Evaluating the partial derivatives and dividing by WF1 to obtain 
the uncertainty as a percentage. 

T to/ ( U C FFTC\ 2 /l u RHOVL03\ 2 

0 WFTC CFFTC / \2 RHOVL03 / 

u %wftc = 10 - 44 * 0/o 
U WFTC : =WFTC-U% W FTC 

U WFTC =0 41 


U% WMFV : - 


1 1 U RMFV\ 


\ 2 RMFV 
U ° /o WMFV = 5 *° /o 


U WMFV - WMFVU% v/mpv 


WMFV := 145.1001 

RMFV := 1.53- 10' 5 
U RMFV : = 01 ° ] RMFV 

RHOVL03 := 2.89- 10" 3 
U RHOVL03 r 0.001 RHOVL03 


WFTC =3.929894 

CFFTC : = 2.025 
U - 0. 1 0* CFFTC 

RHOVL03 := 2.89- 10' 3 
U RHOVL03 -0 06 RHOVL03 
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Rewriting the balance relation and determining the uncertainty 

BAL5 := WMFV + WFTC - WHPFP 

BAL5 =9.4-1 (f 5 

U BAL5 : = WMFV 2 + u WFTC 2 
U BAL5= 7267 

Now write the uncertainty as a percentage of the RHS of the balance relation" 

TT0/ U BAL5 

U /O 13 AT C ■ ; 

| WMFV + WFTC | 

U% BAL5 =4 - 876 * % 
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Balance relation 6: WHPFP*HTVL02+QDOTHPFP=(WMFV+WFTC)*HTVL03 

* WHPFP is an iteration variable 

* H7VL02 is the volume 2 enthalpy and is an iteration variable 

* QDOTHPFP is the power required by the HPFTP pump, calc'd in SPUMP 

* WMFV is the flow through the MFV 

* WFTC is the flow through the fuel turbine cooling line 

* HTVL03 is the volume 3 enthalpy and is an iteration variable 


Now rewrite the balance relation and determine the uncertainty 
BAL6=(WMFV+WFTC)*HTVL03-WHPFP*HTVL02-QDOTHPFP 


HTVL03 := 176.0563 


The required pump power is determined from the SPUMP module and the pump map uncertainties. 


QDOTHPFP=TORQHPFP 


SNRADFH 

RJ 


TORQHPFP is calculated using TAUHPFP from the torque 
parameter vs flow coefficient map. 

TORQHPFPs- TAUHPFP- TORQMHPFP-RHOVL03 • SNRADFH 2 
TORQHPFP =-107789.8 

TRQPFH := | TORQHPFP | 


TAUHPFP := 3.003558 
SNFH = 3.366- 10 4 
SNRADFH = 3.524* 10 3 


Considering an uncertainty associated with the torque coefficient vs flow coefficient map, and an 
uncertainty in density due to the property routine. 

U TAUHPFP = 0.05- TAUHPFP 

U TAUHPFP | / U RHOVLQ3 \ 2 
TAUHPFP / + \ RHOVL03 / 


U% TORQHPFP :: 


U% TORQHPFP “ 7 81 * 0/o 

U TORQHPFP := | U% TORQHPFP TORQHPFP 
U TORQHPFP = 8 -4 1 9- 1 0 3 

U TRQPFH = U TORQHPFP 

u trqpfh =8419 ‘ 1C)3 


QDOTHPFP :=- TORQHPFP 
QDOTHPFP = 4.068* 10 4 


SNRADFH 

RJ 


0 TORQHPFP 


SNRADFH 

RJ 


SNFH:= 33655.01 

SNRADFH : = SNFH-— 

30 

SNRADFH = 3.524344- 10 3 
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u 


QDOTHPFP - J( e TORQHPFP' U TORQHPFp) 2 
U QDOTHPFP =3 177-10 3 


U% QDOTHPFP := - 


U 


QDOTHPFP 
QDOTHPFP 

U% QDOTHPFP = 781 * % 


Returning to the balance relation and its uncertainty 


BAL6 -(WMFV+ WFTC) HTVL03 - WHPFP HTVL02 - QDOTHPFP 
BAL6 =-0.034 

The partial derivatives, 0’s, are 
© WMFV = HTVL03 
0 • = HTVL03 

6 QDOTHPFP := ‘ 1 


so that the uncertainty expression for balance relation 6 becomes 

U BAL6 : = J ( 0 WMFV U WMFv) 2+ ( e WFTC -U WFTc) 2 -*- ( 0 QDOTHPFP U QDOTHPFp) 2 
^ BAL6 = 3-425* 10 3 

Now write the uncertainty as a percentage of the RHS of the balance relation" 

u % BAL6 := UBAL6 

DALO |(WMFV+WFTC)HTVL03| 

U% BAL6 = 13055 
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Balance Relation 7: WMFV=WF2+WF6+WF8+WFIG+WFFI+WFOI 

* WMFV is the main fuel valve flow, calculated in FLOWOO 

* WF2 is the flow through fuel line 2, calculated in FLOWOO 

* WF6 is the flow through fuel line 6, calculated in FLOWOO 

* WF8 is the flow through fuel line 6, calculated in FLOWOO 

* WFIG is the flow through fuel line, calculated in PIPE01 

* WFFI is the flowthrough fuel line, calculated in PIPE01 

* WFOI is the flow through fuel line, calculated in PIPE01 


Now rewrite the balance relation as: BAL7=WMFV-WF2-WF6-WF8-WFIG-WFFI-WFOI 


First, look at the calculation of WF2 and it’s uncertainty 


WF2= 


AP- 


RHQVL04 

RF2 


WF2 =56.02199 


Consider uncertainties in the resistance and 
the density from the property routine. 


U WF2 = J( 0 RF2 U RF2) + ( e RHOVL04 U RHOVL04) 2 


RF2 =4.59- 10' 5 
RF6 = 3.0010' 4 
RF8 = 4.43- 10" 4 
CFFIG : = 0.42 
CFFFI =0.745 
CFFOI: = 0.799 


Evaluating the partial derivatives and dividing by WF2 to obtain 
the uncertainty as a percentage. 


U% 


WF2 - 


1 U RF2\ 2 /l U RHOVL04\ 2 


2 RF2 


+ — 


U j^p2 • = 1 RF2 
RHOVL04 =2.87- 10' 3 

U RHOVL04 :=0 001RH < 


2 RHOVL04 


U % WF2 = 5 '% 


UwF2 =WF2 U%WF2 
U ^p2 = ^-80 ^ 


First, look at the calculation of WF6 and it's uncertainty 

WF6= L* ™™* 

*] RF6 

WF6 := 54.64435 


Consider uncertainties in the resistance and 
the density from the property routine. 


U 


WF6 




) RF6 U RF6)" 


( 0 RHOVLO4 U RHOVL04) 2 
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Evaluating the partial derivatives and dividing by WF6 to obtain 
the uncertainty as a percentage. 


U jyrg -0.1-RF6 


U% 


WF6 


If I U RF6 \ 2 | /l U RHOVLQ4 \ 2 
'J \ 2 RF6 / \2 RHOVL04 / 


u %wF6 =005 
U W p 6 := WF6 U% W p 6 
UwF6= 2-732 

First, look at the calculation of WF8 and it’s uncertainty 

WF8= IaP™?™ 

V RF8 

Consider uncertainties in the resistance and 
the density from the property routine. 


Evaluating the partial derivatives and dividing by WF8 to obtain 
the uncertainty as a percentage. 


U % WF8 "J 
U % WF8 = 


f 1 U RF8 ^ 2 

\ 2 RF8 / 
0.05 


+ 


1 U RHQVL04 

2 RHOVL04 


2 


U WF8 - WF8-U% y/pg 
U WF8 = 1595 


Looking at the flow and uncertainty calculations for the FIG 


WFIG=V AP- RHOVL04 CFFIG 2 -FLOCON 

Consider uncertainties in the resistance and 
the density from the property routine. 

U WFIG = J( e CFFIG U CFFIg) 2+ ( e RHOVL04' U RHOVL04) 2 


Evaluating the partial derivatives and dividing by WFIG to obtain 
the uncertainty as a percentage. 


t to/ _ , U CFFIg\ /l u RHOVL04 i 
wnu ii CFFIG \ 2 RHOVL04 


U% y/pjQ -10*% 


U WFIG - WF 10 u% WFIG 
U WFIG = ® 


WF8 := 31.90515 
U :=0.1RF8 


WFIG =0.85539 
U CFFIG " 01CFFIG 
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Looking at the flow and uncertainty calculations for the FFI 


WFFI=-J AP-RHO VL04 ■ CFFFI 2 FLOCON WFFI =0.85539 

Consider uncertainties in the resistance and ^ CFFFI _0 1 CFFFI 

the density from the property routine. 

U WFF rj ( 6 CFFFI U CFFFl) 2 + ( e RHOVL04 U RHOVL04) 2 


Evaluating the partial derivatives and dividing by WF1 to obtain 
the uncertainty as a percentage. 

tto/ / U CFFFl\ 2 /l U RHOVL04\ 2 

0 WFFI CFFFI / + \2 RHOVL04 / 

U% WFFI = 1 0 * % 

U wpH : = 

u wfpi =0.086 


Looking at the flow and uncertainty calculations for the FOI 

WFOI=a/aP- RHO VL04- CFFOI 2 FLOCON 


WFOI := 0.85539 
U CFFOI = 0-1' CFFOI 


Consider uncertainties in the resistance and 
the density from the property routine. 


wr KJL 


Evaluating the partial derivatives and dividing by WF1 to obtain 
the uncertainty as a percentage. 


U CFFOl 'j 2 (]_ U RHQVL04 
, CFFOI j + \2 RHOVL04 

u % WFOI =01 
U WFOI := WFOI U% y/poi 


U % WFOI := 


2 


U WFOI - 0086 

The uncertainty expression for BAL7 reduces to 


U 


BAL7 


:= U 


WF2 


•+u 


WF6 


-t-U 


WF8 


•+u 


WFIG 


+ U 


WFFI 


•+u 


WFOI 


+ U 


WMFV 


U BAL7= 8 ' 398 

Now write the uncertainty as a percentage of the RHS of the balance relation" 

U % BAL7 

|WF2 + WF6 + WF8 + WFIG + WFFI + WFOl| 

U°/o B AL7 = 5.786-% 
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Balance relation 8: 

WMFV* (HTVL03-UTVL04)= (WF2+ WF6+ WF8+ WFIG+WFFI+ WFOI)*(HTVL04-UTVL04) 

* WMFV, WF2,WF6,WF8 are calculated in FLOWOO and have uncertainties 

* WFIG, WFFI, WFOI are calculated in PIPE01 and have uncertainties 

* HTVL03, HTVL04 are iteration variables 

* UTVL04 is calculated in SSME3.CFG 

HTVL04 and PTVL04 are iteration variables: HTVL04 = 176.0563 

UTVL04 :=HTVL04 - — . PTVL04 PTVL04 =5709.895 

RJ RHOVL04 

UTVL04 =-37.001 


Thus UTVL04 has uncertainty due to the 
uncertainty from the property routine. 



PTVL04 

RHOVL04 2 


"this value for UTVL04 is slightly 
different that the value in the output 
file used of -37.34206 


U UTVL04 : = J (© RHOVL04 U RHOVL04) 2 

U UTVL04 =0 ' 213 

Balance relation 8 is rewritten as: 

BAL8 := (WF2 + WF6 + WF8 + WFIG + WFFI + WFOI)- (HTVL04 - UTVL04) - WMFV (HTVL03 - UTVL04) 
BAL8= 8.002 

0 WF2 =HTVL04 - UTVL04 ©WFIG = HTVL04 - UTVL04 

0 WF6 = HTVL04 - UTVL04 ©WFFI = HT VL04 - UT VL04 

0 WF8 := HTVL 04 - UTVL04 0 WFQl : = HTVL04 - UTVL04 

©UTVL04 : =-WF2- WF6- WF8- WFIG- WFFI- WFOI+ WMFV 
0 WMFV : = - HTVL03 + UT VL04 

U BAL8 = 

U BAL8 = 1789 * 103 

Now write the uncertainty as a percentage of the RHS of the balance relation” 

TT 0/„ - U BAL8 

|( WF2 + WF6 + WF8 + WFIG + WFFI+ WFOI) (HTVL04 - UTVL04) | 

U% BAL8 = 5786 * % 
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Balance Relation 9: WF2=WF3 
* WF2 and WF3 are calculated in FLOWOO 

First, look at the calculation of WF3 and if s uncertainty 

^ RHQVL05 
RF3 

Consider uncertainties in the resistance and 
the density from the property routine. 


U WF3 = J ( e RF3 U RF3) 2+ ( 0 RHOVLO5 U RHOVLO5) 2 



WF3 := 56.0222 
RF3 =4.89- 10' 5 
U RF3 : = 0.1 RF3 


Evaluating the partial derivatives and dividing by WF2 to obtain 
the uncertainty as a percentage. 

TT0/ / 1 U RF3\ 2 /l U RHOVL05\ 2 

0WF3 'Ji2 RF3 / + \2 RHOVL05 / 

U%wp3 = 5*% 


RHOVL05 :=2.89- 10' 3 
U RHOVL05 = 0.001 RHOVL05 


u WF3 -WF3 U% ^3 
U ■yyp3 =2.801 

Rewriting the balance relation and determining the uncertainty 


BAL9 := WF3 - WF2 


BAL9 =2.1*10 4 


I 2 2 

U BAL9 := V U WF2 - | ' U WF3 
U BAL9 = 3.962 


Now write the uncertainty as a percentage of the RHS of the balance relation 


U% BAL9 :: 
U% BAL9 


_ U BAL9 
’ |WF3| 

= 7.071*% 
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Balance Relation 10: WF2*HTVL04+QDOTVL05=WF3*HTVL05 

* WF2, WF3 are calculated in FLOWOO and uncertainties as above 

* HTVL04, HTVL05 are iteration variables 

* QDOTVL05 is defined as zero for this relation 

Rewriting the balance relation and determining the uncertainty: 

BAL10 := WF3-HTVL05 - WF2HTVL04- QDOTVL05 
BAL10 =0.043 


®WF3 := HTVL05 
9 WF2 ■ ~ ■ HTVL04 

U BAL10 := J( 0 WF3' U WF3) 2+ ( e WF2' U WF2) 2 
U BAL10= 697 ' 458 


Now write the uncertainty as a percentage of the RHS of the balance relation' 


U% BAL10 : 
U% BAL10 


U BAL10 
| WF3HTVL05 

= 7.071 •% 


HTVL05 := 176.0564 
QDOTVL05 =0 


A5-22 



Balance Relation 1 1 : WF3=WF4+WNLK1 

* WF3, WF4 are calculated in FLOWOO 

* WNLK1 is the nozzle leakage flow and is currently zeroed out WNLK1 : = 0 

Rewrite the balance relation as: BAL1 1 =WF4+WNLK1-WF3 


First, look at the calculation of WF3 and it’s uncertainty 

^ RHQVL06 
RF4 

WF4 =56.02217 

Consider uncertainties in the resistance and 
the density from the property routine. 

U WF3 = J ( e RF4' U RF4) 2+ ( 0 RHOVLO6 U RHOVLO6) 2 


WF4= 


RF4 =1.05 10 -4 
u PJP 4 := 0.1 RF4 


RHOVL06 = 2.84- 10' 3 
U RHOVL06 = 0.001RHOVL06 


Evaluating the partial derivatives and dividing by WF2 to obtain 
the uncertainty as a percentage. 


U% 


WF4 •= 


1 U RF4| / 1 U RHOVL06\ 2 


2 RF4 


+ 


2 RHOVL06 


U% y/F4 - 5 •% 

U WF4 :=WF4U% WF4 
U =2.801 

Rewriting the balance relation and determining the uncertainty 

BAL11 := WF4 +• WNLK1 - WF3 

BAL11 =-3-10~ 5 


U 


BALI 1 


= U 


WF4 


•+u 


WF3 


U BALI 1 “ 3962 


Now write the uncertainty as a percentage of the RHS of the balance relation" 

v% BALll ■, UEAL11 

DALU | WF4+ WNLKl| 

U% BAL11 = 7071 ’° /o 
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Balance Relation 12: WF3*(HTVL05-UTVL06)=(WF4+WNLK1)*(HTVL06-UTVL06) 

* WF3, WF4 are calculated in FLOWOO, as shown above 

* WNLK1 is defined as zero at this point 

* HTVL05,HTVL06 are iteration variables 

* UTVL06 is calculated in the configuration file 


HTVL06 and PTVL06 are iteration variables: 


HTVL06 := 176.0591 


UTVL06 :=HTVL06 - — — I— 06 
RJ RHOVL06 


Thus UTN/L06 has uncertainty due to the 
uncertainty from the property routine. 


e RHOVL06 


1 PTVL06 
RJ RHOVL06 2 


PTVL06 := 5605.869 


UTVL06 =-35.326 

**this value for UTVL06 is slightly 
different that the value in the output 
file used of -35.09946 


U 


UTVL06 


: = /]( 0 RHOVLO6 U RHOVLO6) 2 


U UTVL06 * 0 21 1 


Now rewriting the balance relation and determining the uncertainty 

BAL12 := ( WF4 + WNLK1 ) (HTVL06 - UTVL06) - WF3 (HTVL05 - UTVL06) 
BALI 2 =0.145 
6 W4 : = HTVL06 - UT VL06 

e W3 :=-HTVL05 + UTVL06 
e UTVL06 --WF4- WNLK1 + WF3 

U BAL12 := J( e WF4' U WF4) 2 + (® WF3 U WF3) 2+ ( e UTVL06 U UTVL06) 2 
U BAL12 = 837411 

Now write the uncertainty as a percentage of the RHS of the balance relation" 

T to/. ._ U BAL12 

U/OpAT 10 

j ( WF4 - WNLK1 ) (HTVL06 - UTVL06)| 

U% BAL12 =7071 *° /o 


A5-24 



Balance Relation 13: WF4=WF5 
* WF4.WF5 are calculated in FLOWOO 

First, look at the calculation of WF5 and it's uncertainty 

._ RHOVL07 

or- 

RF4 

Consider uncertainties in the resistance and 
the density from the property routine. 

U WF5 = J( 0 RF4' U RF4) 2+ ( 0 RHOVLO7 U RHOVL07) 2 

Evaluating the partial derivatives and dividing by WF5 to obtain 
the uncertainty as a percentage. 

TT0/ / 1 U RF5\ 2 [l U RHOVL07\ 2 

0WF5 J\ 2 RF5 / \2 RHOVL07 / 

U% Wp5 = 5.831 •% 

U \^/p5 = WF5U% \^p5 
U WF5 “ 3.267 

Rewriting the balance relation and determining the uncertainty 
BALD := WF5 - WF4 



WF5 := 56.02202 
RF5 := 1.13- 10 -4 
UrF5 := 01RF5 

RHOVL07 =1.46- 10’ 3 
U RHOVL07 : = 0.06 RHOVL07 


BALD =-1.5-10 4 


U BALD : = V U WF5 2+U WF4 2 
U BAL13 =4 - 303 


Now write the uncertainty as a percentage of the RHS of the balance relation” 


U% BAL13 : 


U BALD 
|WF5| 


U ° /o BAL13 “ 7 681 * 0/o 
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Balance Relation 14: WF4*HTVL06+QDOTVL07+QAMBVL07=WF5*HTVL07 

* WF4.WF5 are shown above 

* HTVL06,HTV1.07 are the volume 6 and 7 enthalpies and are iteration variables 

* QHOTVL07 is the heat transfer from the nozzle to the coolant flow 

* QAMBVL07 is the heat transfer from the ambient to the coolant flow 


HTVL06 := 176.0591 
HTVL07 =918.2539 


QDOTVL07 =41579.38 
QAMBVL07 : = 0 


Heat transfers, Q's, are calculated in VOLQ03 by calling the QH2 subroutines, 
which calculates the heat flux, QDOT, through the areas, AHOTVL07 AND AAMBVL07. 


AAMBVL07 : = 28670 


AHOTVL07 := 14490 


QDOT is calculated in QH2 with the heat transfer correlation equation: 


QDOT- 


0.0303 W 
D 18 W 2 



QDOT hot 


QDOTVLQ7 

AHOTVL07 


This equation is highly empirical and trying to base an uncertainty by propagating the estimated 
uncertainties of the individual terms would not be advisable. Some of the areas which contribute to 
the heat flux uncertainty are: 

* There is no way to assess an uncertainty for the constants and exponents. 

* The specific heat, viscosity, and thermal conductivity are difficult to obtain for super-critical 
hydrogen, so they will have considerable uncertainty. 

* Values used for the temperature of the fluid, TF, and metal temperature, TM, are difficult to 
model. 

Thus, the best way to estimate an uncertainty for the heat transfer flux is to examine the 
information upon which the flux equation was developed. There are no comment lines in the 
subroutine discussing the development of the equation and no other information is available at this 
time. Hense, an uncertainty of 20% will be assigned to the heat flux, until information is obtained 
upon which a better estimate can be made. 


U% QDOT - 0 20 u qdoTHOT " U% QDOT Q 001 HOT 

U QDOTHOT = 0574 


AHOTVL07 is the heat transfer area representing the hot wall of the nozzle. The nature of the 
design and manufacture of the nozzle and the multiple tubes which make up the nozzle coolant 
passages make an acurate determination of the heat transfer area difficult. Thus an uncertainty of 
20% is being assigned. 

QDOTVL07 : = QDOT Hor AHOT VL07 U AHOT = 0 - 20 ‘ AHOT VL07 


0 AREA -Q D0T HOT 
0pLUX : = AHOT VL07 
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U QDOTVL07 : ~J( e FLUX' U QDOTHOt) 2+ ( 0 AREA U AHOt) 2 
U QDOTVL07 = 1-176* 10 4 

Rewriting trie balance relation and propagating the uncertainties: 


BALM := WF5 HTVL07 - WF4-HTVL06 - QDOTVL07 - QAMBVL07 
BALM =-0.154 

The partial derivatives needed are: 

0 WF5 := HTVL07 
® QDOTVL07 -- 1 


0 


WF4 


:=-HTVL06 


U BAL14 : = -j( e WF5 U WF5) +( e WF4 U WF4) + (® QDOTVL07' U QDOTVL07 
U BAL14 = 1 - 215 * 104 

Now write the uncertainty as a percentage of the RHS of the balance relation" 

U BALM 


U% 


BALM 


| WF5HTVL07| 
U% BALM = 23 613 - % 
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Balance Relation 15: WF5+WCCV=WF7 

* WF5.WF7 are fuel line flows and calculated in FLOWOO 

* WCCV is the flow through the coolant control valve and is calculated in FLOWOO 


Now rewrite the balance relation as: BAL15=WF7-WF5-WCCV 

First, look at the calculation of WF7 and ifs uncertainty 

RHOVLQ8 
RF7 


WF7= AP 


WF7 =110.667 
RHOVL08 := 1.85- 10' 3 
RF7 =4.02- 10’ 6 


Consider uncertainties in the resistance and the density from the property routine. Estimate 10% 
uncertainty in resistance and 2% uncertainty in density from the property routine. 


U WFrJ ( e RF7 U RF7) 2+ ‘ ( 0 RHOVLO8 U RHOVLO8) 2 


Evaluating the partial derivatives and dividing by WF8 to obtain 
the uncertainty as a percentage. 


U% 


WF7 := 

^1 


1 U RF7 ^ 

2 RF7 


2 


1 U RHOVLQ8 

2 RHOVL08 


U% WF7 — 5 .831 •% 
u WF7 = WF7 -U°/o yfpy 


U y/jrj —6.453 


U Rp7 =0.1 RF7 
U RHOVL08 = 0 06RHOVL08 


And similarly for the Coolant Control Valve. 


wccv= 


fcp RHQVL09 
RCCV 


RHOVL09 := 2.8- 10' 3 
RCCV =6.49- 10' 5 

WCCV =54.64416 


Consider uncertainties in the resistance and the density from the property routine. Estimate 10% 
uncertainty in resistance and 0.2% uncertainty in density from the property routine. 


U WC CV=J ( 0 RCCV U RCCv) 2+ - (® RHOVL09 U RHOVL09) 2 

Evaluating the partial derivatives and dividing by WCCV to obtain U RCCV -O.I RCCV 

the uncertainty as a percentage. U rhqvL09 := 0 RHOVL09 

TT0/ I 1 U RCCV \ 2 /l U RHOVL09\ 2 

°WCCV 2 RCCV j + ^ RHOVLQ9 j 

U% WCCV = 5 * % 

U WCCV : = wcc V U% WCCV 
U WCCV = 2 732 
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Now, rewriting the balance relation as 

BALI 5 := WF7 - WF5 - WCCV 

BAL15 = 8.2*10 4 


©WF5 :=1 
6 WCCV = - ] 

©WF7 :=1 

U BAL15 : = /( e WF7‘ U WF7) 2+ ( e WF5' U WF5) 2+ ( 0 WCCV U WCCv) 2 
U BAL15= 7732 

Now write the uncertainty as a percentage of the RHS of the balance relation 


U% BAL15 :: 
U% BAL15 


U BAL15 
|WF7| 

= 6.986*% 
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Balance Relation 16: WF5*(HTVL07-UTVL08)+WCCV*(HTVL09-UTVL08)=WF7‘(HWL08-UTVL08) 

* WF5,WCCV,WFV are discussed above 

* HTVL07,HTVL08,HTVL09 are iteration variables 

* UTVL08 is calculated from available information. . „ 

HTVL08 =551.7748 

HTVL08 and PTVL08 are iteration variables: HTVL09 = 176.0563 


UTVL08 :=HTVL08 - — PTVL08 
RJ RHOVL08 


Thus UTVL08 has uncertainty due to the 
uncertainty from the property routine. 


PTVL08 := 5328.137 


UTVL08 =243.347 

"this value for UTVL08 is slightly 
different that the value in the output 
file of 243.9922 


® RHOVL08 := 


1 PTVL08 
RJ RHOVL08 2 


U UTVL08 : = V ( 0 RHOVLO8 U RHOVL08) 


U UTVL08 “ 18506 


Now rewriting the balance relation and determining the uncertainty 

BALI 6 := WF7 (HTVL08- UTVL08) - WF5 (HTVL07 - UTVL08) - WCCV-(HTVL09 - UTVL08) 
BALI 6 =0.175 
The partial derivatives are: 

9 wp; : = HTVL08 - UT VL08 
0 ^5 : = - HTVL07 + UTVL08 
9 WCCV =- HTVL09 + UTVL08 
e UTVL08 : =- W7 +• WF5 +- WCCV 

U BAL16 := J( e WF7 U WF7) 2+ ( e WF5' U WF5) 2+ ( e WCCV U WCCv) 2 + ( e UTVL08' U UTVL08) 2 
U BAL16 =2976 * 103 


Now write the uncertainty as a percentage of the RHS of the balance relation" 

U BAL16 


U% 


BAL16 - 


|WF7 (HTVL08- UTVL08)| 
U#/o BAL16 =8.718*% 
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Balance Relation 17. WF6=WCCV 
* WF6.WCCV are calculated in FLOWOO 

First, look at the calculation of WF6 and it’s uncertainty 


^ RHOVLQ4 
RF6 

Consider uncertainties in the resistance and 
the density from the property routine. 

U WF6 = J ( e RF6' U RF6) + ( 0 RHO VL04' U RHOVLOf) 2 

Evaluating the partial derivatives and dividing by WF6 to obtain 
the uncertainty as a percentage. 



U% wF6 ~ 3 ‘ 0// ° 

UWF6 =WF6 ' U %WF6 

UwF6= 2 - 732 

Rewriting the balance relation and determining the uncertainty 


WF6= 


WF6 =54.64435 

RF6 := 1.13- 10" 4 
Urf6 : =0.1RF5 

RHOVL04 := 2.84- 10" 3 
U RHOVL04 : = 0.001 -RHOVL04 


BALI 7 := WCCV- WF6 
BALI 7 = —1 .9* 10 _ 4 


1 2 2 
U BAL17 : = V U WCCV +U WF6 

U BAL17= 3 - 864 

Now write the uncertainty as a percentage of the RHS of the balance relation" 
,_ U BAL17 

U/ ° BAL17 ’"IwccvT 

U% BAL17 =7071 * % 
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Balance Relation 18: WF6*HTVL04+QDOTVL09=WCCV*HTVL09 

* WF6, WCCV are discussed above 

* HTVL04, HTVL09 are iteration variables 

* QDOTVL09 is zeroed in this model version 


Rewriting the balance relation and determining the uncertainty 

BALI 8 := WCCV HTVL09 - WF6 HTVL04- QDOTVL09 
BALI 8 =-0.033 

The partial derivatives are: 

0 WCCV = HTVL09 

0WF 6 :=-HTVL04 

U BALI 8 := J( e WCCV U WCCv) 2+ (®WF6 U WF6) 2 


U BALI 8 = 680.304 


Now write the uncertainty as a percentage of the RHS of the balance relation" 


U% BAL18 

U% BAL18 


u BALI 8 
| WCCVHTVL09 

= 7.071*% 


QDOTVL09 =0 
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WF9 := 3 1 .90514 
RF9 :=10.10' 4 
U RF9 : = 0.1-RF9 

RHOVLIO =2.83- 10" 3 
U RHOVLIO = 0.00 1 -RHOVLIO 

U WF9 = J ( 6 RF9' U RF9 ) 2 + ( 0 RHOVL 1 O’ U RHO VL1 o) 2 

Evaluating the partial derivatives and dividing by WF9 to obtain 
the uncertainty as a percentage. 

U%wp9 

u %WF9 =5 *° /o 
U : = WF9- U% 

u WF 9 = l595 

Rewriting the balance relation and determining the uncertainty 

BALI 9 := WF9 - WF8 

BAL19 =-l*10 _5 



Balance Relation 19: WF8=WF9 
* WF8.WF9 are calculated in FLOWOO 

First, look at the calculation of WF9 and its uncertainty 
WF9= 

Consider uncertainties in the resistance and 
the density from the property routine. 




U BALI 9” 2256 

Now write the uncertainty as a percentage of the RHS of the balance relation" 


U% BAL19 := 


u BAL19 
|WF9| 


U ° /o BAL19= 7071 - % 
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Balance Relation 20: WF8*HTVL04+QDOTVL10=WF9*HTVL10 

* WF8, WF9 are discussed above 

* HTVL04, HTVL1 0 are iteration variables 

* QDOTVLIO is zeroed in this model version 

Rewriting the balance relation and determining the uncertainty 

BAL20 := WF9 HTVL10 - WF8HTVL04 - QDOTVLIO 

BAL20 =-0.005 

The partial derivatives are: 

0 WF9 - HTVL10 

®WF8 ■ = ' HTVL04 

U BAL20 := J( 0 WF9‘ U WF9) 2+ ( e WF8‘ U WF8) 2 
U BAL20= 397209 

Now write the uncertainty as a percentage of the RHS of the balance relation" 

TT0/ U BAL20 

| WF9HTVL10| 

U% BAL20 =7071 * % 


QDOTVLIO =0 
HTVL10 := 176.0562 
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Balance Relation 21 : WF9=WF10+WMLK1 

* WF9.WF10 are calculated in FLOWOO 

* WMLK1 is a leakage flow, zeroed out in this model 


First, look at the calculation of WF9 and ifs uncertainty 

WF10= AP — 0VL11 WF10 =31.90514 

4 RF10 

Consider uncertainties in the resistance and 
the density from the property routine. 


RF10 := 21.7596- 10' 4 
U |^p]Q -0.1 RF10 

RHOVL11 :=2.75- 10' 3 
U RHOVL1 1 = 0.001 RHOVL1 1 

WMLK1 := 0 


U 


WF10 




RF10' U RF10 


) + ( e RHOVL1 1 U RHOVL1 1 ) 2 


Evaluating the partial derivatives and dividing by WF10 to obtain 
the uncertainty as a percentage. 


U% WF10 :_ 


/ 1 U RFlp \ 2 1 1 u RHOVL1 1 \ 
4\ 2 RF10 / + \2 RHOVL11 j 


2 


U% WF10 =5 * % 

U WF10 :=W F1° U% WF10 
U WF10 = 1 - 595 

Rewriting the balance relation and determining the uncertainty 


BAL21 := WF10 + WMLK1 - WF9 


BAL21 =0 


/ 2 2 
U BAL21 : = V U WF10 +U WF9 

U BAL21 =2 - 256 

Now write the uncertainty as a percentage of the RHS of the balance relation" 

U BAL21 
WF10+ WMLK1 1 

7.071 •% 


U% BAL21 : --j 
U% BAL21 = 
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Balance Relation 22: WF9*(HTVL1 0-UTVL1 1 )=(WF1 0+WMLK1 )*(HTVL1 1 -UTVL1 1 ) 

* WF9, WF10 are calculated in FLOWOO, as shown above 

* WMLK1 is defined as zero at this point 

* HTVL10,HTVL11 are iteration variables 

* U7VL1 1 is calculated in the configuration file 


HTVL1 1 and PTVL1 1 are iteration variables: 


HTVL11 := 176.0558 
PTVL11 := 5193.024 


1 PTVI 1 1 

UTVL11 :=HTVL1 1 - — ■ Llltil 
RJ RHOVL11 


Thus UTVL1 1 has uncertainty due to the 
uncertainty from the property routine. 


UTVL11 =-26.17 

**this value for UTVL1 1 is slightly 
different that the value in the output 
file used of -25.93668 


e RHOVLll 


1 PTVL11 
RJ RHOVLll 2 


U 


UTVL11 : -V( 0 RHOVLll U RHOVLll) 2 
U UTVL1 1 =0 - 202 


Now rewriting the balance relation and determining the uncertainty 

BAL22 := ( WF10 + WMLK1 ) (HTVL1 1 - UTVL1 1 ) - WF9-(HTVL10 - UTVL1 1 ) 
BAL22 =-0.013 

e WF10 =HTVL11 - UTVL11 

e w9 : = -HTVL10 +• UTVL1 1 
0 UTVL1 1 =-WF10- WMLK1 + WF9 

U BAL22 := J( 0 WF1O' U WFlo) 2 -*- ( 0 WF9‘ U WF9) 2- '" ( 0 UTVLH' U UTVLll) 2 
U BAL22= 456253 


Now write the uncertainty as a percentage of the RHS of the balance relation" 


U% BAL22 - 


U BAL22 

( WF10 + WMLK1 ) (HTVL1 1 - UTVL1 1 ) 


U% BAL22 ~ 7 071 ‘° /o 
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Balance Relation 23: WF10=WLPFT 

* WF1 0 is the fuel line 1 0 flow and is discussed above 

* WLPFT is the flow through the Low Pressure Fuel Turbine 

and is calculated in TURB02 and uses the flow area map in TBMP04 


W= RHOIN- 5! ? AREA 

4 I RHOMAPGC (sFAREA CLEN 2 ) 

PTVL12 := 4388.501 PTIN:=PTVL12 

PTLTFD := 3288.016 PTOUT :=PTLTFD 

DP :=PTTN- PTOUT 
RHOIN = RHOVL12 


U ° /o XAREA =005 
U% RHOIN :=002 
U XAREA : ” U% XAREA' XAREA 
U RHOIN = U% RHOIN RH0IN 

u% wlpfx = ^ * 0// ° 

U WLPFT :=w LPFT U% WL p FT 
U WLPFT = 1914 

BAL23 := WLPFT - WF10 
BAL23 = 1* 10~ 5 

U BAL23 :=U WLPFT 
U BAL23 =L914 


DP 


XAREA 


W= JRHOIN- 

^ RHOMAP GC ( SFAREA- CLEN 2 ) 


U% 


WLPFT 


XAREA \ / 1 U RHOIN \ 

1 XAREA / + \2 RHOIN 


WLPFT =31.90515 

RHOVL12 :=1.311(T 3 
XAREA =1 


Now write the uncertainty as a percentage of the RHS of the balance relation" 


U% BAL23 " 


U BAL23 
WLPFT 


U ° /o BAL23 ~ 6 * 0/o 
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Balance Relation 24: WF1 0*(HTVL1 1 -UTVL1 2)+QDOTVL12=WLPFT*HTVL1 2 

* WF10, WLPFT discussed above 

* HTVL1 1 ,HTVL12 are iteration variables 

* UTVL12 is calculated using enthalpy, pressure, and density 

* QDOTVL12 is the heat transfer for the chamber wall cooling 

HTVL12 =832.1201 


HTVL12 and PTVL12 are iteration variables: 


PTVL12 =4388.501 


1 PTVI 12 

UTVL12 :=HTVL12 - - — LLlti±. 

RJ RHOVL12 

UTVL12 =473.368 

**this value for UTVL12 is slightly 

Thus UTVL1 1 has uncertainty due to the different that the value in the output 

uncertainty from the property routine. file used of 473.8838 


1 PTVL12 

e RHOVL12 : U RHOVL12 =006RHOVL12 

RHOVL12 


U UTVL12 := J( e RHOVL12 U RHOVL12) 2 
U UTVL12 =21 - 525 


Now look at the uncertainty in QDOTVL12, the chamber wall heat transfer. 

As discussed in balance relation14, the uncertainty in the heat transfer calculations is estimated by 
combining an uncertainty associated with the heat flux calculation and an uncertainty associated 
with the heat transfer area. The heat transfer subroutine, QH2 includes a second heat transfer 
correlation to account for coolant boiling and two-phase flow. Again, no information is available at 
this time to assess whether this heat flux model correctly models the actual physical phenomena, 
so an overall uncertainty estimate, based upon engineering judgment, of 25% will be used. 

The area through which the heat transfer is occurring should be easier to determine for the main 
combustion chamber, so an uncertainty of 10% will be used. 


AHOTVL12 := 3912 
QDOTVL12 := 20931.83 

U% AHOTVL12 :=01 
U% QHOTVL12 - Q2S 

U QDOTVL12 = Q D0TVL12 U% QDOTVL12 
U QDOTVL12 = 1067 * 1q4 


QDOTVL12=AHOTVL12Q H0TVL12 

Determining the uncertainty expression on a percentage basis: 

U% QDOTVL12 = J U% AHOTVL12 2 + U ° /o QHOTVL12 
U% QDOTVL12 =051 
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Rewriting the balance relation and propagating the uncertainties: 

BAL24 := WLPFT- (HTVL 12 - UTVL12) - WF10-(HTVL11 - UTVL12) - QDOTVL12 
BAL24 =-0.003 


® WLPFT - HTVL12 
©WF10 =-HTVLll +UTVL12 

0 UTVL12 :=WF1 ° 

® QDOTVL12 ;= ‘ 1 


U BAL24 


( e WLPFT' U WLPFT) 2 + (°WF10 U WFlo) 2+ ( 0 UTVL12 U UTVL12) 2 - 
-J + ( 0 QDOTVL12 U QD0TVL12) 


U BAL24 = 1 082 * 104 


Now write the uncertainty as a percentage of the RHS of the balance relation" 


U% BAL24 : 
U% BAL24 


- U BAL24 
WLPFT- HTVL 12 
= 40.769*% 
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Balance Relation 25: WLPFT=WF11+WF13+WFRPR 

* WLPFT, WF11 are discussed above 

* WF13 is the flowthrough fuel line 13 and is calculated in FLOWOO 

* WFRPR is the fuel repressurization flow and is set at 0.2 Ib/s 

RF1 1 : = 2.5510‘ 4 
U jyr] | =0.1 RF1 1 

RHOVL13 = 1.09- 10" 3 
U RHOVL13 = 006RHOVL13 

WFRPR =0.2 


Evaluating the partial derivatives and dividing by WF1 1 to obtain 
the uncertainty as a percentage. 


First, look at the calculation of WF1 1 and it’s uncertainty 

WF1 1 =13.43799 


WF1 1= AP 


RHOVL13 


RF1 1 

Consider uncertainties in the resistance and 
the density from the property routine. 


U WF1 ( e RFl l' U RF1 1) 2+ ( e RHOVL13' U RHOVL13) 2 


U %WF11 


I 1 U RFll \ 2 | / 1 u RHOVL13 \ 2 
4\ 2 RF1 1 / \2 RHOVL13 / 


U% WF11 =5.831*% 

U WF1 1 := WFH U% WF u 
U WF11 = 0 - 784 


Now, look at the calculation of WF1 3 and it’s uncertainty 

RHOVL13 
RF13 

Consider uncertainties in the resistance and 
the density from the property routine. 

U WF13 ~J ( e RF13' U RF13) 2+ ( e RHOVL13 U RHOVL13) 2 

Evaluating the partial derivatives and dividing by WF13 to obtain 
the uncertainty as a percentage. 

TT0/ / 1 U RF13\ 2 / 1 U RHOVL13\ 2 

0WF13 J\ 2 RF13 / + \2 RHOVL13 / 

U% Y/p 1 3 = 0.058 
U WF13 : -WF13-U% wl3 


WF13= 


WF13 := 18.26689 

RF13 := 1.38 10" 4 
u Rpi 3 =0.1 -RF13 

RHOVL13 := 1.09- 10“ 3 
U RHOVL13 ~ 0.06 RHOVL1 3 


U \vpi3 - 1.065 


Including a 10% uncertainty for the fuel repress flow. 


U wprpr ’-0.1 0- WFRPR 
U WFRPR = 0 02 
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Rewriting the balance relation and determining the uncertainty 


BAL25 = WF11 + WF13 + WFRPR- WLPFT 
BAL25 = -2.7*lCf 4 


U 


2 2 2 2 
BAL25 : "^ U WF11 +U WF13 +u WLPFT + U WFRPR 

U BAL25 = 2 - 327 


Now write the uncertainty as a percentage of the RHS of the balance relation" 


U% BAL25 : 
U% BAL25 


. U BAL25 

|WF11 + WF13 + WFRPR) 

= 7.293*% 
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Balance Relation 26: WLPFTHTLTFD+QDOTVL1 3=(WF1 1 +WF1 3+WFRPR)*HTVL1 3 

* WLPFT, WF11 ,WF1 3.WFRPR are discussed above 

* HTLTFD, HTVL1 3 are iteration variables 

* QDOTVL1 3 is defined as zero 


Rewriting the balance relation and determining the uncertainty 

QDOTVL13 =0 
HTVL13 := 781 .1489 
HTLTFD : = 781 . 1 362 


BAL26 =(WF11 + WF13 + WFRPR) HTVL13 - WLPFT HTLTFD- QDOTVL13 
BAL26 =0.194 


® WF1 1 _ HTVL13 
e WF13 "HTVL13 
e WFRPR =HTVL13 
©WLPFT “ - HTLTFD 


U BAL26 “ J( e WFll' U WFll) + ( e WF13 U WF13) + ( e WFRPR u WFRP r) 2+ (®WLPFT u WLPFt) 2 


U BAL26 — 1 -817* 10 3 


Now write the uncertainty as a percentage of the RHS of the balance relation" 

U% - UbaL26 


BAL26 |(WF11 + WF13 + WFRPR)-HTVL13| 
U% BAL26 =7292 * % 
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Balance Relation 27: WF1 1 +WF1 3=WFSLV+WF1 5 

* WF11.WF13 are discussed above 

* WFSLV is the flow through the fuel line to the main fuel injector 

* WF15 is the flowthrough fuel line 15. 

First, look at the calculation of WFSLV and it’s uncertainty 

WFSLVs AP • RH0VL16 WFSLV := 4.179799 

V RFSLV 


Consider uncertainties in the resistance and 
the density from the property routine. 


RFSLV =23.5 10 -4 
U RFSLV -0.1RFSLV 

RHOVL16 := 9.97-10 -4 
U RHOVL16 =006RHOVL16 


U 


WFSLV 




RFSLV U RFSLV 


)M' 


0 


RHOVL16 U RHOVL16 


Evaluating the partial derivatives and dividing by WFSLV to obtain 
the uncertainty as a percentage. 


U% WFSLV : = 


i 1 U RFSLv \ 2 / 1 U RHOVL16 \ 2 
4 \ 2 RFSLV / + \2 RHOVL16 / 


U% y/FSLV ~ ~ C// ° 


U WFSLV : -WFSLV U% ^psly 
U WFSLV = 0 244 


Now, look at the calculation of WF15 and it's uncertainty 


WF15= 


AP 


RHOVL16 

RF15 


WF15 - 27.52525 


Consider uncertainties in the resistance and 
the density from the property routine. 


U WF15“J ( 0 RF15 U RF15) 2+ ( e RHOVL16 U RHOVL16) 2 

Evaluating the partial derivatives and dividing by WF13 to obtain 
the uncertainty as a percentage. 


RF15 : = 3.16- 10" 4 
Urp^^o.irfis 


no/ [j 1 U RF15\ 2 /l u RHOVL16\ 2 

0WF15 J\ 2 RF15 / + \2 RHOVL16 / 

U%wfi5 =5.831*% 

u WF15 : =WF15-U% WF15 
U YVF15 = 


A5-43 



Rewriting the balance relation and determining the uncertainty 
BAL27 := WFSLV + WF15 - WF1 1 - WF13 

BAL27 = 1.69*1 (f 4 

( 2 2 2 2 
U BAL27 := V U WFSLV +U WF15 + U WF11 +U WF13 

u Ft AT. 7 7 =2.094 


Now write the uncertainty as a percentage of the RHS of the balance relation: 


U% BAL27 : 
U% BAL27 


U BAL27 

’ | WFSLV + WF15| 
= 6.604*% 
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Balance Relation 28: (WF11+WF13)*HTVL13+QDOTOMCI+QDOTFMCI=(WFSLV+WF15)*HTVL16 

* WF1 1 ,WF1 3,WFSLV,WF1 5 are discussed previously 

* HTVL13.WHVL16 are iteration variables 

* QDOTOMCI, QDOTFMCI are manifold cooling heat transfers 

From the configuration file 

TKMCO =0.00167 

TKMCF:= 0.001 31 
AHTMCO =728.0 

AHTMCF =872.0 

QDOTFMCI = TKMCO AHTMCO-(TTMFI - TTVL13) 

QDOTOMCO : = - QDOTFMCI 
QDOTOMCI : =TKMCF- AHTMCF -(TTMFI- TTVL13)- 

QDOTOMCO := QDOTOMCI 

QDOTOMCI = 1. 363- 10 3 
QDOTFMCI = 1.976- 10 3 

Using the same logic and justification as discussed in balance relations 16 & 24, the heat transfer 
uncertainty is estimated by estimating an uncertainty for the heat transfer area and an uncertainty 
for the heat flux. 

QDOTOMCI= AHTMCF- QDOT qmCI 

Determining the uncertainty expression on a percentage basis: U% AHTMCF = 0 1 

U% QDOTOMCI := J U% AHTMCF 2 + U% QOMCI U/o QOMCI 

U% QDOTOMCI = 45826 * 0/o 

U QDOTOMCI := Q DOTOMCI U ° /o QDOTOMCI 
U QDOTOMCI =624 - 806 
QDOTFMCI= AHTMCO QDOT FMCI 
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HTVL16 := 886.4781 
TTMFI := 1549.379 

TTVL13 := 267.5007 



Determining the uncertainty expression on a percentage basis: U% AHTMCO := 0. 1 

U% QDOTFMCI : ' J U% AHTMCO 2 + U% QFMCI U /0 QFMCI °‘ 2 ° 

u% QDOTFMCI = 45826 * 0/o 


U QDOTFMCI - QDOTFMCI- U% q DO tfmCI 
U QDOTFMCI = 905 623 

Rewriting the balance relation and determining the uncertainty: 

BAL28 := (WFSLV + WF15) HTVL16 - (WF1 1 + WF13) HTVL13 - QDOTOMCI - QDOTFMCI 

BAL28 =-0.071 


0 WFSLV =HTVL16 
e WF15 =HTVL16 
6 WF1 1 = - HTVL13 
© WF13 = ■ HTVL13 
0 QDOTOMCI = - ! 

0 QDOTFMCI := - 1 


U BAL28 


( 0 WFSLV' U WFSLv) 2 -*- ( e WF15' U WFls) 2 + ( e WFll' U WFll) 2+ ( 0 WF13' U WF13) 2 - 
^l + ( @ QDOTOMCI U QDOTOMCl) + ( e QDOTFMCI 13 QDOTFMCl) 


U BAL28 =2085 * 103 


Now write the uncertainty as a percentage of the RHS of the balance relation: 


U% BAL28 " 
U% BAL28 = 


U BAL28 

| (WFSLV + WF15)HTVL16 

: 7.419-% 
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Balance Relation 29: WF7=WFFPB+WFOPB+WOTC 

* WF7 is discussed above 

* WFFPB is the fuel flow through the line to the fuel prebumer 

* WFOPB is the fuel flow through the line to the oxidizer preburner 

RFFPB := 1.243- 10' 4 
U rffpb =0.1 RFFPB 

RHOPBSF :=9.97-10" 4 
U rhoPBSF " = 0 06- RHOPBSF 

the density from the property routine. 


* WOTC \s the fuel flow through the oxidizer cooling line 

First, look at the calculation of WFFPB and it's uncertainty 

WFFPB= L ™°1 BSF WFFPB := 75.21584 

<\/ RFFPB 

Consider uncertainties in the resistance and 


U 


S J( € 


WFFPB = J ( 0 RFFPB’ U RFFPB 


)V(e 


RHOPBSF U RHOPBSF) 


Evaluating the partial derivatives and dividing by WFSLV to 
the uncertainty as a percentage. 


U% 


WFFPB ' 


l 1 U RFFPB 
1 2 RFFPB 


/ 1 U RHOPBSF 
+ \2 RHOPBSF 


2 


obtain 


U% wpppB - 5.83 1 *% 


U WFFPB WFFPB- U% WFFPB 
U WFFPB =4 - 386 


Now, look at the calculation of WFOPB and ifs uncertainty 

WFOPB= L P- RHQPBSF WFOPB : = 34.52161 

RFOPB 

Consider uncertainties in the resistance and 
the density from the property routine. 


RFOPB := 5.383- 10' 4 
U rpoPB = ^ RFOPB 

RHOPBSF 1.85- 10“ 3 
U RHOPBSF - 0-06- RHOPBSF 


U WFOPB = a/ ( e RFOPB U RFOPB) 2 + ( e RHOPBSF' U RHOPBSf) 2 

Evaluating the partial derivatives and dividing by WFSLV to obtain 
the uncertainty as a percentage. 


TT0/ _ / i u rfopb\ 2 1 1 u rhopbsf\ 2 

u/o WFOPB" — * + — 

4\ 2 RFOPB / \2 RHOPBSF j 

U% wpoPB = 5.831 •% 

U WFOPB = WFOPB- U% wfopb 
^ WFOPB ~ 2.01 3 
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Looking at the flow and uncertainty calculations for the oxidizer turbine cooling line, WOTC 

WOTC=-J AP- RHOPBSF- CFOTC 2 FLOCON WOTC : = 0.929067 

Consider uncertainties in the resistance and CFOTC = 0.661 

the density from the property routine. 

U CFOTC ~ 01 CFOTC 

U WOTC rj ( e CFOTC U CFOTC) 2 + RHOPBSF' U RHOPBSf) 2 


Evaluating the partial derivatives and dividing by WFIG to obtain the uncertainty as a percentage. 
U%woTC 

U% ^qtc = 10.44 •% 

U WOTC : = W0TC U % WOTC 

U WOTC = 0097 

Rewriting the balance relation and determining the uncertainty 

BAL29 : = WFFPB + WFOPB + WOTC - WF7 

BAL29 =-4.83*10 4 

I 2 2 2 2 

U BAL29 := V U WFFPB +U WFOPB +U WOTC +- U WF7 
U BAL29 =8058 

Now write the uncertainty as a percentage of the RHS of the balance relation: 



U% BAL29 " 


U BAL29 

WFFPB + WFOPB +• WOTC 


U% BAL29 =7 - 282 ‘ % 
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Balance Relation 30: WF7*(HTVL08-UTPBSF)=(WFFPB+WFOPB+WOTC)*(HTPBSF-UTPBSF) 

* WF7, WFFPB.WFOPB, WOTC are discussed above 

* HTVL08, HTPBSF are iteration variables 

* UTPBSF is determined from HTPBSF, PTPBSF, RHOPBSF 


HTVL12 and PTVL12 are iteration variables: 


HTPBSF := 551.7749 


i PTPRSF 

UTPBSF := HTPBSF - - — — - 

RJ RHOPBSF 


Thus UTPBSF has uncertainty due to the 
uncertainty from the property routine. 


PTPBSF =5301.58 


UTPBSF =244.884 

**this value for UTPBSF is slightly 
different that the value in the output 
file used of 244.5773 


0 RHOPBSF - 


1 PTPBSF 
RJ RHOPBSF 2 


U UTPBSF -*j( 0 RHOPBSF U RHOPBSf)' 


U UTPBSF - 18413 


Rewriting the balance relation and determining the uncertainty 

BAL30 := ( WFFPB+ WFOPB + WOTC) (HTPBSF - UTPBSF) - WF7-(HTVL08 - UTPBSF) 
BAL30 =-0.137 

0 y /FFP B * r BTPBSF — UTPBSF 
0 WFOPB ^HTPBSF- UTPBSF 
0 WOTC • ~ HTPBSF — UTPBSF 
0 WF7 : = . HTVL08 -h UTPBSF 
0 UTPBSF —WFFPB- WFOPB - WOTC+ WF7 

U BAL30 := ( 0 WFFPB U WFFPb) 2+ ( e WFOPB U WFOPb) 2 + ( e WOTC U WOTc) 2 • 
4 + ( 0 WF7 U WF7) + ( 6 UTPBSF' U UTPBSf) 

^ BAL30 = 2.473* 10 3 


Now write the uncertainty as a percentage of the RHS of the balance relation: 

uo/o , = U BAL30 

°BAL30 j (WFFPB+ WFOPB + WOTC)- (HTPBSF- UTPBSF)| 

U% BAL30 =7 - 282 ‘ % 
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Balance Relation 31 : WPDIPOGO=WPDOPOGO 

recirc line continuity WPDIPOGO 

BAL3 1 : = WPDOPOGO - WPDIPOGO WPDOPOGO 

BAL31 =0 

U WPDOPOGO = 0.10- WPDOPOGO 
U WPDIPOGO : = °- 1 ' WPDIPOGO 


U BAL3 1 : = WPDOPOGO + U WPDIPOGO 2 

U BAL31 = 0 - 385 

Now write the uncertainty as a percentage of the RHS of the balance relation: 
U% U BAL31 


BAL31 | WPDOPOGO | 
U ° /o BAL31 = 14142 * % 


= 2.723936 
= 2.723936 
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Balance Relation 32: QINPOGO=QOUTPOGO 
energy 

BAL32 = QOUTPOGO - QINPOGO 
BAL32 =0.005 

u qoutpogo : =o.ioqoutpogo 

U QINPOGO : = 0.1 QINPOGO 


QINPOGO =21.79906 
QOUTPOGO =21.80359 


[ 2 2 
U BAL32 = V U QOUTPOGO + U QINPOGO 

U BAL32 = 3083 

Now write the uncertainty as a percentage of the RHS of the balance relation: 

no/ - U BAL32 

“BAL32 | QOUTPOGO | 

U% BAL32 = 14141 * % 
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Balance Relation 33: WLXIPOGO=WLXOPOGO 
LOX continuity 

BAL33 := WLXOPOGO - WLXIPOGO 
BAL33 =0 

U WLXOPOGO =0- 10- ■ WLXOPOGO 
U WLXIPOGO :=01 WLXIP 000 


WLXIPOGO 

WLXOPOGO 


2 2 
U BAL33 := V U WLXOPOGO + U WLXIPOGO 

U BAL33 = 0 - 339 

Now write the uncertainty as a percentage of the RHS of the balance relation: 

TT0/ .. U BAL33 

°BAL33 | wlxoPOGO | 

U% BAL33 = 14 142 - % 


= 2.398129 
= 2.398129 
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Balance Relation 34: WGXIPOGO=WGXOPOGO 
GOX continuity 

BAL34 := WGXOPOGO- WGXIPOGO 
BAL34 = 0 

U WGXOPOGO = 0 1 °' WGXOPOGO 
U WGXIPOGO =0.1' WGXIPOGO 


WGXIPOGO -0.3511104 
WGXOPOGO =0.351 1104 


U BAL34 : = J U WGXOPOGO 2 + U WGXIPOGO 2 


U BAL34 _0 05 


Now write the uncertainty as a percentage of the RHS of the balance relation: 
U BAL34 


U% 


BAL34 


| WLXOPOGO| 
U% BAL34= 2071 - % 
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Balance Relation 35: WHEIPOGO=WHEOPOGO 
He continuity 

BAL35 : = WHEOPOGO - WHEIPOGO 
BAL35 =3.222*10 f 9 
U WHEOPOGO =0.10- WHEOPOGO 

U WHEIPOGO : = 0.1 WHEIPOGO 


WHEIPOGO =0 
WHEOPOGO = 8.7810" 9 


/ 2 2 
U BAL35 = 'J U WHEOPOGO + U WHEIPOGO 

U BAL35 = 3.222* 10~ 10 

Now write the uncertainty as a percentage of the RHS of the balance relation: 

tto/. U BAL35 

U/oqaL 35 ' — i T 

| WHEOPOGO | 

U% BAL35= 10 * % 
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Balance Relation 36: W01+WRIV=W02 

* WOI is the flow through oxidizer line 1 

* WRIV is the recirculation flow, somehow related to the pogo 

* W02 is the flow through oxidizer line 2 


First, look at the calculation of WOI and ifs uncertainty 
RHOOTNK 
ROl 

Consider uncertainties in the resistance and 
the density from the property routine. 



WOI =903.2136 
RHOOTNK =4.09- 10' 2 

U RHOOTNK := 0 002 RHOOTNK 
ROl :=1.0010' 6 
U roi :=0.1ROl 


U WO 1“ J ( 0 RO 1 U RO 1 ) 2 - 1 - ( 0 RHOOTNK’ U RHOOTNk) 2 

Evaluating the partial derivatives and dividing by WOI to obtain 
the uncertainty as a percentage. 


u % woi 


/ 1 U ROl \ 2 1 U RHOOTNK ^ 2 
J 2 ROl / + \2 RHOOTNK / 


U% WQ1 = 5.001 •% 


U 


WOI 


:= WOI U% 


WOI 


U woi =45-17 


Now, look at the calculation of W02 and it's uncertainty 


W02= (AP 


RHOVL17 

R02 


Consider uncertainties in the resistance and 
the density from the property routine. 


U W02® J ( 6 R02' U R02) 2 + (® RHOVL17 U RHOVLn) 2 

Evaluating the partial derivatives and dividing by W02 to obtain 
the uncertainty as a percentage. 


W02:= 905.9374 
RHOVL17 =4.09- 10' 2 

U RHOVL17 =0 002 RHOOTNK 
R02 =2.23 10 -7 

U R02 : =0.1RO2 


U% W02 := 

U ^ 2 R02 / \2 RHOVL17 / 

U% W02 =5001 * 0/o 
U W02 : = W °2 U% W02 

U W02 = 45 306 
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Sufficient information to determine an uncertainty for WRIV was not found, so an uncertainty will be 
estimated as the same percentage uncertainty as determined for the majority of the other ducts. 


U% v • - 0-05 

U WRIV :=U% WRIV WRJV 


WRIV := 2.723989 


U WRIV ~ 0136 


Rewriting the balance relation and determining the uncertainty 


BAL36 : = W02- WOl - WRIV 


BAL36 =-1.89-10 


.-4 


U BAL36 := J U W02 2 + U WOl 2 +U WRIV 2 


U BAL36 = 63.976 


Now write the uncertainty as a percentage of the RHS of the balance relation: 
U BAL36 


U% 


BAL36 


W02 


U% BAL36 =7062 * % 
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Balance Relation 37: (W01+WRIV)*HT0TNK+QD0TVL17=W02*HTVL17 

* WOI ,WRIV,W02 are discussed above 

* HTOTNK, HTVL17 are iteration variables 

* QDOTVL1 7 is zeroed out in this version 


Rewriting the balance relation and determining the uncertainty 

BAL37 := W02HTVL17 - (WOI + WRIV) HTOTNK+ QDOTVL17 
BAL37 =-0.003 

0 WO2 :=HTVL17 
©WOI • ~ - HTOTNK 
0 wRiv ~ * HTOTNK 

U BAL37 := J ( 0 W02' U WO2) ^ (®W01 U WOl) 2+ ( e WRIV U WRiv) 2 
U BAL37 = 4.033* 10 3 


Now write the uncertainty as a percentage of the RHS of the balance relation: 


U% BAL37 
U% BAL37 


U BAL37 
’ | W02 HTVL17 

= 7.062’% 


QDOTVL17 =0 
HTOTNK =63.04315 
HTVL17 =63.04316 
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Balance Relation 38: W02=WLPOP 

* W02 is discussed above 

* WLPOP is the flow through the Low Pressure Oxidizer Turbine, and is an iteration variable 
Rewriting the balance relation and determining the uncertainty: 


BAL38 := WLPOP -W02 
BAL38 = 1*1 0 _ 4 

U BAL38 ~ U W02 
U BAL38 = 45306 


WLPOP =905.9375 


Now write the uncertainty as a percentage of the RHS of the balance relation: 

TT0/ U BAL38 

°BAL38 | wlPop | 

U% BAL38 = 5 001 *° /o 
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Balance Relation 39: W02*HTVL17+QD0TVL18=WLP0P*HTVL18 

* W02.WLP0P are discussed above 

* HTVL17, HTVL18 are iteration variables 

* QDOTVL1 8 is heat transfer and is zeroed out in this model version 


QDOTVL18 :=0 
HTVL18 :=63. 04316 


Rewriting the balance relation and determining the uncertainty 

BAL39 :=WLPOP HTVL18- W02 HTVL17- QDOTVL18 
BAL39 =0.006 


0 WO2 --HTVL 17 


U 


BAL39 


4 


0 WO2' U W02)' 


U BAL39 = 2-856* 10 3 

Now write the uncertainty as a percentage of the RHS of the balance relation: 

us ° BAL39 


BAL39 |WLPOPHTVL18| 
U% BAL39 =5001 ‘ % 
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Balance Relation 40: WLP0P+WLP0T=W03 

* WLPOP is discussed previously 

* WLPOT is the flowthrough the Low Pressure Oxidizer Turbine (LPOT), 
uses TURB02 and the map in TBMP06 

* W03 is the flow through oxidizer line 3 


The LPOT flow is calculated using the equation 

SFW 

WLPOT=WMAP — 

CMASS-VSFRHO 


where WMAP is the value from the flow parameter vs pressure ratio and speed map 
SFW is the turbine area scale factor 
CMASS is a mass conversion parameter 
SFRHO is a density scale factor 


AREAMAP =1.0 
RHOMAP =0.040509 


WLPOT=WMAP 


I AREAMAP \ 
\AREA-CLEN 2 / 


1 

1 

CMASS 

1 

RHOMAP 

J 

RHOIN- CMASS- CLEN 3 


Consider uncertainties for the hardware characteristics WMAP and AREA, and for the inlet density 
density uncertainty from the physical property routine. 

Taking the partial derivatives with respect to each uncertainty source and dividing by WLPOT 


e% 


WMAP - 


1 _ 

WMAP 


-WMAP— 

2 [. 

00/0 RHODST 

WMAP 


AREAMAP 


AREA 


VcLEN- [ VcMASS- (^RHOMAP ^RHOIn) ] ] ] 


AREAMAP 


l AREA- CLEN 2 / CMASS 


1 


RHOMAP 


RHOIN- CMASS CLEN^ 


simplifying 


e% 


RHOIN (2-RHOIN) 
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-WMAP- 


AREAMAP 


0 % 


AREA- 


AREA 2 [^CLEN- (^CMASS- ^RHOMAp) ]] 


4 


RHOIN 


WMAP- 


AREAMAP 


AREA-CLEN 2 / CMASS 


1 


RHOMAP 


RHOIN- CMASS- CLEN 


simplifying 


RHOVL19 = 4.04-10' 


0 % 


AREA = 


-1 

AREA 


U RHOVL19 =0-002-RHOOTNK 
WLPOT := 178.7272 


0 WMAP = 1 

6 AREALPOT AREALPOT 
e RHOVL19 (2-RHOVL19) 


AREALPOT := 1. 


U% WMAP ;_0 05 
U ° /o AREALPOT : " 0 0 1 


U% RHOVL19 “ 


° RHOVL19 
RHOVL19 


The uncertainty expression then becomes 

no/ - I /t to/ \2 / TT0/ \2 / U% RHOVL19 

u/o WLPOT "K u/o WMAP) + (' u/o AREALPOT ) + 


U% WLPOT = 51 *° /o 
U WLPOT = WLPOT WLPOT 
U WLPOT = 9 115 


Now, look at the calculation of W03 and ifs uncertainty W03 : = 1084.665 

RHOVL19 
R03 

Consider uncertainties in the resistance and R03 :=2.23 10' 7 

the density from the property routine. 

u R03 :=01RO3 

U W03 = J (®R03' U R03) 2+ ( 0 RHOVL19 U RHOVL19) 2 


W03= 
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Evaluating the partial derivatives and dividing by W03 to obtain 
the uncertainty as a percentage. 


U % W03 


I 1 U RQ3 \ 2 /l U RHOVL19 \ 2 
J\ 2 R03 / \2 RHOVL19 / 


U% — 5.001 *% 
U W03 =W03- U %w03 
U W03 = 54 244 


Now, rewrite the balance relation and determine the uncertainty 

BAL40 := W03 - WLPOP - WLPOT 
BAL40 =3*1 0 _ 4 


U 


BAL40 : 'J U W03 2 + U 
U BAL40 = 55005 


WLPOT 2 


Now write the uncertainty as a percentage of the RHS of the balance relation: 


U% BAL40 : 
U% BAL40 


U BAL40 
|W03| 

= 5.071*% 
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Balance Relation 41: WLPOP*HTVL18+WLPOT*HTLTOD+QDOTLPOP=W03‘HTVL19 

* WLPOP, WLPOT, W03 are flowrates discussed above 

* HTVL18, HTVL19, HTLTOD are iteration variables 

* QDOTLPOP is the power required to drive the pump 

The required pump power is determined from the SPUMP module and the pump map uncertainties. 

The LPOP torque is calculated in SPUMP and uses a map for TAU vs PHI TAULPOP : = 1 .658089 

SNOL = 5.129*10 3 

TORQLPOP=- TAULPOP- TORQMLPOP-RHOVLO 1 • SNRADOL 2 SNRADOL = 537. 1 35 

TORQLPOP =-10521.79 

TRQPOL = | TORQLPOP | 

Considering an uncertainty associated with the torque coefficient vs flow coefficient map, and an 
uncertainty in density due to the property routine. 

U taulpop 3 0.05 TAULPOP 

U% TORQLPOP = 5004 * % 

U TORQLPOP = | u% TORQLPOP' TORQLPOP | 

U TORQLPOP =52652] 

U TRQPOL :=U TORQLPOP U TRQPOL = 526 521 
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Rewriting the balance relation and determining the uncertainty 

BAL41 := W03HTVL19 - WLPOPHTVL18- WLPOT- HTLTOD - QDOTLPOP 


HTLTOD : = 74.48289 
HTVL19 := 65.95411 


BAL41 =507.601 

0 W03 := HTVL19 
0 WLPOT : = ‘ HTLTOD 
0 QDOTLPOP =1 

U BAL41 : = J( 0 W03' U W03) 2 + ( 0 WLPOT' U WLPOt) 2 + ( 0 QDOTLPOP U QDOTLPOP) 2 
U BAL41 =3 - 642 * 1 ° 3 

Now write the uncertainty as a percentage of the RHS of the balance relation: 

ttoz. - U BAL41 

U/o rat 41 r 

1 |WQ3HTVL19| 

U% BAL41 “5.09*% 
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Balance Relation 42: W03+W06+WP0GL+WP0GG+W04=WHP0P 

* W03 is discussed previously 

* W04.W06 are the flows through oxidizer lines four and six, cal'd in PIPE01 

* WPOGG, WPOGL are pogo related flows, WPOGG is zeroed out, WPOGL has 
a value. A large uncertainty for WPOGL will be assigned. 

* WHPOP is the flowthrough the High Pressure Oxidizer Pump and is an iteration variable 

Looking at the flow and uncertainty calculations for oxidizer line four, W04 

WHPOP := 1094.372 
WPOGG =0 
WPOGL -- 2.38062 
W04:= 1.99514 
W06 := 10.08544 
CF04 := 0.221 
U cpo4 = 0-l-CFO4 

CF06 := 0.8416 
U CF06 =°1CF06 
RHOVL21 =4.05- 10* 2 
U RHOVL21 :::: 0.002-RHOVL2 1 


U W04 ;= W ° 4 U % W04 
U W04 = °' 2 

Assigning a 50% uncertainty to the WPOGL flowrate because no information is available. 

U wpoGL := 0-5- WPOGL 

U WPOGL = -1 - 19 

Looking at the flow and uncertainty calculations for oxidizer line six, W06 

RHOPBSO := 4.08- 10’ 2 
U RHOPBSO = 0.002- RHOPBSO 


U W06 = J ( 0 CF06 U CF06) 2 + ( 0 RHOPBSO U RHOPBSo) 2 


WO 6=V APRHOPBSO- CF06 2 FL0C0N 

Consider uncertainties in the resistance and 
the density from the property routine. 


W04=VAPRH0VL21 CF04 2 -FLOCON 

Consider uncertainties in the resistance and 
the density from the property routine. 

U W04 S J ( 0 CFO4 U CF04) 2+ ( 0 RHOVL21' U RHOVL2l) 2 

Evaluating the partial derivatives and dividing by W04 to obtain 
the uncertainty as a percentage. 

TT0 . ( U CF04\ 2 / 1 u RHOVL2l\ 2 

U/owo4 = 1 + — 

CF04 / \2 RHOVL21 / 

U% WO4 = 10 * % 
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Evaluating the partial derivatives and dividing by W06 to obtain 
the uncertainty as a percentage. 

T to/ / U CF 06\ 2 / 1 u RHOPBSo\ 2 

° W06 "J\ CF06 / + \2 RHOPBSO / 

U% W06 = 10 ‘° /o 

U W06 - W06 U% W q6 

U W06 = 1009 

Rewriting the balance relation and determining the uncertainty 

BAL42 := WHPOP - W03- W06- WPOGL- WPOGG- W04 
BAL42 = 0.007 

I 2 2 2 2 

U BAL42 := W03 + U W06 + U WPOGL + U W04 

U BAL42 = 54 267 


Now write the uncertainty as a percentage of the RHS of the balance relation 
U BAL42 


U% 


BAL42 


WHPOP 


U% BAL42 =4 - 959 * % 
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Balance Relation 43: 

W03*HTVL1 9+ W06*HTPBS0+ (WPOG L+ WPOGG)*HTVL20+ W04*HTVL2 1 +QDOTVL20 
=WHPOP*HTVL20 


* W03,W04,W06, WHPOP, WPOGL, WPOGG are discussed above 

* IHTVL19,HTVL20,HTVL21 .HTPBSO are iteration variables 

* QDOTVL20 is zeroed out in this version 


HTVL1 9 =65.954 
HTVL20 := 66.20525 


HTVL21 := 80.70924 
HTPBSO =90.34464 
QDOTVL20 =0 


BAL43 := WHPOP HTVL20 - W03HTVL19 - W06- HTPBSO- (WPOGL +■ WPOGG ) HTVL20 ... 
+ - W04 HTVL2 1 - QDOTVL20 


BAL43 =0.475 


0 WO3 =-HTVL19 
© W06 ~ ■ HTPBSO 
e WPOGL =-HTVL20 
© W04 := -HTVL21 


U BAL43 




I^WOS^WOj) 2 -*- ( e W06’ U 
+ ( 0 WPOGL' U WPOGL) +( 0 


W06) 2 - 

W04' U W04) 


2 


U BAL43 — 3.58-10 3 


Now write the uncertainty as a percentage of the RHS of the balance relation: 
U BAL43 

’ | WHPOP HT VL20 1 
= 4.941-% 


BAL43 

U ° /o BAL43 
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Balance Relation 44: WHP0P=W04+W05+WHX+WM0V+WPRBP+W0RPR 

* WHPOP is an iteration variable 

* W04.W05 are discussed above. 

* WHX is a flow to the POGO accumulator (?) No information is available in the configuration file 
upon which to calculate an uncertianty. Estimate an uncertainty of 25% 

* WMOV is the flow through the main oxidizer valve. 

* WPRBP is the flow through the Prebumer Pump, and is an iteration variable. 

* WORPR is the oxygen repressurization flow, and is a setpoint. 


Now, look at the calculation of WMOV and it's uncertainty 


WMOV* 


U RHOVL21 

Ar 

RMOV 


Consider uncertainties in the resistance and 
the density from the property routine. 


U WMOV® J ( e RMOV U RMOv) 2+ ( e RHOVL21 U RHOVL2l) 2 

Evaluating the partial derivatives and dividing by WMOV to obtain 
the uncertainty as a percentage. 

TT0 . I 1 U RMOv\ 2 /l u RHOVL2l\ 2 

0 WMOV 2 rj^ov j + RHQVL21 / 

U% WMOV = 005 
U WMOV ;= WMOV U% WMOV 
U WMOV = 40 2 13 


WMOV =804.09 
RMOV =7.0- 10' 5 
U RMOV =0.1 RMOV 

WHX =0.3511104 
WPRBP := 108.1089 
WORPR := 1.1 
W05 := 178.7268 


Now, look at the calculation of W05 and it's uncertainty 

W05=f^^i 

V R05 

Consider uncertainties in the resistance and 
the density from the property routine. 


R05 =0.86- 10' 4 
U rq5 :=0.1RO5 


U WOSrJ (® R05' U R05) 2+ ( 0 RHOVL21 U RHOVL2l) 2 

Evaluating the partial derivatives and dividing by WMOV to obtain 
the uncertainty as a percentage. 

U% W05 

u % wo5 =5.001 •% 

U W05 = W05 U% W05 
U W05 =8938 
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Include a 10% uncertainty for the oxygen repressurization flow 
U WORPR : = °- 1 °' WORPR 
u WORPR =011 

Include a 25% uncertainty for the POGO flow 

U WHX := °- 25 WHX 
U WHX = 0088 

Rewriting the balance relation and determining the uncertainty. 

BAL44 := W04 + W05 + WHX + WMOV +. WPRBP + WORPR - WHPOP 
BAL44 = -4.96* lO"" 5 


U 


12 2 2 2 
BAL44 := V U W04 - | - U W05 + U WHX +U WMOV + U 

U BAL44 = 41195 


WORPR 2 


Now write the uncertainty as a percentage of the RHS of the balance relation: 


U% BAL44 


U BAL44 

W04 + W05 + WHX + WMOV + WPRBP + WORPR| 


U% BAL44= 3 - 764 *% 
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Balance Relation 45: 

WHPOP*HTVL20+QDOTHPOP=(WO4+WO5+WHX+WMOV+WPRBP+WORPR)*HTVL21 

* WHPOP, W04,W05 I WHX,WM0V,WPRBP,W0RPR are discussed above 

* HTVL20,HTVL21 are iteration variables 

* QDOTHPOP is the power required by the High Pressure Oxidizer Pump, calculated in SPUMP 
The required pump power is determined from the SPUMP module and the pump map uncertainties. 

The HPOP torque is calculated in SPUMP and uses a map for TAU vs PHI 

TORQHPOP*- TAUHPOP- TORQMHPOP RHOVL2 1 • SNRADOH 2 TAUHPOP : = 0. 1 524294 

TORQHPOP : = - 51369.37 SNOH =2.755* 10 4 

SNRADOH = 2.885* 10 3 


Considering an uncertainty associated with the torque coefficient vs flow coefficient map, and an 
uncertainty in density due to the property routine. 


U% TORQHPOP := 


/ U TAUHPQP \ 2 / U RHOVL21 \ 2 
J [ TAUHPOP / + \ RHOVL21 / 


U% TORQHPOP “ 5 004 ’° /o 


U XAUHPOP TAUHPOP 


U TORQHPOP " | U% TORQHPOP T ORQHPOP | 
U TORQHPOP = 2 - 571 ‘ 1 0 3 


QDOTHPOP :=- TORQHPOP 


SNRADOH 


RJ 


QDOTHPOP = 1.587279* 10 4 
SNRADOH 


6 TORQHPOP 


RJ 


U QDOTHPOP := J( 0 TORQHPOP 11 TORQHPOP) ' 

U QDOTHPOP = 794274 

tto/ _ U QDOTHPOP 

u ^ QDOTHPOP QDOTHpop 

U% QDOTHPOP = 5 004 *° /o 
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Rewriting the balance relation and determining the uncertainty 

BAL45 := ( W04 -h W05 + WHX + WMOV + WPRBP +• WORPR) -HTVL21 - WHPOP HTVL20 - QDOTHPOP 
BAL45 =-0.03 


6 W04 -HTVL21 
6 W05 := HTVL21 
6 WHX =HTVL21 
® WMOV =HTVL21 
e WORPR :=HTVL21 
6 QDOTHPOP = ' ! 


U 


BAL45 - 


(® W04'° W04) 2 + " (® W05° W05) 2 + ( e WHX U WHx) 2+ ( e WMOV U WMOv) 2 - 
+ ( e WORPR' ° WORPr) + ( e QDOTHPOP 0 QDOTHPOp) 


°BAL45 -3 - 418 ‘ 1C)3 


Now write the uncertainty as a percentage of the RHS of the balance relation: 


U% BAL45 : 
U% BAL45 


°BAL45 

| ( W04 + W05 - 1 - WHX + WMOV + WPRBP +■ WORPR)HT VL2 1 1 

: 3.87*% 
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Balance Relation 46: W05=WLP0T 
* W05, WLPOT are discussed above 

Rewriting the balance relation and determining the uncertainty 


BAL46 := WLPOT- W05 
BAL46 = 4-10~ 4 
U 


2 2 
BAL46 := J U WLPOT + U W05 


U BAL46 = 12766 


Now write the uncertainty as a percentage of the RHS of the balance relation 


u 


u% 


BAL46 


BAL46 


| WLPOT| 
U% BAL46 =7 143 * % 
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Balance Relation 47: W05*HTVL21 +QDOTVL22=WLPOTH7VL22 

* W05, WLPOT are discussed above 

* HTVL21 , HTVL22 are iteration variables 

* QDOTVL22 is zeroed out in this model version 


HTVL22 =80.70924 
QDOTVL22 :=0 


BAL47 := WLPOT HTVL22 - W05 HTVL21 - QDOTVL22 
BAL47 =0.032 

0 W05 := -HTVL21 

U BAL47 ; -J( 0 WO5' U W05) 2 
U BAL47 = 721 - 389 


Now write the uncertainty as a percentage of the RHS of the balance relation: 


U% BAL47 : 
U ° /o BAL47 


U BAL47 
| WLPOT- HTVL22 

= 5.001*% 
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Balance Relation 48: WPRBP=W06+WOPOV+WFPOV 

* WPRBP is the flow through the Prebumer Pump and is an iteration variable 

* W06 is the flowthrough oxidizer line six, and calculated in PIPE01 

* WOPOV is the flow through the oxidizer prebumer oxidizer valve, cal'd in FLOWOO 

* WFPOV is the flow through the fuel prebumer fuel valve, cal'd in FLOWOO 


Looking at the flow and uncertainty calculations for oxidizer line six, W06 
W06=V APRHOPBSO CF06 2 FLOCON 


W06 := 10.08544 
CF06 = 0.8416 
RHOPBSO =0.041 


Consider uncertainties in the resistance and 
the density from the property routine. 


U wcxrj ( e CF06 ° CF06) + ( e RHOPBSO 0 RHOPBSo) 2 

Evaluating the partial derivatives and dividing by W04 to obtain the uncertainty as a percentage. 


u% 


W06 


'u 


CF06 


CF06 


’ (1 U RHOPBSO j 

+ 2 RHOPBSO / 


U% W06 " 10 *° /o 


U W06 - W06 U% WQ6 
U W06 = 1009 

Now, look at the calculation of WOPOV and it’s uncertainty 


WOPOVs AP 


RHOPBSO 


RLINOPOV 

Consider uncertainties in the resistance and 
the density from the property routine. 


WOPOV =28.04103 
RLINOPOV =0.0635 
u RLINOPOV ' °- 1 RLINOPOV 


U WOPOV® J ( e RLINOPOV' U RLINOPOV) + ( e RHOPBSO' 0 RHOPBSo) 2 


Evaluating the partial derivatives and dividing by WOPOV to obtain 
the uncertainty as a percentage. 


TT0/ / 1 °RLINOPOv\ 2 / 1 ° RHOPBSO \ 2 

0 WOPOV 2 rlkvjopov / + \2 RHOPBSO / 

U% WOPOV = * 0/ ^° 

° WOPOV - WOPOV- U% WOPOV 
° WOPOV =1 402 
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Now, look at the calculation of WFPOV and it’s uncertainty 


WFPOV =69.98241 


WFPOVs 


RHOPBSO 

RLENFPOV 


Consider uncertainties in the resistance and 
the density from the property routine. 


RLINFPOV =0.00905 
U RLINFPOV = °- 1 RLINFPOV 


U WFPOV® J ( e RLINFPOV U RLINFPOv) 2 + ( 0 RHOPBSO U RHOPBSo) 2 


Evaluating the partial derivatives and dividing by WOPOV to obtain 
the uncertainty as a percentage. 


u % WFPOV 


1 U RLINFPOV | 

2 RLINFPOV / 


/ 1 U RHOPBSO 
+ \2 RHOPBSO J 


U % WFPOV — 


U WFPOV := WFPOV- ^ U% WFPOV 


U WFPOV - 35 


Rewriting the balance relation and determining the uncertainty 


BAL48 := W06 + WOPOV + WFPOV - WPRBP 
BAL48 =-2*10 _5 


/ 2 2 2 
U BAL48 : = V U W06 +JJ WOPOV +u WFPOV 

U BAL48 =3 - 903 


Now write the uncertainty as a percentage of the RHS of the balance relation: 


U% BAL48 :: 
U% BAL48 


. U BAL48 

| W06 + WOPOV + WFPOV | 

= 3.61*% 
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Balance Relation 49: WPRBP*HTVL21+QD0TPRBP=(W06+W0P0V+WFP0V)‘HTPBS0 

* WPRBP, W06,W0P0V,WFP0V are discussed above 

* HTVL21 , HTPBSO are iteration variables 

* QDOTPRBP is the power required by the prebumer pump, cal'd in SPUMP. 


The required pump power is determined from the SPUMP module and the 
The PRBP torque is calculated in SPUMP and uses a map for TAU vs PHI 
TORQPRBPs- TAUPRBP TORQMHPOPRHOVLPBSO- SNRADOH 2 

TORQPRBP : = - 337 1.18 


pump map uncertainties. 

TAUPRBP =0.00992 
SNOH= 2.755- 10 4 
SNRADOH = 2.885- 10 3 


Considering an uncertainty associated with the torque coefficient vs flow coefficient map, and an 
uncertainty in density due to the property routine. 

U TAUPRBP = 0.05 T AUP RBP 

U TAUPRBP \ / U RHQPBSo \ 2 

TAUPRBP / + \ RHOPBSO / 

U% XORQPRBP = 5004 *° /o 

U TORQPRBP | U% TORQPRBP' T0R Q PRBP | 

U TORQPRBP = 168 694 


U% TORQPRBP = 


U TRQPOH ; “4 U TORQHPOP + U TORQPRBP 


U TRQPOH ~ 2 - 576 * 


10 J 


QDOTPRBP : = - TORQPRBP 


SNRADOH 


RJ 


QDOTPRBP = 1.0417-10 3 

SNRADOH 


RJ 


6 TORQPRBP := 

U QDOTPRBP := J ( e TORQPRBP' U TORQPRBp) 2 
U QDOTPRBP = 52 125 


U ° /o QDOTPRBP 
U% QDOTPRBP 


U QDOTPRBP 
QDOTPRBP 

= 5.004-% 
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Rewriting the balance relation and determining the uncertainty 

BAL49 := ( W06 + WOPOV + WFPOV)HTPBSO- WPRBP HTVL21 - QDOTPRBP 

BAL49 =-0.001 

0 yvo6 = ^ TPBSO 
9 WOPOV = HTPBSO 
0 WFPOV =HTPBSO 
6 QDOTPRBP 

U BAL49 := ( 0 W06' U W06) 2+ ( 0 WOPOV U WOPOv) 2 + ( e WFPOV U WFPOv ) 2 
( e QDOTPRBP U QDOTPRBP) 

U BAL49= 356437 

Now write the uncertainty as a percentage of the RHS of the balance relation: 

t to/. - U BAL49 

u/o 'RAT49 _ “i r 

|(W06+ WOPOV + WFPOV) HTPBSO| 

U% BAL49 = 3 - 649 * % 
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Balance Relation 50: WOOPB+WOPBF+WHE2=WHG2 

* WOOPB is the oxidizer to the OPB 

* WOPBF is the fuel to the OPB 

* WHE2 is the He dilutant to the OPB 

* WHG2 is the hot gas out from the OPB 


BAL50 : = WHG2 - WOOPB - WOPBF - WHE2 
BAL50 = —3* 10~ 5 


WOOPB 
WOPBF 
WHE2 := 
WHG2 := 


= 28.04103 
= 35.37508 
0 

63.41608 
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BAL51=WHG2*(GAMA(PTOPRB,OFROPRB,HFROPRB,TTHETK,TTOPBF,TTOPBI)-1)*TTOPRB 
-(WOOPB+WOPBF+WHE2)*(GAMA(PTOPRB,OFROPRB,HFROPRB,TTHETK, 
TTOPBF.TT OPBI)*TTCBST (PT OPRB ,TTHETK,TTOPBF,TTOPBI)-TTOPRB) 


Balance Relation 52: WOOPB=OFROPRB*(WOOPB+WOPBF+WHE2) 

BAL52=OFROPRB (WOOPB + WOPBF + WHE2) - WOOPB 


Balance Relation 53: WHE2=HFROPRB*(WOOPB+WOPBF+WHE2) 

BAL53=HFROPRB ( WOOPB + WOPBF + WHE2) - WHE2 


Balance Relation 54: 

HFRHTOD*WHPOT+HFROTBP*WHG4+HFRPBSF*WOTC=HFROSF*(WHPOT+WHG4+WOTC) 

BAL54=HFROSF • ( WHPOT + WHG4 + WOTC) - HFRHTOD- WHPOT - HFROTBP WHG4 - HFRPBSF- WOTC 


Balance Relation 55: 

OFRHT OD*WHPOT +OFROTBP*WHG4+OFRPBSF*WOTC=OFROSF*(WHPOT + WHG4+ WOT C) 

BAL55=OFROSF- ( WHPOT + WHG4+ WOTC) - OFRHTOD- WHPOT- OFROTBP WHG4- OFRPBSF WOTC 


Balance Relation 56: WHPOT+WHG4+WOTC=WHG6+WOLK 

BAL56= WHG6 + WOLK - WHPOT- WHG4- WOTC 
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Balance Relation 57: 

WHPOT*GAMAHTOD*(TTHT OD-TT OSF)+ WHG4*GAMAOTBP*(TT OTBP-TTOSF) 
+WOTC*GAMAPBSF‘(TTPBSF-TTOSF)+QDOTOSF*GAMAOSF/CPOSF 
=(WHG6+WOLK)*(GAMAOSF-1 )*TTOSF 


BAL57s( WHG6 + WOLK) • ( GAMAOSF - l) TTOSF- WHPOT- GAMAHTOD- ( TTHTOD - TTOSF) . 
+ - WHG4- GAMAOTBP- ( TT OTBP - TTOSF) - WOTC-GAMAPBSF (TTPBSF - TTOSF) ... 


+ - 


QDOTOSF- 


GAMAOSF 


CPOSF 


Balance Relation 58: WOFPB+WFPBF+WHE1=WHG1 

BAL58=WHG1 - WOFPB - WFPBF- WHE1 


Balance Relation 59: 

(WOFPB+WFPBF+WHE1)*(GAMA(PTFPRB,OFRFPRB,HFRFPRB, 

TTHETK,TTFPBF,TTFPBI)*TTCBST(PTFPRB,OFRFPRB,HFRFPRB,TTHETK,TTFPBF, 

TTFPBI)-TTFPRB)=WHG1*(GAMA(PTFPRB,OFRFPRB,HFRFPRB,TTHETK,TTFPBF, 

TTFPBI)-1)*TTFPRB 


Balance Relation 60: WOFPB=OFRFPRB*(WOFPB+WFPBF+WHE1) 

BAL60=OFRFPRB ( WOFPB + WFPBF + WHE1 ) - WOFPB 
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Balance Relation 61: WHE1+HFRFPRB*(WOFPB+WFPBF+WHE1) 


BAL61=HFRFPRB(WOFPB+ WFPBF + WHE1)- WHE1 


Balance Relation 62: 

HFRHTFD‘WHPFT+HFRFTBP*WHG3+HFRVL03*WFTC=HFRFSF*(WHPFT+WHG3+WFTC) 

BAL62=HFRFSF-( WHPFT + WHG3 + WFTC) - HFRHTFD- WHPFT - HFRFTBP- WHG3 - HFRVL03- WFTC 


Balance Relation 63: OFRHTFD*WHPFT+OFRFTBP*WHG3+OFRVL03*WFTC 

=OFRFSF*(WHPFT +WHG3+WFTC) 

BAL63=OFRFSF-( WHPFT + WHG3 + WFTC) - OFRHTFD WHPFT - OFRFTBP- WHG3 - OFRVL03- WFTC 


Balance Relation 64: WHPFT+WHG3+WFTC=WHG5 

BAL64=WHG5 - WHPFT - WHG3 - WFTC 


Balance Relation 65: 

WHPFT*GAMAHTFD*(TTHTFD-TTFSF)+WHG3*GAMAFTBP*(TTFTBP-TTFSF) 
+WFTC*GAMAVL03*(TTVL03-TTFSF)+QDOTFSF*GAMAFSF/CPFSF 
= WHG5* (GAMAFSF-1 )*TTFSF 
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Balance Relation 66: 

HFRFSF*WHG5+HFROSF*WHG6+HFRVL16*WFSLV=HFRMFI*(WHG5+WHG6+WFSLV) 

BAL66=HFRMFI(WHG5+ WHG6+ WFSLV) - HFRFSF WHG5 - HFROSF WHG6- HFRVL 16- WFSLV 


Balance Relation 67: OFRFSF*WHG5+OFROSF*OFROSF*WHG6+OFRVL16*WFSLV= 

OFRMFI‘(WHG5+WHG6+ WFSLV) 

BAL67=OFRMFI(WHG5+ WHG6+ WFSLV) - OFRFSF WHG5 + OFROSF WHG6 - OFRVL 1 6- WFSLV 


Balance Relation 68: WHG5+WHG6+WFSLV=WFINJ 

BAL68=WFINJ - WHG5 - WHG6 - WFSLV 


Balance Relation 69: 

WHG5*GAMAFSP(TTFSF-TTMFI)+WHG6*GAMAOSF*(TTOSF-TTMFI) 
+WFSLV*GAMAVL1 6‘(TTVL1 6-TTMFI)+(QDOTFMCO+QDOTOMCO)*GAMAMFI/CPMFI 
=WFIN J*(GAMAMFIR-1 )*TTMFI 
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Balance Relation 70: TRQPFLTRQNFL 

* TRQPFL is the torque required by the Low Pressure Fuel Pump 

* TRQNFL is the torque produced by the Low Pressure Fuel Turbine 
Both routines use speed as an iteration variable 


TAULPFP := 1.971925 

The LPFP torque is calculated in SPUMP and uses a map forTAU vs PHI SNFL = 1.379- 10 4 
TORQLPFP=- TAULPFP- TORQMLPFP- RHOVLO 1 • SNRADFL 2 SNRADFL = 1 .444- 1 0 3 

TORQLPFP =-10521.79 

TRQPFL — TORQLPFP 

Considering an uncertainty associated with the torque coefficient vs flow coefficient map, and an 
uncertainty in density due to the property routine. 

^ TAULPFP ■ = 0.05- TAULPFP 

U TAULPFP \ 2 / u RHQVL0l \ 2 
, TAULPFP / + \ RHOVLO 1 / 

U% TORQLPFP = 005 


U% TORQLPFP : 


U TORQLPFP = | u% TORQLPFP TOR Q LPFP | 

U TORQLPFP = 526 195 
U TRQPFL := U TORQLPFP 
U TRQPFL =526 - 195 

Now consider the torque calculation for the LPFT, modules TURB02 and TBMP04. 

Uncertainties are introduced through the turbine flow, WLPFT, as discussed in balance relation 23, 
and the turbine efficiency map. The DH is determined in TBMP04 by calling the H2 property 
routine for the exit pressure (an iteration variable) and inlet entropy, since this value is so closely 
related to the iteration variable no uncertainty will be assigne to DHLPFT. 


A 2.5% uncertainty for ETALPFT is being assigned based upon engineering judgment since no 
information was obtained. 


TORQLPFT=DHLPFT- WLPFT- — 

SNRADFL 


U#/o ETALPFT -° 025 
TORQLPFT := 10521.78 


DHLPFT=ETALPFT- DHIDLLPFT 


TORQLPFT =ETALPFT DHDDLLPFT • WLPFT • — 

SNRADFL 


TRQNFL := TORQLPFT 
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U% TORQLPFT : -J u% ETALPFT 2 + u% WLPFT 2 
u % TORQLPFT = 0 065 
U TORQLPFT = TORQLPFT- U% TORQLPFT 
U TORQLPFT = 683 916 
U TRQNFL =U TORQLPFT 
U TRQNFL = 683 916 

Rewriting the balance relation and determining the uncertainty 

BAL70 : = TRQNFL - TRQPFL 
BAL70 =-0.01 

U BAL70 : = TRQNFL 2 + U TRQPFL 2 
U BAL70 = 862915 


Now write the uncertainty as a percentage of the RHS of the balance relation: 


U% BAL70 - 


U BAL70 

|trqnfl| 


U% BAL70 “ 8 201 *° /o 
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Balance Relation 71: TRQPFH:TRQNFH 

* TRQPFH is the torque required by the High Pressure Fuel Pump 

* TRQNFH is the torque produced by the High Pressure Fuel Turbine 
Both routines use speed as an iteration variable 


The HPFP torque is calculated in SPUMP and uses a map for TAU vs PHI 
TORQHPFPss- TAUHPFP- TORQMHPFP RHO VL03 SNRADFH 2 

TORQHPFP : = - 107789.8 

TRQPFH := |TORQHPFP| 


TAUHPFP =3.003558 
SNFH= 3.366- 10 4 
SNRADFH = 3.524- 10 3 


Considering an uncertainty associated with the torque coefficient vs flow coefficient map, and an 
uncertainty in density due to the property routine. 


U% 


TORQHPFP := 


/ U TAUHPFP \ 2 / u RHOVLQ3 
TAUHPFP j + \ RHOVL03 


U TAUHPFP - TAUHPFP 


U% TORQHPFP " 0 078 


U TORQHPFP | U% TORQHPFP T0R QHPFP | 

U TORQHPFP = 8.419- 10 3 

U TRQPFH :=U TORQHPFP U TRQPFH = 8 4 190 ° 3 

Now consider the torque calculation for the HPFT, modules TURB01 and TBMP03. 

Uncertainties are introduced through the turbine flow, WHPFT, as discussed in balance relation 64, 
and the turbine efficiency map. The DH is calculated in TBMP03 as a function of the inlet and exit 
pressures, the inlet temperature, and the gas constant. The pressures are iteration variables, 
since DH is so closely related to the iteration variable no uncertainty will be assigne to DHLPFT. 

A 5% uncertainty for ETAHPFT is being assigned based upon engineering judgment since no 
information was obtained. 


T ORQHPFT =DHHPFT WHPFT — 

SNRADFH 


U% ETAHPFT :=005 
TORQHPFT := 107789.8 


DHHPFT = ET AHPFT • DHIDLHPFT 


T ORQHPFT =ET AHPFT • DHIDLHPFT • WHPFT — 

SNRADFH 


TRQNFH = TORQHPFT 


Since this balance relation is 
being worked prior to working #64 
uncertainties for WHPFT is 
being estimated 
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U% TORQHPFT • = Ju% ETAHPFT 2 + U% WHPFT 2 U /0 WHPFT U /0 WLPFT 

U% TORQHPFT =0.078 
U TORQHPFT = TORQHPFT U% TORQHPFT 

U TORQHPFT = 8-419* 10 3 
U TRQNFH :=U TORQHPFT 

u trqnfh = 84 19 * 1C|3 

Rewriting the balance relation and determining the uncertainty 
BAL71 := TRQNFH- TRQPFH 
BAL71 =0 

U BAL7 1 : = J U TRQNFH 2 +" U TRQPFH 2 
U BAL71 =1191 ‘ 104 


Now write the uncertainty as a percentage of the RHS of the balance relation: 


U% BAL71 

U% BAL71 


U BAL71 
~ | TRQNFH | 

= 11.045 •% 
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Balance Relation 72: TRQPOLTRQNOL 

* TRQPOL is the torque required by the Low Pressure Oxidizer Pump 

* TRQNOL is the torque produced by the Low Pressure Oxidizer Turbine 
Both routines use speed as an iteration variable 

The LPOP torque is calculated in SPUMP and uses a map for TAU vs PHI TAULPOP = 1 .658089 

SNOL =5.129*10 3 

TORQLPOP=- TAULPOP- TORQMLPOP- RHOVLO 1 SNRADOL 2 SNRADOL = 537.135 

TORQLPOP =-10521.79 
TRQPOL := | TORQLPOP | 


Considering an uncertainty associated with the torque coefficient vs flow coefficient map, and an 
uncertainty in density due to the property routine. 

I u TAULPOP = ° 05 ' TAULPOP 

TT0/ |/ U TAULPOP \ 2 / U RHOVLOl\ 2 

0 TORQLPOP J ( TAULPOp J + ( rhOVLO 1 / 

U% TORQLPOP = 005 

U TORQLPOP := | U% TORQLPOP' TORQLPOP j 

U TORQLPOP =5 26. 195 
U TRQPOL :=U TORQLPOP 

U TRQPOL = 526 195 

Now consider the torque calculation for the LPOT, modules TURB02 and TBMP06. 


Uncertainties are introduced through the turbine flow, WLPOT, as discussed in balance relation 46, 
and the the torque correlation in TBMP06. The ideal "delta-enthalpy” is determined in TBMP06 by 
calling the 02 property routine for the exit pressure (an iteration variable) and inlet entropy, since 
this value is so closely related to the iteration variable no uncertainty will be assigned to DHIDL. 


The torque produced by the LPOT is calculated in TBMP04 using a "torque correlation" equation 


where AOT1 ,BOT1 , ctqotl are defined constants 

PHIOT1 is a function of turbine speed and turbine flowrate 
and CLEN.CFORCE, CFORCE are unit conversion constants 


AOT1 =1.0998 
BOT1 := 0.1 6443 


PHIOT1 := 


SNRADOL 

(wLPOT- SFW- CMASS +M0’ 10 ) 


PfflOTl =3.005 


SFW := 1 


CTQOT1 := 1 


TORQLPOT := (AOT1 - BOT1 -PfflOTl )• WLPOT- 1 WLPOTl -CTQOT1 • — CMASS 2 

(CLEN- CFORCE) 


TORQLPOT 


AOT1 - BOT1- 


SNRADOL 


+ CTQOT1-- 


( WLPOT- SFW- CMASS + 110 10 ) 
SFW 


-WLPOT- 1 WLPOT| ... 


(CLEN- CFORCE) 


-CMASS 
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e AOTl = WLPOT- 1 WLPOT| • CTQOT1 ■ 


SFW 

(CLENCFORCE) 


•CMASS 2 


0 


BOT1 - 


-SNRADOL 


WLPOT- SFW • CMASS + - 


1 


10000000000 


WLPOT I WLPOTl CTQOT 1 - 


SI ‘ W —CMASS 2 


(CLENCFORCE) 


0 WLPOT ~ B0T1 ' 


SNRADOL 


WLPOT- SFW- CMASS - 


1 


-SFW 2 -CMASS 3 -WLP0T I WLPOTl ... 


+ CTQOT 1 | 2 


CLENCFORCE 


10000000000/ 
AOT1 - BOT1- 


SNRADOL 


+ CTQOT 1 


SFW 


WLPOT- SFW- CMASS + L_ 

\ 10000000000 


| WLPOTl ... 


(CLENCFORCE) 


•CMASS 


U AOT1 : =005AOT1 U B0T1 :=0.05-BOT1 


U TORQLPOT -J( e AOT1 U AOTl) 2 ^ ( e BOTl U BOTl) 2+ ( e WLPOT U WLPOt) 2 
U TORQLPOT = 3 - 395 * 1q3 


U% TORQLPOT " 
U% TORQLPOT 


U TORQLPOT 
TORQLPOT 

= 17.551 *% 


TORQLPOT = 1.934699* 10 4 
TRQNOL : = TORQLPOT 
U TRQNOL =U TORQLPOT 

Rewriting the balance relation and determining the uncertainty 
BAL72 := TRQNOL - TRQPOL 
BAL72 = 8.825* 10 3 

U BAL72 = TRQNOL 2 + U TRQPOL 2 
U BAL72 = 3 - 436 * 103 
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Now write the uncertainty as a percentage of the RHS of the balance relation 


U 


U% 


BAL72 


BAL72 
U% 


BAL72 


TRQNOLl 
= 17.76-% 
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Balance Relation 73: TRQPOH :TRQNOH 

* TRQPOH is the torque required by the High Pressure Oxidizer Pump and the Prebumer Pump 
TRQNOH is the torque produced by the High Pressure Oxidizer Turbine 
Both routines use speed as an iteration variable 


The HPOP torque is calculated in SPUMP and uses a map for TAU vs PHI 
TORQHPOP=- TAUHPOP- TORQMHPOP RHOVL2 1 • SNRADOH 2 

TORQHPOP :=- 51369.37 


TAUHPOP : = 0. 1524294 
SNOH= 2.755* 10 4 
SNRADOH = 2.885- 10 3 


Considering an uncertainty associated with the torque coefficient vs flow coefficient map, and an 
uncertainty in density due to the property routine. 


U% 


TORQHPOP := 


/ U TAUHPQP \ 2 
J\ TAUHPOP / + 


/ U RHOVL2l \ 2 
\ RHOVL21 / 


U TAUHPOP = 0 05- TAUHPOP 


U% TORQHPOP “ 005 


U TORQHPOP " | U% TORQHPOP' TORQHPOP 
U TORQHPOP = 2-57 1-1 0 3 


The PRBP torque is calculated in SPUMP and uses a map for TAU vs PHI 

TORQPRBP*- TAUPRBP- TORQMHPOPRHO VLPBSO SNRADOH 2 

TORQPRBP : = - 3371.18 


TAUPRBP := 0.00992 
SNOH= 2.755- 10 4 
SNRADOH = 2.885- 10 3 


Considering an uncertainty associated with the torque coefficient vs flow coefficient map, and an 
uncertainty in density due to the property routine. 


U% 


TORQPRBP - 


/ U TAUPRBP \ 2 / u RHOPBSO 
J \ TAUPRBP / + \ RHOPBSO 


U% TORQPRBP “ 5 004 *° /o 


U TAUPRBP =0 .05- TAUPRBP 


U TORQPRBP := | U% TORQPRBP' TORQPRBP j 
U TORQPRBP = 168 (594 
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TRQPOH := |TORQHPOP+ TORQPRBP 
TRQPOH = 5.474* 10 4 


U TRQPOH - 


2 

U TORQHPOP +U 


TORQPRBP 2 


U TRQPOH — 2.5 76 * 10 3 


Now consider the torque calculation for the HPOT, modules TURB01 and TBMP05. 

Uncertainties are introduced through the turbine flow, WHPOT, as discussed in balance relation 54, 
and the turbine efficiency map. The DH is calculated in TBMP05 as a function of the inlet and exit 
pressures, the inlet temperature, and the gas constant. The pressures are iteration variables, 
since DH is so closely related to the iteration variable no uncertainty will be assigne to DHLPOT. 

A 5% uncertainty for ETAHPOT is being assigned based upon engineering judgment since no 
information was obtained. 


T ORQHPOT=DHHPOT ■ WHPOT — 

SNRADOH 


U% ETAHPOT -° 05 
TORQHPOT =54740.53 


DHHPOT=ETAHPOT DHEDLHPOT 


TORQHPOT=ETAHPOT DHDLHPOT- WHPOT— — — 

SNRADOH 


TRQNOH := TORQHPOT 


Since this balance relation is 
being worked prior to working #64 
uncertainty for WHPOT is 
being estimated 


U% TORQHPOT : " J U% ETAHPOT 2 + U ° /o WHPOT 2 U /0 WHP0T U /o WL? 01 

U% TORQHPOT = 007 1 
U TORQHPOT : = TORQHPOT- U% TORQHPOT 

U TORQHPOT = 3.91- 10 3 

U TRQNOH := U TORQHPOT 

U TRQNOH =3 - 91 * 1C)3 
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Rewriting the balance relation and determining the uncertainty 

BAL73 : = TRQNOH- TRQPOH 
BAL73 =-0.02 

I 2 2 

U BAL73 : = a] u TRQNFH + U TRQPFH 
U BAL73= 1191 * 104 


Now write the uncertainty as a percentage of the RHS of the balance relation: 


U% BAL73 - 


U BAL73 
|TRQNOH| 


U% BAL73 =21.749*% 
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Balance Relation 74: PTLPFP=PTVL01 

* PTLPFP is calculated in SPUMP and utilizes the LPFP head coefficient map 

* PTVL01 is an iteration variable 

From the SPUMP module, the exit pressure and uncertainty are determined from 

PTHTNK =39.01 


PTLPFP := PTHTNK + PSILPFP- RHOVLOl - SNRADFL 2 

PTLPFP =229.477... . _. . . .. ..„ 

Note: This value is slightly different 

than the model calculated value of 
PTLPFP=229.5579 


PSILPFP := 0.0357 
PTVL01 := 229.5579 


Since PTHTNK is an input value and SNFL is an iteration variable, 

PSILPFP and RHOVL01 are the terms which contain uncertainties. 

Using an uncertainty of 3% for the head coefficient from the „ . ... pFp 

head coefficient vs flow coefficient map. Estimate based upon PSILPFP ~ 

information provided by Rocketdyne. 


PSILPFP- SNRADFL 
8 psiLPFP ' ~ RHOVLO j . SNRADFL 

U PTLPFP = J( e PSILPFP' U PSILPFP) 2+ ( 0 RHOVLO1 U RHOVLOl) 2 
U PTLPFP = 57 17 

Rewriting the balance relation and determining the uncertainty 
BAL74 =PTVL01 - PTLPFP 
BAL74 =0.081 


0 PTLPFP := ~ 1 


U BAL74 ~J( e PTLPFP' U PTLPFp) 2 
U BAL74 =5 - 717 

Now write the uncertainty as a percentage of the RHS of the balance relation 

no/„ U BAL74 

BAL74 |p TVL01 | 

U% BAL74 =2 - 491 ‘ % 
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Balance Relation 75: PTHPFP=PTVL03 
PTHPFP is calculated in SPUMP and utilizes the HPFP head coefficient map 
* P7VL03 is an iteration variable 

From the SPUMP module, the exit pressure and uncertainty are determined from 

PTVL02 := 196.507 


PTHPFP :=PTVL02 + PSIHPFP RHOVL03- SNRADFH 2 

PTHPFP = 5.8227 1 • 10 3 ... . 

Note: This value is slightly different 

than the model calculated value of 
PTHPFP=5821 .24 


PSIHPFP =0.1567331 
PTVL03 := 5821.243 


Since PTVL03 is an input value and SNFH are iteration variables, 

PSIHPFP and RHOVL03 are the terms which contain uncertainties. 

Using an uncertainty of 3.5% for the head coefficient from the 

head coefficient vs flow coefficient map. Estimate based upon U PSIHPFP ~ 0 035 PSIHPFP 
information provided by Rocketdyne. 

©HHOVL03 = PSIHPFP' SNRADFH 2 
© psiHPFP : = RHOVL03- SNRADFH 2 

U PTHPFP := J( e PSIHPFP U PSIHPFp) 2 ( e RHOVL03 u RHOVL03) 2 
U pxhPFP ~ 390.808 

Rewriting the balance relation and determining the uncertainty 
BAL75 : = PTVL03 - PTHPFP 
BAL75 =-1.465 

6 PTHPFP = * 1 

U BAL75 := J( d PTHPFP' U PTHPFp) 2 


U BAL75 = 390.808 


Uncertainty as a percentage of the right-hand side of the balance relation: 
U BAL75 


U% 


BAL75 - 


-■100 


PTVL03 
U% BAL75= 6713 

Now write the uncertainty as a percentage of the RHS of the balance relation 


U% 


U 


BAL75 


BAL75 


U% 


BAL75 


|PTVL03| 
= 6.713*% 
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Balance Relation 76: PTLPOP=PTVL19 

* PTLPOP is calculated in SPUMP and utilizes the LPOP head coefficient map 

* PTVL1 9 is an iteration variable 

From the SPUMP module, the exit pressure and uncertainty are determined from 

PSILPOP : = 0.0276 PTVL18 =74.07178 

PTLPOP :=PTVL18 + PSILPOP RHOVL19 SNRADOL 2 RHOVL19 :=00404PTVL19 :=396.0983 

PTLPOP = 395.776. Note; Jhjs vg|ue fe s|jght | y djfferent 

than the model calculated value of 
PTLPOP=396.0983 

Since PTVL18 and SNOL are iteration variables, 

PSILPOP and RHOVL1 9 are the terms which contain uncertainties. 

U RHOVL19 =0.002 RHOVL 19 

Using an uncertainty of 5.0% for the head coefficient from the -n^PQnpnp 

head coefficient vs flow coefficient map. Estimate based upon u PSILPOP ~ • 1L 
information provided by Rocketdyne. 

0 RHOVL 19 = PSILPOP- SNRADOL 2 
6 PSILPOP : = RHOVL 1 9- SNRADOL 2 

U PTLPOP : = j( 0 PSILPOP U PSILPOP) 2 + ( e RHOVL19' U RHOVL19) 2 
U PTLPOP = 16098 

Rewriting the balance relation and determining the uncertainty 

BAL76 :=PTVL1 9- PTLPOP 
BAL76 =0.322 


0 PTLPOP — 1 


U BAL76 -J( & PTLPOP' U PTLPOp) 2 
U BAL76 = 16098 

Now write the uncertainty as a percentage of the RHS of the balance relation 


U% BAL76 : 


U BAL76 
I PTVL19 1 


U% BAL76 = 4 064 
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Balance Relation 77: PTHPOP=PTVL21 

* PTHPOP is calculated in SPUMP and utilizes the HPOP head coefficient map 

* PTVL21 is an iteration variable 

From the SPUMP module, the exit pressure and uncertainty are determined from 

PTVL20 := 362.0607 PTVL21 :=4123.193 

PTHPOP =PTVL20 + PSIHPOPRHOVL21SNRADOH 2 PSIHPOP =0.0112 RHOVL21 =0.0403 

PTHPOP =4.1 1975-10 3 #M , , 

* Note: This value is slightly different 

than the model calculated value of 
PTHPOP=4123.194 

Since PTVL20 and SNOH are iteration variables, 

PSIHPOP and RHOVL21 are the terms which contain uncertainties. 

U RHQVL21 := 0-002- RHOVL21 

Using an uncertainty of 6% for the head coefficient from the 

head coefficient vs flow coefficient map. Estimate based upon u PSIHPOP ' ihpop 
information provided by Rocketdyne. 

©RHOVL21 - PSIHPOP SNRADOH 2 

8 psiHPOP ' " PHOVL2 1 • SNRADOH^ 

U PTHPOP := J( @ PSIHPOP' U PSIHPOp) 2 + (® RHOVL2 1 U RHOVL21 ) 2 
U PTHPOP = 225587 

Rewriting the balance relation and determining the uncertainty 

BAL77 :=PTVL21 - PTHPOP 
BAL77 =3.44 

® PTHPOP = - 1 

U BAL77 = J( e PTHPOP - U PTHPOp) 2 
U BAL77 = 225.587 


Now write the uncertainty as a percentage of the RHS of the balance relation" 


U% BAL77 :: 
U% BAL77 


. U BAL77 
|PTVL21 1 
= 5.471*% 
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Balance Relation 78: PTPRBP=PTPBSO 

* PTPRBP is the preburner pump discharge pressure and calculated in SPUMP 

* PTPBSO is the oxidizer splitter duct and is an iteration variable 

The preburner pump inlet pressure is calculated in the configuration file using the equation: 


PTPBPI :=PTVL21 - 


4.47-10' 
(RHOVL21) 


- WPRBP 2 


PTPBPI = 3.99356-10* 


Note, this value is slightly different than 
the value in the output of 3994.133 


WPRBP := 108.1089 
PTPBSO := 6932.563 


RHOPBSO = 0.0408 
U RHOPBSO ' 0.002RHOPBSO 


From the SPUMP module, the exit pressure and uncertainty are determined from 

PSIPRBP =0.00864 


PTPRBP := PTPBPI +- PSIPRBPRH0PBS0-SNRAD0H 2 


PTPRBP = 6.9283 1-10 3 


* Note: This value is slightly different 
than the model calculated value of 
PTPRBP=6932.563 


Since PTPBSO and SNOH are iteration variables, 

PSIPRBP and RHOPBSO are the terms which contain uncertainties. 

Using an uncertainty of 5% for the head coefficient from the head coefficient vs flow coefficient 
map. Estimate based upon engineering judgment, since no information was provided by 
Rocketdyne for the Prebumer Pump. 

U psiprbp = 0.05- PSIPRBP 


e RHOPBSO = PSIPRBP SNRADOH 2 
0 PSIPRBP = RHOPBSO SNRADOH 2 

U PTPRBP := J( 0 PSIPRBP' U PSIPRBP) 2 + ( 0 RHOPBSO U RHOPBSo) 2 
U ptprbp ~ 146.855 
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Rewriting the balance relation and determining the uncertainty 
BAL78 :=PTPBSO- PTPRBP 
BAL78 =4.249 

®PTPRBP = " 1 

U BAL78 : = J ( e PTPRBP’ U PTPRBp) 2 
U BAL78 = 146855 

Now write the uncertainty as a percentage of the RHS of the balance relation 

U BAL78 
BAL78 | PTpBSO | 

U% BAL78 =2118 '° /o 
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Balance Relation 79: XROPV=XROPOV 


BAL79=XROPO V - XROPV 


Balance Relation 80: PTHG2=PTOPRB 
BAL80=PTOPRB - PTHG2 


Balance Relation 81 : PTHG6=PTOSF 

BAL81=PTOSF- PTHG6 


Balance Relation 82: PTHG1=PTFPRB 

BAL82=PTFPRB - PTHG1 


Balance Relation 83: PTHG5=PTFSF 

BAL83=PTFSF - PTHG5 


Balance Relation 84: PTFINJ=PTMFI 
BAL84=PTMFI - PTFINJ 
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Balance Relation 85: WMCHB=WNOZL 


BAL85=WNOZL- WMCHB 


Balance Relation 86: WHG1=WFTBP 

BAL86=WFTBP - WHG1 


Balance Relation 87: WHG2=WOTBP 

BAL87=WOTBP - WHG2 


Balance Relation 88: PTMCHB=PCREQ 
BAL88=PCREQ- PTMCHM 


Balance Relation 89: EMRFCNTR=OFREQ 

BAL89=OFREQ - EMRFCNTR 


Balance Relation 90: QNETNOZX=0 
Balance Relation 91 : QNETNOZU=0 
Balance Relation 92: QNETCHML=0 


Balance Relation 93: QNETCHMU=0 



Summary of balance relation uncertainties 


U% BAL1 = 5 *° /o 

U ° /o BAL2 = 5 - 032 " 0/o 

U% BAL3 =5 *° /o 
U% BALA = 5,0/0 

770/0 BAL5 =4.876*% 

U% BAL6 = 13055 •% 

u %baL7 = 5 -786*% 
770/0 BAL8 =5786*% 
U% BAL9 =7071 * % 
U ° /o BAL10 =7071 *° /o 


U% BALI 1 = 7071 *° /o 

U% BAL12 = 7071 * % 

U% BAL13= 7 - 681 * % 
U ° /o BAL14 = 23613 ,% 

U% BAL15= 6 - 986 * % 

U% BAL16 = 8718 * % 

U% BAL17 =7071 * % 

U% BAL18= 7071 * % 

U% BAL19 =7071 * % 

U% BAL20= 7071 ‘ % 


U% BAL21 = 7071 * % 

170/0 BAL22 = 7071 *° /o 

U% BAL23 =6 * % 
U %BAL24 =4a769 * % 

U% BAL25 = 7-293 •% 

U% BAL26 =7 - 292 *° /o 

U ° /o BAL27 =6.604*% 
170/0 BAL28 = 7 -4 1 9 *% 
170/0 BAL29 = 7 282 *° /o 
U% BAL30 =7 - 282 *° /o 


U ° /o BAL31= 1 4- 1 42*% 
U% BAL32 = 14141 *° /o 

U ° /o BAL33 = 14142 ,0/o 
UO/o BAL 34= 2071 ‘ % 
U% BAL35 = 10 *° /o 
U% BAL36 =7062 * % 

U% BAL37 = 7 -06 2 *% 
U ° /o BAL38= 5 - 001 * 0/o 
170/0 BAL39 = 5(3<3 1 *% 
U ° /o BAL40= 5071 * O/o 


U% BAL41 = 509 *° /o 
770/0 BAL42 =4 - 93 9 *% 

U ° /o BAL43 = 4 - 941 *° /o 

770/0 BAL44 = 3764 *% 

770/0 BAL45 = 3.8 7 *% 

UO/o BA1A6 =7143 *° /o 

770/0 BAL47= 5001 *° /o 
UO/o BAL48= 361 - % 

770/0 BAU9= 3 - 649 * 0/o 

■■■■= *% 
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Abstract 

Providing accurate performance predictions for 
a liquid rocket engine is essential to providing a safe 
and effective space launch system, since test 
programs cannot, in a cost efficient manner, provide 
all the necessary information. This paper discusses a 
new methodology which includes knowledge of the 
uncertainties associated with the test program as well 
as the uncertainties associated with the analytical 
model and the uncertainty associated with the 
solution of the model. This methodology will be 
used during the analysis of test data through a test 
data reconciliation strategy to incorporate the test 
data and its uncertainty within the model solution to 
provide the most plausible region of engine 
performance conditions. 

This work is based around e nhancing the 
analysis tools for the continued evolution of the 
Space Shuttle Main Engine (SSME) and its on-going 
development at the Technology Test Bed (TTB) 
facility at Marshall Space Flight Center. Assessment 
of the experimental uncertainties associated with 
TTB SSME test data and assessment of the 
uncertainties associated with the analytical model 
are also discussed. 
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Introduction 

The conflicting demands placed upon current 
and future liquid rocket engine systems, maximum 
performance at minimum cost, require the 
development of new analysis tools, tools which 
provide more information upon which to base 
decisions. The uncertainties in the test results and 
the uncertainties in the analytical model predictions 
have often been neglected. Hence, a research 
program has been in progress, led by NASA 
Marshall Space Flight Center (MSFC), to develop 
analysis tools for use with the Space Shuttle Main 
Engine (SSME) which considers these uncertainties. 
The intent of these tools is to provide information 
which identifies the region of most plausible engine 
performance conditions. 

This paper describes the combined effort of 
researchers from NASA/MSFC, Christian Brothers 
University 1 , and the University of Alabama in 
Huntsville 2 to develop these tools. This paper 
presents the background and overall objective of the 
effort, aspects of the experimental uncertainty 
analysis, aspects of the data integration and 
modeling work, and a discussion of future work. 

Background and Objective 

Analytical performance models simulate the 
internal flow environment of a rocket engine by 
modeling subsystem components and integrating 
these components into an engine system. These 
performance models are essential in the support of 
SSME hardware design, evaluating hardware 
performance, predicting engine operating for off- 
nominal hardware conditions, and anomaly 
investigations. The primary purpose of these 
performance models is to accurately predict the 
performance characteristics of the individual engine 
components. These individual components are 
modeled using analytical relationships which 
describe tuibopump characteristics, pressure losses. 
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and heat tr ansf er effects. These analytical 
relationships are adjusted to correlate to test data. 
Thus, the quality and validity of the performance 
predictions generated rely primarily on the process 
in which test data is integrated with the analytical 
relationships. If this process is flawed, the hardware 
characteristics and performance predictions 
generated from these models will be inaccurate, 
unreliable, and could steer the rocket engine 
development program on a costly and undesirable 
course. 

Accurately predicting the performance 
characteristics of the individual engine components 
is extremely difficult because of three factors. First, 
limited instrumentation exist on most rocket engines. 
There are approximately 700 parameters within the 
current SSME performance model which represent 
all of the predicted pressures, temperatures and 
flowrates within the engine. The instrumentation on 
the Technology Test Bed SSME provides data for 
only 12% of these internal flow parameters. 
Secondly, rocket engine testing is extremely costly. 
This limits the amount of data that can be attained 
for a specific rocket engine design. Third, most 
rocket engine performance analysis models integrate 
test information wi thin the physical relationships by 
assuming that the measured test data is absolutely 
accurate. This, so-called, "pristine data" assumption 
is flawed since all measurement systems have 
inherent inaccuracies. In addition, rocket engine 
models consist of analytical relationships based upon 
simplified assumptions. Individual engine 
components are modeled assuming one-dimensional, 
uniform, flow-averaged conditions. Local 
measurements in this complicated flow network are 
incapable of characterizing these flow-averaged 
conditions with accuracy. Computational 
predictions are forced to agree with the data at 
instrumented locations, often at the expense of 
balancing fundamental physical relationships (mass 
conservation, momentum, and energy conservation). 
Consequently, this leads to inaccurate hardware 
characterization and engine system performance 
predictions. The first two factors are inherent within 
all rocket engine technology. The third factor can be 
significantly improved by developing new methods 
to integrate uncertain test information within 
simplifi ed analytical relationships. 

A new generalized performance analysis model, 
termed the Performance Reconciliation Model 
(FRM), is currently being developed which contains 
a test data integration strategy that eliminates the 
pristine data assumption. The solution procedure 
incorporates a system level nonlinear mathematical 


optimization strategy. The optimization strategy 
requires uncertainty estimates for both test 
measurements and the balance of the fundamental 
physical relationships. The optimization procedure 
selects values for engine variables (flow rates, 
pressures, temperatures, speeds, etc.) within defined 
measurement uncertainty limits, wliile satisfying the 
fundamental physical relationships within defined 
balance uncertainty limits. The objective of the 
optimization procedure is to minimize adjustments 
in the engine variables from the initial test values. 
The PRM generates a set of performance predictions 
by reconciling test data with fundamental physical 
relationships within specified uncertainty limits. In 
effect, PRM predictions provide the most plausible 
set of engine performance conditions, recognizing 
the inherent uncertainties in both fundamental 
physical relations and test measurements. Figure 1 
schematically shows the fundamental aspects of this 
new computational strategy. Test data obtained from 
the SSME Technology Test-Bed facility at MSFC is 
being used to support this model development. The 
goal of this strategy is to get the test information and 
fundamental physical relationships to reinforce each 
other and thus, more accurately predict the 
performance characteristics of the individual engine 
components. 


Engine Test 
* Uncertain 
Measurements 



Engine Model 

* Modeling Uncertainties 

* Uncertainties from 
Previous Information 



Performance 
Reconcilation Model 


I 

“Most Plausible" Engine 
Performance Conditions 


Figure 1 Schematic of Research Effort 


Uncertainty Analysis in TTB Testing 
The focus in this task is to identify the uncertainty 
that should be associated with measured variables such 
as temperature or pressure or with a determined result 
such as flowrate that is calculated using a number of 
measured variables. 
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Desired Variable 


l r ^ g- Environment 

(InstalMai&ccnceptual 
biases; unsteaciness) 


Sensor 

(Calibration biases) 


Etta Acqusiticn System, DAS 

^ 1 (Calibration biases) 

Masured Value of 
triable 

Figure 2 Schematic of TTB Measurement Process from an 
Error Sources Viewpoint 

Figure 2 shows a schematic of the viewpoint used 
in identifying error sources that contribute to the overall 
uncertainty. The desired variable is taken to be the one 
with which a model output will be compared - a cross- 
section averaged temperature, for example, that would 
usually be referred to as "the temperature" of the flow at 
a particular location in the engine. If the sensor 
responds to temperature at a point, then an installation 
or conceptual bias exists due to the sensor not actually 
responding to the desired variable (the average 
temperature). This is an elemental source that must be 
included, and it is potentially one of the dominant 
elemental sources in temperature and pressure data in 
TTB testing The traditional "measurement 
uncer tainty " sources are shown as biases in the sensor 
calibration and biases in the calibration of the data 
acquisition system (DAS). Additionally, the effect of 
unsteadiness in, and due to, the operating environment 
must be considered since the sensor calibrations and 
DAS pre-test calibration checks are not done with the 
engine operating, and the unsteadiness certainly can 
have an effect on the final system output - the measured 
value of the variable. 

Choice of the appropriate precision uncertainty 
estimate (precision limit) to use with TTB data needs to 
be carefully done. A precision limit determined using a 
standard deviation from a time slice during one test 
gives information about the steadiness of the "steady 
state" at that operating condition during that particular 
test, but includes no effects of the test-to-test variation of 
the variable at that operating condition. Computing a 
standar d deviation of a variable or result from multiple 
tests, all of which were at the same operating condition. 


gives the appropriate precision estimate for use in 
discussing the uncertainty in a measured TTB variable. 
It is also the appropriate precision limit to consider 
when comparing the results from one test to results from 
another test in an effort to determine if a change in 
component, for instance, had any discernible effect on 
the value of the result 

An investigation to determine the experimental 
uncertainties associated with test measurements from 
the SSME Engine 3001 installed in the Technology 
Test Bed facility was conducted. This investigation 
consisted of reviewing previous analyses, discussions 
with NASA personnel, review of other technical 
literature, and new analyses. Since the 
thermodynamic performance analysis of the SSME 
was the motivation behind this contractual effort, the 
pressure, differential pressure, temperature, and 
mass flow rate measurements were the focus of the 
initial investigation. 

Analysis of Previous TTB Test Data 

To achieve a set of data to review which could 
be defined as from the “same ' hardware, tests 
TTB039 through TTB 051 were chosen. These tests 
were conducted with Engine 3001 with the large 
throat combustion chamber and a consistent set of 
other hardware. 

The mass flow rate uncertainties for the venturi 
flowmeters were determined using the methodology 
in references [3], [4], [5] and [6]. The systematic 
uncertainty estimates (bias limits) used in the 
uncertainty propagation equations were based upon 
the information gathered from the available 
documentation, discussions with TTB personnel, and 
engineering experience. 

Potentially significant bias uncertainties yet to 
be evaluated in the TTB measurements are the 
conceptual bias uncertainties. The conceptual bias 
uncertainties in the temperature and pressure 
measurements are particularly important in this 
effort because of the interest in comparing the 
experimental results with the analytical predictions. 
In many of the SSME measurements the flowfield is 
highly complex due to the sharp turns and bends, 
valves, pump and turbine inlets and discharges, and 
other complicating factors. These factors accentuate 
the difference between the physical quantity at the 
sensor and the quantity for which the measurement 
is desired, typically an average value at a cross- 
section. These assessments will require extensive 
review of the measurement the sensor and its 
installation, the thermodynamic and fluid dynamic 
flowfield, and their interaction. This detailed 
analysis of the measurements is currently underway. 
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The precision limits for the mass flowrate 
uncertainties are dependent upon the question being 
asked, or rather, what is purpose for the information. 
Precision limits can be calculated in many different 
ways, but the interpretation of the precision limit and 
the use of it depends upon data used to calculate it. 
The variables which must be considered for the 
precision limit calculation include, but are not 
limited to: 

• engine number 

• specific engine component configuration 

• engine test(s) 

• power level 

• test profile 

• specific engine adjustments 

• time slice within the test 

• data sample rate (data points used for 
s tandar d deviation calculation) 

The precision limits for the flowrate 
uncertainties were based upon review of the flowrate 
data for the chosen time slice. Precision limits were 
estimated in two primary ways. First, one was based 
upon the full sample-rate (25 or 50 samples/sec) data 
within each test The second precision limit was 
estimated based upon averaging the full sample-rate 
data over a given time slice, for a given test at a 
chosen power level to provide a single data point for 
that test condition and using similar points from 
other tests to form a sample from which a precision 
limit estimate could be calculated. 

The uncertainty in the thermophysical property 
data must be included in the uncertainty analysis 
since there was some experimental uncertainty in the 
original experiments upon which the property tables 
were developed. If the property data is represented 
by a curve-fit of the experimental data an uncertainty 
associated with the data and with the curve-fit must 
be included. The venturi flowrate data reduction 
utilises the thermophysical property routine GASP 7 . 
The GASP documentation shows foe uncertainty in 
density for liquid hydrogen and liquid oxygen to be 
within approximately 0.25% outside of foe critical 
region and within approximately 2.0% near foe 
critical region. The GASP program is based upon 
National Bureau of Standards data for foe 
thermodynamic properties of hydrogen and oxygen. 

Calibration of Mass Flowmeters, Determination 
of Discharge Coefficient, C D 

The accurate determination of foe mass flow 
rate for cryogenic rocket propellants is difficult 
because of the special problems presented. These 
problems arise because of foe cryogenic 


temperatures, the property data uncertainty, the 
differences in calibration fluids (water or air versus 
cryogens), and calibration limitations due to safety 
considerations Each of these can introduce an 
uncertainty in foe mass flow rate determination. 

The accurate determination of a mass flow rate 
depends upon foe ability to trace the output from a 
given device back to a standard certified by the 
National Institute of Standards and Technology 
(NIST) or some other respected standard. The basic 
problem with cryogen flowmeters is that there are 
limited facilities which can produce an accurate 
standard using foe actual cryogenic fluid. The 
facilities which do exist are limited to foe calibration 
of liquid nitrogen flows, or have limited capacity. 

The principle of dynamic similarity, matching 
foe Reynolds number of foe test fluid with foe 
Reynolds number of foe calibration fluid, has been 
relied upon, with corrections for the dimensional 
changes due to thermal contraction and other 
miscellaneous effects being made. A literature 
survey indicates that very little work has taken place 
to assess the accuracy of this procedure, primarily 
due to foe cost and complexity involved with 
developing the necessary experiments with liquid 
oxygen and liquid hydrogen. The systematic 
uncertainty due to not calibrating foe venturi 
flowmeters in foe appropriate cryogenic fluid has 
been estimated to be as much as 1%. 9>0 Other 
assumptions and approximations are made to 
achieve the venturi mass flowrate data reduction 
equations. Two of the fundamental assumptions are 
that the flow through foe venturi is adiabatic and 
one-dimensional, which is not true during engine 
operation. Another significant factor is that none of 
foe venturis are installed with the recommended 
length of upstream straight duct. An attempt was 
made 8 to account for this problem by calibrating foe 
venturis with the “as-installed” ducts connected, 
however they could not reproduce foe other flow 
characteristics, such as flow swirl, turbulence, 
oscillations, etc. 

The venturi calibration report 8 develops a 
polynomial curve-fit to obtain the value of foe 
discharge coefficients where foe curve-fit is based 
upon the Reynolds number of foe flow through foe 
venturi. For all of the venturi flowmeters calibrated 
with water dynamic similarity could not be obtained, 
that is, foe flowmeters could not be calibrated at foe 
Reynolds number expected during an engine test. In 
fact, most of the calibrations were performed at a 
fraction of foe operating Reynolds number, from 
1.7% to 10%, and foe curve-fits were extrapolated 
out to foe operating condition. 
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The extrapolation of the calibration curve-fits 
and the use of a different fluid for the calibration 
introduce an uncertainty which must be associated 
with the venturi discharge coefficient. In the 
absence of performing detailed experiments to assess 
these uncertainties, estimates must be made based 
upon the existing information. In this study, 
systematic uncertainties of 2-4% were estimated for 
the discharge coefficients. 

Uncertainty in Mass Flowrate Determinations 

A computer program was developed to calculate 
the flowrates for the venturi flowmeters and to 
perform the uncertainty calculations. This program 
accesses the raw, full sample-rate data for all of the 
engine measurements. The program uses the 
measured pressure, differential pressure, and 
temperature for each venturi, combined with the 
specific dimensional data and other constants (such 
as discharge coefficient), averages the data for the 
chosen time slice and calculates the mass flowrate. 
The bias limi ts and precision limits for the 
uncertainty sources for each venturi are stored in the 
program as constants and propagated through the 
data reduction equation using the uncertainty 
propagation equations. The only correlated bias 
uncertainty included in the uncertainty propagation 
equation is the uncertainty in the venturi 
dimensions. All other possible bias correlations 
were considered negligible, particularly with respect 
to the magnitude of the discharge coefficient bias 
uncertainty. The partial derivatives are determined 
numerically using a finite difference technique, often 
referred to as a jitter routine. 

The mass flow rate uncertainties were 
determined to be in the 4% to 8% range, with the 
two dominant uncertainty contributions being the 
systematic uncertainty from the discharge coefficient 
and the precision uncertainty. Further review of the 
test data indicated that most of the precision 
uncertainty was due to the variation in the engine 
balance from test-to-test for a given power level. A 
factor which complicates the analysis of the SSME is 
that the engine will balance at different operating 
points for a given power level. Hence, the test-to-test 
uncertainty estimate contains this information and is 
the appropriate interval to use within the PRM. 

Uncertainty Analysis in SSME Modeling 

When comparing output of a model with 
experimental data, the uncertainties that should be 
associated with the model predictions must be 
considered for proper conclusions to be drawn. In the 
past, most of the work reported in this area has simply 


considered the sensitivity of the model output to 
uncertainties in the input data. This obviously does not 
include any uncertainties in the model itself and thus is 
not a satisfactory approach. In this research effort, we 
have divided the sources that cause uncertainty in the 
model output into three categories: ( 1 ) uncertainties due 
to assumptions and approximations in the model. (2) 
uncertainties due to the incorporation of previous 
experimental data into the model, and (3) uncertainties 
due to the numerical solution algorithm. Consideration 
of the third category is considered in the next section. 

The first category, uncertainties due to assumptions 
and approximations in the model, does not include the 
installation and/or conceptual bias source shown in 
Figure 2 and discussed above since that uncertainty is 
associated with the measured value. Consider the 
temperature at a particular position in the flow. The 
uncertainty associated with the measured value of the 
temperature includes the effect of making a point 
measurement but desiring a cross-sectional averaged 
value. The inability of the model to calculate a correct 
average temperature at a particular location because the 
one-dimensional flow approximation has been made 
results in an uncertainty in the predicted temperature. 
(Stated another way, if the model predicts the correct 
average temperature at a particular location, then the 
one-dimensional flow approximation has caused no 
uncertainty in the model output.) 

The uncertainties due to the incorporation of 
previous experimental data in the model arise when 
material property data is used, when valve resistance 
characteristics are used, when pump maps are used, etc. 
These are all instances in which previous experimental 
data has been used by replacing the data with curvefits. 
The original data contained uncertainties, but the 
curvefit equations used in the predictive models have 
been treated as the "truth" in most previous 
considerations of uncertainty in model outputs. Adding 
further complication, there is no accepted way of 
estimating the influence of systematic uncertainties on 
the uncertainty that should be associated with a 
regression. This aspect has been investigated as a part 
of this program and is discussed later. 

Modeling Assumptions and Approximations 

When a model of a physical sy stem is developed, 
assumptions and approximations about the system 
are made to simplify the system to one which 
mathematical expressions can describe. By making 
these simplifications an error is introduced and the 
model cannot exactly describe the physical system. 

Some of the primary assumptions and 
approximations made within the SSME model 
include: 
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• 1-dimensional 

• fully developed 

• steady-state 

• adiabatic 

• ideal gas 

If the uncertainty to associate with a particular 
assumption or approximation can be estimated, then 
sometimes the model should be improved to include 
this estimate instead of trying to estimate the 
uncertainty. For example, it was determined that a 
turbine exit temperature was being predicted using 
an ideal gas, constant specific heat approximation, 
which for the specific temperature range of interest 
was a poor approximation. Instead of trying to 
estimate an uncertainty to associate with that 
approximation, the model was altered to include a 
better thermodynamic description of the process. 

If the models are developed based upon the 
soundest assumptions and approximations available, 
assessing the uncertainty due to these assumptions 
and approximations would require an effort beyond 
the scope of the current research program 

Uncertainties From Using Previous Experimental 
Information 

In all of the component modules, information 
from previous testing is used. For example, the 
model of the liquid oxygen flow through a given duct 
or through a given valve is based upon its component 
testing, which provides an equation for the resistance 
through the duct as function of the flowrate. This 
test information is often reduced to the form of a line 
or curve. The values of the polynomial constants in 
the thermodynamic property routines are also 
examples of using previously obtained uncertain test 
information in a model. 

The methodology to assess the uncertainty in the 
coefficients of a linear regression were developed as 
part of this effort The uncertainty analysis 
methodology previously discussed was applied to the 
expressions for the regression coefficients to develop 
the technique. The details of this methodology were 
presented at the 1995 AIAA Aerospace Sciences 
conference". The work presented in that paper 
demonstrates that this technique provides a method 
to determine the uncertainty in linear regression 
coefficients, and incorporates the effects of correlated 
systematic uncertainties. 

Recently, the methodology has been 
extended to provide the appropriate uncertainty 
interval for the predicted value using the regression 
model equation. Work is continuing to develop the 
specific techniques needed to apply this methodology 


to the types of models and previous experimental 
information found in the SSME program. 

Performance Analysis 

Data integration strategies are defined by 
formulation, including both physical and 
mathematical components, and computational 
solution procedure. A reliable strategy requires a 
computationally tractable and physically consistent 
performance model as a fundamental starting point. 
In mathematical terms, the computational steady- 
state performance analysis problem can be stated as 
follows: 

Problem 1 Computational Performance Analysis 
Problem 

determine X 

such that F ( X ; Y ) < Tol 

where 

F - the column vector [Fj , F 2 , . . . , F n ] T 
composed of n independent balance relations that 
govern modeled engine performance 

X - the column vector [ X| . x ; , . . . , x n ] T 
composed of n independent physical characteristics 
that define modeled engine performance 

Y - the column vector [ y[ . y 2 y m ] T 

composed of m performance parameters, assumed 
constant, for the modeled engine 

Tol - the column vector [ toll . tol 2 , . . . , tol n ] T 
composed of the n solution tolerances associated 
with the n independent balance relations governing 
modeled engine performance 

Acceptable solutions to Problem 1 are located 
within the intersection of tolerance bands about the 
individual balance solution contours. The method 
for selecting individual tolerances is an important 
consideration which can significantly affect the 
solution obtained. Despite the critical nature of 
tolerance selection, heuristic guidelines are generally 
used to specify individual relation tolerances. A 
more systematic procedure for specifying the 
acceptability threshold is to define tolerance limits 
in terms of balance relation uncertainties. Balance 
relation ( F,(X;Y)=0 ) uncertainty arises from a 
variety of sources including model 
assumptions/approximations, scarce/ inaccurate 


6 



thermodynamic property data in ranges of interest, 
scatter/extrapoiation of experimental data used to 
estimate parameters in semi-empirical relations, and 
secondary effects such as small scale leakage often 
ignored in large scale performance models. 

Techniques for estimating experimental 
uncertainties and their propagation are well 
developed, as discussed previously. Although not 
uniformly applicable to tolerance estimation, these 
techniques can be expanded to provide a 
fundamental approach to balance uncertainty 
estimation. Once defined, the balance uncertainty 
vector U=[ Ui , U 2 , . . . , U n ] T can be used in place 
of the tolerance vector Tol in Problem 1, yielding a 
modified computational performance analysis 
problem in which uncertainty band intersections 
define the solution feasible region as depicted 
conceptually in Figure 3. This type approach is 
logically consistent and provides consistent 
utilization of available information in performance 
modeling. 



Figure 3 Feasible Region, Physical Relations 


Y„ - a subset of Y composed of k 
hardware parameters that are treated as independent 
data reduction solution variables 

Yf - the complementary subset (Yt=Y- 
Y v ) of m-k hardware parameters that remain at 
assumed constant values 


Data Reduction 

The term "data reduction" indicates a specific 
procedure for integrating experimental test data 
within performance prediction models. 
Traditionally, data reduction has referred to data 
integration based on a single test firing of a specified 
engine. This does not, however, preclude the use of 
data reduction as a method for ensemble average 
data integration. 

The data reduction procedure is a strategy for 
modifying engine hardware performance 
characteristics in order to enforce model agreement 
with a specified set of experimental observations. To 
facilitate mathematical description of the data 
reduction process, it is useful to define the following 
set partitions: 

x = [ X, , x 2 Xk I XfcH xJ T = [ X f I X T ] T 
Y = [ yi , y 2 yk I yw,..,ym] T = [Y V |Y,] T 
where 

Xf - a subset of X composed of the k physical 
characteristics that are fixed by test data 

X y - the complementary subset (X v =X-X f )of n-k 
physical characteristics that remain as independent 
data reduction solution variables 


Using the above set partition definitions, the data 
reduction analysis problem can be stated as follows: 

Problem 2 Data Reduction Analysis Problem 

fix k physical characteristics X f at test data 
indicated values 

select k hardware parameters Y v to treat as 
variable 

determine X v and Y v 

such that F(X„,Y V ; X f ,Y f ) < Tol 

Because there is some freedom in the selection 
of the elements of Y v , the solution of the data 
reduction analysis problem is generally not unique^ 
This troublesome characteristics results from the fact 
that the system F(X,Y)=0 is underprescribed. It 
places a significant burden on the model user in 
selecting the variable hardware components of Y v . 
A logical, although not a priori correct, procedure is 
to select as elements of Y v those hardware 
characteristics which most significantly impact 
corresponding fixed elements of X f . In practice, the 
variable hardware characteristics are selected based 
on experience and related observation. This ripe of 
subjective selection provides unfortunate 
opportunities for systematic prediction bias. 
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A conceptual view of consistent data reduction 
analysis solutions, assuming absolutely accurate 
(pristine) test data, is presented in Figure 4. It is 
obvious that the set of acceptable physical solutions 
becomes more restricted as the dimension of X r 
increases (i.e., the number of test fixed physical 
parameters increases). 



Figure 4 Data Reduction, Consistent Pristine Data Solutions 

In addition, a significant problem arises when 
experimental observation falls outside the physically 
consistent solution region as displayed in Figures 5. 
It is an unfortunate consequence of experimental 
uncertainty that consistent data reduction becomes 
more difficult as the number of physical 
measurements increases. In the limi t, as the number 
of test fixed physical characteristics approaches the 
total number of prediction variables, a physically 
consistent performance solution is unobtainable due 
to test data measurement uncertainty. Traditional 
data reduction methods do not address this 
fundamental problem. A logical strategy is needed 
to resolve the integration anomalies depicted in 
Figure 5. 

Data Reconciliation 

Three approaches can be taken to resolve the 
data integration problems described in the previous 
section. The first, and perhaps simplest, technique is 
to retain the pristine data assumption and increase 
the physical relation tolerance bands to include test 
values. Uncertainty bandwidths would both guide 
and limit the allowable tolerance relaxations. A 
second approach is to consider measurement 
uncertainty as well as balance uncertainty in data 
integration. The intersection of balance and 
measurement uncertainty bands would then define 
all feasible solutions. The pristine data assumption 
can be retained if test measurements fall within the 
feasible region. In addition, the second approach 
provides a strategy for determining feasible 



Figure 5 Data Reduction, No Consistent Pristine Data Solution 

reconciliation solutions while relaxing the pristine 
data assumption . The third, and perhaps least 
desirable, approach is to simply ignore a portion of 
the flow physics in order to retain the pristine data 
assumption. This method is often used in current 
data reduction strategies when data restrictions 
preclude a consistent physical solution. This 
approach fails to recognize inherent measurement 
uncertainty and reduces physical consistency. This 
type of rigid adherence to the pristine data 
assumption is difficult to defend as a prediction 
enhancement strategy. 

The fundamental question and problem related 
to test data integration is simply stated. What is the 
best estimate of engine performance in the presence 
of uncertain, non-physical test data? The available 
flow physics and test data are insufficient to provide 
a rigorous answer to this question. Therefore, the 
quality of performance prediction depends in large 
measure on effective problem formulation/solution. 
Since best estimates are sought for an 
underprescribed engine system, it is logical to 
assume that an effective data integration model will 
take the form of a mathematical programming 
problem. Several such formulations based on 
squared residual sums have been developed during 
the interim phase of this study. These formulations, 
termed data reconciliation models, are described 
below. 

Reconciliation Formulation 1 

Assumes pristine data and physical relations 
uncertainty band satisfaction. 


n 

minimize f( X, Y ) = [F, (X v , Y v ; X f , Y f )] 

i=l 

by selection of X v , Y v 
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subject to 


|F(X t ,Y v ; Xf , Yf ) | < U 

Xf = test value 

Formulation 1 is the mathematical representation of 
data integration approach 1. No feasible solution 
exists if Xf does not fall within the physical relation 
uncertainty band intersection. 


In addition, it eliminates constraint complexity and 
utilizes both physical relation and test data 
uncertainty weighing to identify the "best' 1 operating 
state prediction. 

Reconciliation Formulation 4 

Assume uncertainty weighing of test data and 
physical relations. 


Reconciliation Formulation 2 

Assumes minimum deviation from test data 
and physical relations uncertainty band satisfaction. 

minimize f(X,Y)= 

xj sX r 

by selection of X , Y v 

subject to | F(X , Y* ; Yf ) | < U 

Formulation 2 is one example of a mathematical 
model reflecting data integration approach 2. It 
seeks the physically consistent performance state 
closest to test data The feasible region is the 
intersection of physical relation and test data 
uncertainty bands. In this formulation, the 
reconciliation solution is the point within the feasible 
region that deviates the least from test data. 

Reconciliation Formulation 3 

Assumes pristine data and minimum violation of 
physical relations. 


minimize f ( X , Y ) 



F,(X v ,Y v ;X f ,Y f ) 

U, 


-|2 


by selection of X», Y v 


subject to Xr = test value 

Formulation 3 presents a mathematical model 
related to data integration approach 2. It seeks the 
performance state with the highest level of physical 
consis ten cy that can be attained while enforcing 
absolute test data agreement. The reconciliation 
solution selected by formulation 3 is the operating 
state closest to the feasible region that agrees with 
test data. 

It is evident that reconciliation formulations 2 
and 3 return data integration solutions that differ 
only in the perceived quality of test data relative to 
physical consistency. Formulation 4 described below 
includes both formulations 2 and 3 as specific cases. 


minimize 


1=1 


f(X,Y) = 

F,(X,Y v ;Y f ) n2 

u, 


+ I 

X; £\ f 


(x,-x ute J 

U; 


by selection of X , Y v 

The effect of uncertainty weighing in the above 
squared residuals relation is obvious. Reconciliation 
formulation 4, as displayed in Figure 6. provides a 
logical and consistent platform for combining 
information from the historical data base (Y and U) 
with single test experimental data ( X f ) and physical 
model relations (F). It admits plausible "solutions" 
that approach physical and test consistency in an 
"optimum" fashion. In addition, it can be used to 
relax the requirement that the number of adjustable 
hardware performance parameters (in Y v ) match the 
number of available test measurements (in X f ) in the 
data integration process. 

The L 2 or squared residuals formulation 
presented above is but one of many potential 
extremum formulations that can be used for effective 
data integration. The formulation of choice will be 
the one that admits the most reliable and efficient 
solution. 
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Summary 

An uncertainty analysis of the SSME test 
program at the Technology Test Bed facility has 
shown that non-negiigible uncertainties exist in the 
determination of the mass flowrates, ranging from 
2% to up to 10%, and even larger in some cases. 
The systematic uncertainty is dominated by the bias 
limit in the experimentally determined venturi 
discharge coefficient. The precision uncertainty is 
dominated by the test-to-test variation in the power 
balance within the engine, causing a variation in the 
actual flowrate in addition to the inherent flow 
randomness. 

The uncertainty intervals provided in this report 
can be used for comparison to performance 
prediction models; however the engine configuration 
being modeled must be representative of the engine 
configuration for which the precision uncertainty 
was determined. 

A primary source of uncertainty in the 
analytical modeling of the SSME is from the use of 
previous experimental information. This 
information is usually utilized in the form of 
regressions or curvefits. Initial research efforts to 
determine a methodology to properly account for the 
uncertainty in the regression coefficients has shown 
to be promising, and work is co ntinuing to develop 
the methodology. 

Mathematical prog rammin g formulations 
provide a natural analytical basis for test data 
integration in performance models. Plausible 
formulations have been developed and implemented. 
Computational experience with these 
implementations is currently being extended. 
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Appendix 7 

Miscellaneous Background Information 

The information provided in this appendix is included to assist in 
reviewing and updating the uncertainty estimates. They are also included as 
background, to used to resolve currently open issues. 

The pump performance maps contained here were obtained from 
Rocketdyne, the proper correlation between these maps and the performance 
maps in ROCETS is currently not known. When the proper relationship is 
determined the transfer of the standard deviation bands can also be 
transferred. 

Portions of three reports from water testing of pump inlets and exits are 
provided. This information can be used to estimate the conceptual bias 
uncertainty for the pressure measurements located near the pump inlets and 
exits. This information was obtained too late in the performance period 
include in the report. 
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- Static Pressure 

- Total Pressure 

- Static Pressure Differentials 
* Total Pressure Differentials 

- Flow Angles 


"ThBlocations of these probes is indicated in figure A5. There was a number of specialized 
probes implemented in the model; these probes are described below. 

Total Pressure Kiel Probes 7^. 

-r , J&. 

^ .!? h SSUr9 - ki f ' ° r l beS ar9 USeti at ,hree >°“>i°"soWjobes are located at the ' 

test at the"imi>el’l.T o' 6 lea0l " a ® d 9 e ^- ssto a^»«rguide vanes, and for the dynamic 
est at the impeller discharge (See (figur^s_C63c7 and C8, respectively) The total 

ht° beS h at the m ° del inl6t Wer ® r ° tated fr ° m ° de 9 rees to +25.7. +52.4 and 
-25 7 to obtain the maximum number of points allowed by the model geometry of the flow 
profile exiting the SSME fuel supply duct. 


Wedge Probes 

Wedge probed are direction-sensingnrdBes; as the^me implies these'brobes have a 
we ge ortriangular shaped hea^fhere is a pressure port on each^ide of the wedge 

rotate^ho I" 9 Sd9e and l?^ reSSUr9 P^Tthe lea Oing edgdhe probe can be 
h S! 3n aX ' S per B^ icular to the^vedge allowing assignment with the flov 
wh^ he pressure differ^al across theZrts on the sides^e probe is zero the pSe 

znr e y^ na me ■*, >» P r* s, 

beS v ^ 1 ° Cated slightl^brward of the inlet aSde vane leading edgesbetween 
thJ ? ? , ne p Pf' x,ons and were^emoved when not ih use. The wedge ordbes provided 
the inlet volute^iischarge^ngl/ o verify inlet guide/vane angle of attadd 

Cobra Proves 



Uke the wepfde probes the Cobras were usdd to measure flow angle add total pressure' 

IOCateCl 90 de j^ apart With three Prao^wrow the span ot 
h J 9: dlSChar9e annulus ,n close Prp^imity to the impeller inlePptene during lhe sta fJc 

indued test. Although wedge pro^ are a more accurat^w angle XuS, 
in^rument they require rotation/Sapability and could ncrfbe accommodated this 
JT' a T 13 m ® asured us jP9 th ® curves obtained frp^ probe calibratipdcalibration 
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ure C6. Instrumented Inlet Ring 
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K nees can be observed from reviewing. the two test procepufes for static and 
First it witfoe noted that ipsfructions to ciearytfte inlet screens are 
lough, this was common practice during the static test it became cnjtea! 
as testing was performep / debris would build/rfp in screens and evemually 
ditior> / were inlet pressure would be lost duphg a data point (the Fluia Wafer 
irep aprox. 2 min, to scan alt channels) resulting in an unacceptable 
condition. /Also different from the static procedure is the 
aforementioned valve callouts /affected by the/dew flowpath required /to maintain the 
highest possible inlet pressures. The addition^ an oil cooled bearino/box to the system 
resulting/m the addition of/nstructions to ©Keck oil levels as a maintenance function as 
well as a set of instructions to ensure^oil tank heat exchanger was operating. The 
attention given to the bleed procedure has also been enhanced since the hook-up to the 
transducer panel was not as straight forward as during the static test. 


Results and Discussion 

Figure B-t shows the total and wall static pressure profiles versus circumferential position 
measured at the inlet flange. The instrumentation ring was clocked around the 
circumference of the flange to map the incoming flowfield. Total pressure measurements 
were taken at three different radial locations* each of which displayed similar trends, 
indicating peak pressures at 270 degrees and minimum pressures at 90 degrees 
circumferentially. The turning bend in the SSME duct upstream of the inlet flange, 
producing the inlet profile shown, results in higher flowrate and total pressure toward the 
forward side of the inlet duct and lower flowrate and pressure on the rear side near the 
inner wall of the bend. Radially, the data shows the lowest pressures near the center of 
the inlet flange with pressures increasing radially outward. 

The wall static pressures were reasonably uniform around the perimeter of the inlet 
flange. A calculation of the resultant velocity profiles from the data is presented in Figure 
B-2 at three radial locations around the circumference of the inlet flange. Significant inlet 
flow distortion is indicated entering the model, with variation between maximum and 
minimum velocities of up to approximately 45 percent. A subsequent test to assess the 
impact of the distortion on model performance was attempted employing flow 
straighteners upstream of the inlet flange. Negligible changes in the inlet profile forced by 
the inlet duct were affected however, as indicated by the inlet instrumentation 
measurements. v 




Figure^ B^3^6hows th^top, left, and right side wall static^pr66sures around/the 
circumference of t&e volute. /A/diff user section was jerdployeo in the de^n just 
downstream of the^nlet flana^to initially redupe the hjgt^ncon^Kg velocity^frpm the inlet 
jet. The resujtdnt diffusions apparent from the static pressure increase /Native to inlet 
^flange vaiues/^tatic pressures continue .to increase along the oujefwajkfo approximately 
210 degrejd^ircumfefential location, ^Xft oHhe diffuser^ection the volute was desigp^^f 
for a corfstant angu m momentum fjdwfjetd, which produces a gradual acceleratiop^df the 
flOY^Stfe^ms as s^H^eamline radius/s/educed. The ap^etdration of the flowfield ^evidenced 
pressures, is^ikewise appare^p^om the data through thp^remainder of 
the inlet scpdlK No flow/anpmalies or flow separations were observed^within the volute 
section . either by test dafaor flow visualization. Overall the flowfiejd^ppeared steady and 
reason^ly/well behaved. 
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SECTION III 
INSTRUMENTATION 


To quantify the flow phenomena observed in the Alternate Turbopump Development inlet 
flow models, a variety of probes and ports were used. The types of measurements obtained were 
as follows: 


• Static pressure 

• Total pressure 

• Static pressure ditierentials 

• Total pressure differentials 

• Flow angles. 


The locations of these probes are indicated in Figure III-l. The number of specialized 
probes implemented in the model are described below. 


A. TOTAL PRESSURE KIEL PROBES 






Total pressure kiel probes are used at three locations. These probes are located at the inlet^ ? 
Jo the model (Figure IH-2),j"at the leading edge (LE) of selected inlet guide vanes (IGVs) 
1-3), and for the dynamic test at the inducer discharge. The total pressure probes at the 
model inlet were rotated from 0 to +30 degrees, and -30 degrees to obtain an adequate number of 
points to map the flow profile exiting the Space Shuttle Main Engine LOX supply duct. During 
the dynamic test, the inducer total pressure discharge ring was also used to map the flowfield at 
three span -wise locations every 30 degrees. 


B. WEDGE PROBES 

Wedge probes are direction-sensing probes. These probes have a wedge or triangular shaped 
head. There is a pressure port on each side of the wedge near the trailing edge (TE) and a total 
pressure port at the LE. The probe can be rotated about an axis perpendicular to the wedge, 
allowing its alignment with the flow. When the pressure differential across the ports on the sides 
of the probe is zero, the probe is aligned with the flow. Wedge probes were located slightly 
forward of the IG V LE between selected vane positions and were removed when not in use. The 
wedge probes provided the inlet volute discharge angle to verify IGV angle-of-attack. 

C. WAKE RAKES 

A wake rake is a series of total pressure impact tubes, in close proximity, used to measure 
the total pressure profile at the TE of a fixed body in a flow. These rakes were used in the model 
behind one IGV on each side of the volute to examine the characteristics of the total pressure 
profile at the IGV TE. 

D. COBRA PROBES 

Like the wedge probes, the cobras were used to measure flow angle and total pressure. 
These instruments were located 90 degrees apart with three probes across the span of the IGV 
discharge annulus in close proximity to the inducer inlet plane during the static test. Although 
wedge probes are a more accurate flow angle measurement instrument, they require rotation 
capability and cannot be accommodated in this location; angle is measured using the curves 
obtained from probe calibration. Calibration of the cobra probes used in these tests was 
performed at Marshall Space Flight Center before testing. 
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SECTION V 

TEST DATA AND RESULTS 

Figure(v - l^gnows the total pressure and wall static pressure profiles measured at the inlet 
flange versus circumferential position. The instrumentation ring was clocked around the 
circumference of the inlet flange to map the incoming flowfield. Total pressure measurements 
were taken at five different radial locations, each of which displayed similar trends. To 
investigate the impact of the bypass duct and the inlet bushing (employed to isolate vibrations 
from the inlet housing) on measured inlet profiles, a subsequent test was conducted with the - 
model and bushing removed and straight discharge duct installed. The results indicated Last*. 
negligible effects on total pressures, but apparent influences from the inlet bushing on wall static y*- 

readings, as well as bypass duct influences on the side opposite the inlet model (see Figure V-2). £>0 

The decrease in total pressure on the model side (0 to 180 degrees) is due to losses and flow £ 

disturbances from the Pogo suppressor flange. Figure V-2 is considered to be most representative 
of the High-Pressure Oxidizer Turbopump (HPOTP) inlet duct static pressures. Significant inlet 
duct flow distortion entering the model is indicated by the data. Figure V-3 shows the velocity 
profiles derived from the data at varied circumferential and radial position, showing variations 
between maximum and minimum velocities of up to approximately 60 percent. 


V-4 shows the front, top, and rear volute wall staticpresSuie distributions at varied 
circumfeieTitial4ocations around both sides of the voluterKapid acceleration and diffusion is 
apparent around tKe>oh|te corners along^tbfi-^op'wall, with no indication of flow separation. The 
measured data indicates a mo^ja^ely^uniform pressure distribution along the walls on both sides 
of the volute from theipletrtothe cusp^rea^A moderate acceleration of flow along the side walls 
near the cusp^nd^deceleration along the top"wall^are apparent within the vane channels 
adjoining thtfcusp. 

Figures V-5 through V-7 show a comparison^if-the^tatic pressures on each half of the 
volute for each walL The^progles are ^nailar^nboth sides along each wall with slightly higher 
pressures measiired^n^thej^ 

Flow visualization tests revealed no flow separations or anomalies within thya^h3te~6xcept 
for a twin vortex fonfcecj^ff the inlet flow splitter that separated the^tfist-txibdel from the bypass 
duct. The vortices cairieaihtQthe fourth and fifth vanec^axmelfupi from the cusp on both sides 
of the volute where they rapi&y-^sipated^Jo^ of the vortices could be observed 

downstream of the inlet guide vanes^JQ^VsTwithin the discharge annulus. 


Figure V-8 shows t^e'^span wise total pressutae^at the leading edge (LE) of 4 of the 21 
turning vanes. The nr^ssures were reasonably uniformdccQss the span for all four vanes, with 
slightly lower values for the two vanes located at 70 degrees fr&m^the top of the volute. A wedge 
probe was inserted upstream of two vanes on each side of the volute^agd traversed span wise 
across the vane LE to measure the inlet swirl angle. Figure V-9 shows the measured swirl angles 
copa^ared to predicted values around the volute calculated from a three-din^nsionaTpotex^ti^ 
low analysis. The measured swirl angles compare favorably with the predicted values as shown. 
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pressure. The higher static pressure was the result of flow diffusing through the 
lower half of the nozzle. As will be shown later in Figure 8.19 most of the flow 
was being pushed out through the upper half of the nozzle. The pressure profiles 
were similar for all three Q/N's, but with slightly more circumferential pressure 
variation with increasing flows. At higher flows the main discharge leg includes 
more of the cutwater channel flow and more positive incidence on the splitter. 

The inner and outer walls along the branch leg were instrumented to measure the 
amount of diffusion through the branch leg. Two static pressure taps were posi- 
tioned along the outer wall and four pressure taps along the inner wall. The two 
outer wall taps were located at the entrance and exit of the branch leg. The four 
inner wall taps were located upstream of the entrance, at the entrance, part way 
through, and at the exit of the branch leg. Unfortunately, the static tap at the 
exit of the branch leg along the inner wall, where a majority of the diffusion 
would occur was plugged. Therefore the measurement of the complete diffusion 
through the branch leg along the inner wall was unobtainable. Figure 8.18 shows 
the static pressure along the inner and outer wall for three different Q/N's. A 
moderate amount of diffusion takes place along the outer wall. Once inside the 
branch leg, the inner wall shows a greater amount of diffusion as expected. The 
calculated wall loadings for design and off-design flow conditions were low. 
There were no signs of any flow separation on either wall. 

The trends for the off-design Q/N's were very similar, although there was more 
diffusion along the inner wall and slightly less diffusion along the outer wall as 
less flow traveled through the branch leg. Higher inlet Q/N's resulted in slightly 
lower flows through the branch discharge leg, due to the engine cycle flow split 
requirements. Also at higher Q/N's less of the cutwater channel flow traveled 
through the branch leg, resulting in higher positive incidence angles on the main 
- branch leg splitter. 

^ The discharge flanges of the main and branch legs were instrumented with a keil 
probe ring. Each ring contained four keil probes 90 degrees apart, to measure the 
total pressure profile exiting into the downstream ducting. Plotted in Figure 8.19 
are the total pressures versus circumferential position, for the main and branch 
leg discharge flanges at the design Q/N. Approximately 84 percent of the total 
flow through the pump passes through the main discharge leg at the design point. 
The outer wall for the main discharge flange, opposite from the splitter, was lo- 
cated at 270 degrees. The outer wall for the branch discharge flange was located 
at 60 degrees. The total pressure profile at the branch leg discharge flange was 
quite uniform. The main discharge flange showed higher pressures on the outside 
wall away from the splitter. The measured profile was consistent with test data 
of similarly designed volutes which displayed higher total pressures and flows 
pushed out toward the outer wall. The total pressures measured at the branch 
leg discharge were slightly higher than those measured at the main leg discharge. 
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